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ABSTRACT 



Wc study particle production and its density evolution and equilibration in hot dense 
medium, such as hadronic gas after quark gluon plasma hadronization and relativistic 
electron positron photon plasma. For this study we use kinetic momentum integrated 
equations for particles density evolution with Lorentz invariant reaction rates. We extend 
these equations, used before for two-to-two particles reactions (1-1-2 ^ 3-1-4), to the 
case of two-to-one and backward reactions (1-1-2 <-> 3). One type of hot dense medium, 
which we study, is hadronic gas produced at quark gluon plasma hadronization in heavy 
ions collisions in SPS, RHIC and LHC experiments. We study hadron production at quark 
gluon plasma hadronization and their evolution in thermal hadronic gas phase. We consider 
non-equilibrium hadronization model, for which the yields of the light quark hadrons are 
defined by entropy conservation. Yields of hadrons containing heavier (strange, charm, 
bottom) quarks are mainly controlled by flavor conservation. We predict yields of charm 
and bottom hadrons within this non-equilibrium statistical hadronization model. Then 
we use this non-equilibrium hadronization as the initial condition in the study of hadronic 
kinetic phase. During this time period some hadronic resonances can be produced in lighter 
hadrons fusion. This reaction is opposite to resonance decay. Production of resonances is 
dominant over decay if there is non-equilibrium excess of decay products. Within this model 
we explain apparently contradictory experimental results reported in RHIC experiments: 
S(1385) yield is enhanced while A(1520) yield is suppressed compared to the statistical 
hadronization model expectation obtained without kinetic phase. We also predict A(1232) 
enhancement. The second type of plasma medium we consider is the relativistic electron 
positron photon plasma (EP^) drop. This plasma is expected to be produced in decay of 
supercritical fleld created in ultrashort laser pulse. We study at what conditions this plasma 
drop is opaque for photons and therefore may reach thermal and chemical equilibrium. 
Further we consider muon and pion production in this plasma also as a diagnostic tool. 
Such heavy particles can be diagnostic tool to study the properties of EP'^ plasma, similar 
to the role taken by heavy hadrons production in heavy ions collisions. Finally all these 
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theoretical developments can be applied to begin a study of particles evolution in early 
universe in temperatures domain from QGP hadronization (160 MeV) to nucleosynthesis 
(0.1 MeV). The first results on pion equilibration are presented here. 
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CHAPTER 1 

PARTICLE PRODUCTION IN MATTER AT EXTREME CONDITIONS 
1.1 OUTLINE 

In this dissertation I present particle production and equilibration in different types of 
plasma medium. In this chapter [1] I overview the challenges and earlier developments, 
related to quark gluon plasma (QGP) in relativistic heavy ions collisions and electron- 
positron-photon plasma production in strong laser field that provide motivation for our 
research. In this chapter I discuss the hadronization process, in which quark gluon plasma 
breaks (freeze-out) into hadrons. The hadronization conditions have dominant influence 
on almost all final hadrons yields, even though these yields change during kinetic phase. 
For electron-positron-photon plasma we consider possibilities of its creation in strong laser 
field in future experiments. 

In chapter [2] I discuss the statistical hadronization model (SHM) in greater depth 
and show how to estimate numerically 7^ from entropy conservation an 7j i = s, c, b 
from corresponding flavor conservation strangeness conservation. In chapter [3] we con- 
sidered heavy flavor (charm, bottom) hadron production within statistical hadronization 
model. The new feature compared to the others studies is that we assume entropy and 
strangeness conservation during hadronization, accounting in this way for higher light 
quark and strangeness content in QGP. We study how this model improvement influences 
the yields of heavy flavor hadrons. We studied in depth how the (relative) yields of strange 
and non-strange charmed mesons vary with strangeness content. 

As the result of high strangeness environment we find also a relative suppression 
of the multi-heavy hadrons, except when they contain strangeness. The degree of this 
suppression depends on both, strangeness and light quarks content. When phase space 
occupancy of light and strange quark is relatively high the probability for charm quarks to 
make hadrons with strange quarks increases and probability to find the second charm quark 
among light and strange quarks decreases. These results have been published in [1,2]. 
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In chapter H] I derive equations for Lorentz invariant rates and particle density 
evolution for decay reaction l-to-2 particles and backward reaction 2-to-l (particle fusion). 
In this approach we connected the particle decay time in vacuum with kinetically modified 
particle decay time in medium, and with relaxation time for the backward reaction for 
the resonance production in two particles fusion. We calculated the relaxation times in 
medium for reactions: p ^ tt + vr; S(1385) ^ A + vr; vr*^ ^ 7 + 7. We are going to publish 
these results in [3]. 

In chapter [5] I apply equations derived in chapter H] to baryon resonance densities 
evolution in thermal hadron gas after quark gluon plasma hadronization. The goal is 
to explain ratios S(1385)/A° and A(1520)/A° reported by RHIC-experiments and also to 
predict A(1232)/A^ ratio. In this chapter I also take into account non-equilibrium condition 
at hadronization, defined by entropy conservation, used also in the model presented in 
chapter [3l I find that a significant additional yields of A(1232), S(1385) can be produced 
by the back-reaction of the over- abundance of the decay products of resonances: vr + — > 
A(1231), TT + A ^ S(1385). A more complex situation arises for a relatively narrow 
resonance such as A (1520), which can be also seen as a stable state, which is depopulated 
to increase the heavier resonance yield (A(1520) + tt ^ S*). I find that a suppression of 
the yield of such resonances, as compared to statistical hadronization model, is possible. 
The pattern of deviation of hadron resonance yields from expectations based on statistical 
hadronization model is another characteristic signature for a fast hadronization of entropy 
rich QGP. The total yield of the ground state baryons used in analysis of data (such as N, 
A) is not affected. The results are in agreement with yields of these resonances reported 
by RHIC experiments. This part of thesis is published in references [4] and [5]. 

In chapter [6] I consider e+e~7 plasma. I investigate the size and temperature limits 
of thermally and chemically equilibrated plasma drops, created by sub-optical wavelength 
high energy light-laser pulses. The plasma to become equilibrated must be opaque to elec- 
tron and photon interactions. Opaqueness condition is determined by comparing plasma 
size with the free electron and photon paths, which are calculated using thermal Lorentz 
invariant reaction rates for pair production and Compton scattering. These results are in 
preparation [6]. 

In this chapter I also study heavy particles (pion, muon) production in this plasma 
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at a temperature T ^ m^, 111^^. I argue that the observation of pions and muons can be 
a diagnostic tool in the study of the initial properties of such a plasma formed by means 
of strong laser fields. Conversely, properties of muons and pions in thermal environment 
become accessible to precise experimental study. 

In this chapter [7] I consider the pion equilibration in early universe. This chapter 
is part of reference [3]. In chapter [8] I present summary of results from each chapter, and 
conclusions. 

List of publications related to dissertation, including works in prepara- 
tion: 

1.) Chemical Equilibration Involving Decaying Particles at Finite Temperature 
Inga Kuznetsova, Takeshi Kodama, and Johann Rafelski; (in preparation) 



We study kinetic master equations for chemical reactions involving the formation 
and the natural decay of particles in a thermal hath. We consider decay channel into 
two particles and the inverse process, the fusion of two thermal particles into one. We 
derive chemical equilibrium condition for the particle density. We evaluate the thermal 
invariant rate using as input the free space (vacuum) decay time. A particularly interesting 
application of our formalism is the 7r° ^ 7 + 7 evolution in the early Universe. 



2.) Equilibration size limit of e-, e+, gamma plasma, accessible to high energy 
light pulse. 

Inga Kuznetsova, Johann Rafelski; (in preparation) 



We investigate the size and temperature limits of thermally and chemically 
equilibrated e"'"e~7 plasma drops, created by sub-optical wavelength laser pulses. For the 
plasma to be equilibrated it must be opaque to electrons and photons interactions. Opaque 
condition is determined by comparing plasma size with the free electron and photon paths, 
which are calculated using thermal Lorentz invariant reaction rates for pair production 
and Compton scattering 
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3.) Charmed hadrons from strangeness-rich QGP. 
Inga Kuznetsova and Johann Rafelski. May 2006. 6pp. 

Contributed to International Conference on Strangeness in Quark Matter (SQM2006), 
Los Angeles, California, 26-31 Mar 2006. 
J.Phys.G32:S499-S504,2006; e-Print: |hep-ph/0605307' 



4.) Heavy flavor hadrons in statistical hadronization of strangeness-rich QGP. 
Inga Kuznetsova and Johann Rafelski . Jun 2006. 18pp. 



Eur.Phys.J.C51:113-133,2007; hep-ph/0607203 



5. ) Heavy Ion Collisions at the LHC - Last Call for Predictions. N. Armesto, 
I. Kuznetsova (ed.) et al. Nov 2007. 185pp. 

Presented at Workshop on Heavy Ion Collisions at the LHC: Last Call for Predictions, 
Geneva, Switzerland, 14 May - 8 Jun 2007. 

Published in J.Phys.G35:054001,2008. e-Print: larXiv:0711.0"974l [hep-ph] 

6. ) Non-Equilibrium Heavy Flavored Hadron Yields from Chemical Equilibrium 
Strangeness-Rich QGP. Inga Kuznetsova, Johann Rafelski Jan 2008. 6pp. 

Presented at International Conference on Strangeness in Quark Matter (SQM 2007), 
Levoca, Slovakia, 24-29 Jun 2007. 
J.Phys.G35:044043,2008; larXiv:0801.0788l [hep-ph] 

In above 4 papers we study h, c quark hadronization from QGP. We obtain the 
yields of charm and bottom flavored hadrons within the statistical hadronization model. 
The important novel feature of this study is that we take into account the high strangeness 
and entropy content of QGP, conserving strangeness and entropy yields at hadronization. 



7.) Pion and muon production in e-, e-|-, gamma plasma. 
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Inga Kuznetsova, Dietrich Habs, Johann Rafelski. Mar 2008. 14pp. 
Pliys.Rev.D78:014027,2008. larXiv: 0803. 15881 [hep-ph] 

We study production and equilibration of pions and muons in relativistic electron- 
positron-photon plasma at a temperature T <^ m^, m^. We argue that the observation of 
pions and muons can be a diagnostic tool in the study of the initial properties of such a 
plasma formed by means of strong laser fields. Conversely, properties of muons and pions 
in thermal environment become accessible to precise experimental study. 



8.) Enhanced Production of Delta and Sigma(1385) Resonances. 
Inga Kuznetsova, Johann Rafelski (Arizona U. and Munich U. and Munich, Tech. U.) . 
Apr 2008. 6pp. 

Phys.Lett.B668:105-110,2008;lM2£ivi08QO352|[nucl-th] 

Yields of A(1230), S(1385) resonances produced in heavy ion collisions are 
studied within the framework of a kinetic master equation. The time evolution is driven 
by the process A ^ Nn, S(1385) ^ Avr. We obtain resonance yield both below and 
above chemical equilibrium, depending on initial hadronization condition and separation 
of kinetic and chemical freeze-out. 



9. ) Resonance Production in Heavy Ion Collisions: Suppression of A (1520) and 
Enhancement of ZI(1385) 

Inga Kuznetsova, Johann Rafelski, Nov 2008, (Arizona U. and Munich U), 13 pp. 
Phys. Rev. C 79, 014903 (2009) larXiv:0811. 14091 [nucl-th] 

10. ) Resonances Do Not Equilibrate. 
I. Kuznetsova , J. Letessier, J. Rafelski 

in Fourth Workshop on Particle Correlations and Femtoscopy (WPCF2008), Krakow, 
September 11-14, 2008, |M2£iYiQ902.2550vl [nucl-th] 
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In last two works we investigate the yield o/A(1520) resonance in heavy ion colli- 
sions within the framework of a kinetic master equation without the assumption of chemical 
equilibrium. We show that reactions such as A(1520) + vr ^ E* can favor S* production, 
thereby decreasing the A(1520) yield. Within the same approach we thus find a yield en- 
hancement for 'EllSSb) and a yield suppression for A{1520) . 

In this dissertation we consider two types of extreme matter condition: hadronic 
matter in relativistic heavy ions coUisions and e+e~7 plasma created in strong laser field. 

1.2 Relativistic Heavy Ions Collisions 
1.2.1 Overview 

The quarks differ by quantum number called fiavor. Six fiavors are known, u and d are 
light quarks, s is strange, c is charm and b is bottom quark. These quarks are arranged 
into doublets: 

The quarks on the top of doublets {u, c, t) have charge +2/3, the quarks on the bottom 
of doublets (d, s, b) have charge -1/3. Quarks u and d are the lightest. It is known that 
THu ~ 0.5md, both are in the range 2-8 MeV, and their average mass m — {ruu + Tnd}/2 ^ 
2.5 — bMeV. All stable matter is made of only u and d quarks, is between 75-125 MeV, 
rUc is about 1.3 GeV, 4.2 GeV and rut ^ 174 GeV. 

In nature quarks are always confined in hadrons. This is one of the reasons that 
mass of quarks is not well defined. Even if a lot of energy is applied to separate qq pair, 
at some point of separation it becomes energetically preferred to create one more confined 
qq pair. Therefore in the result we have two confined quark-antiquark pairs. 

In relativistic heavy ions collisions hot and dense fireball of nuclear matter is 
created. We believe that at very high energy the deconfinement of partons is expected and 
results to new phase of matter, quark gluon plasma (QGP). The temperature of the fireball 
drops with expansion. At critical hadronization temperature the hadrons are formed again. 
Observing final yields of hadrons we try to learn about physics of QGP formed in heavy 
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ions collisions. The difference in hadrons yields between pp collisions and heavy nuclear 
collisions can justify the existence of QGP state. [8] 

In part of this dissertation we develop models that explain (or predict) yields 
of some baryon resonances and heavy flavor hadrons, reported in heavy ions collision 
experiments. We use experimental data reported in Relativistic Hadronic Collider (RHIC) 
experiments and predict yield of heavy flavor hadrons at Large Hadronic Collider (LHC) 
energies. RHIC is built in Brookhaven National Laboratory to create and search the new 
phase of matter (QGP), colliding Au ions with energy in center of mass frame ranging 
from a/s — 20 to 200 AGeV. The much higher energies are expected at the Large Hadronic 
CoUider (LHC), which began to operate recently at CERN. Here we will also refer to Super 
Photon Synchrotron (SPS) accelerator with fixed target results at CERN, where energies 
of accelerated Pb ions up to 158 ^GeV were achieved. 

In our approach we assume that the following evolution stages are present in 
heavy-ions collisions: 

1. Primary partons collide, practically all heavy c, b quarks are produced; 

2. A thermalized parton state within r — Tth — 0.25 — 1 fm/c is formed. By the end of 
this stage nearly all entropy is produced. 

3. The subsequent chemical equilibration: diverse thermal particle production reactions 
occur, allowing first the approach to chemical equilibrium by gluons g and light non- 
strange quarks q — u,d. 

4. The strangeness chemical equilibration within r ~ 5 fm/c). 

5. Chemical freeze-out (hadronization) near r ~ 10 fm/c), when quarks become confined 
in hadrons. 

6. Kinetic phase, where hadrons can scatter and be regenerated. 

7. Kinetic freeze-out. 
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1.2.2 Statistical hadronization model 

The transition from quark gluon plasma to hadrons can be described by statistical 
hadronization model. SHM arises from the Fermi multi-particle production model [9]. In 
Fermi model it was micro-canonical ensemble the flavor and energy are conserved exactly. 
It was developed further by Landau [10] and Hagedorn [11]. In [11] the infinite number 
of hadrons with increasing mass results to exponential mass spectrum, which diverges at 
critical temperature. This was the first evidence toward the phase transition from confined 
hadrons to QGP. There was a transition towards the finite size of hadrons, composed of 
quarks, and phase transition between QGP and hadronic gas in Hagedorn and Rafelski 
work [12,13]. 

The transition from micro-canonical to canonical, and grand-canonical ensembles, 
where averaged flavor is conserved, simplifies the computational effort considerably [14]. 
This important step does not in our context introduce the hadron phase, although before 
the understanding of QGP this of course was the reaction picture: a highly compressed 
hadron gas matter evaporates particles. Today, it is the highly compressed hot quark-gluon 
matter that evaporates particles. 

Yields of most hadrons arc described by statistical hadronization model without 
kinetic phase. For resonances with strong decay there are the deviations of their yields 
ratios reported in experiment from predicted by statistical hadronization only. This makes 
necessary to introduce the kinetic phase, where yields of hadrons can be changed by reac- 
tions between hadrons. 

The SHM is related to assumption that system also described by hydrodynamic 
model. This includes assumption that for strong interactions relaxation times are small 
enough that prchadronic matter (QGP) is in thermal equilibrium and the system can be 
treated relativistic fluid. 

In a hydrodynamic description with flow in the longitudinal and transverse direc- 
tions the current of particle with mass m is 

-^J^. j d'pnv" - m^)^(poK/(p) = / ^P"/(P> 7); (1-1) 

p is 4- vector momentum of particle: 

p^^{E,p) (1.2) 
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In this equation we use Lorentz invariant distribution functions for bosons and 
fermions: 

where the Fermi '(+)' and Bose '(— )' distributions are indicated, g is the degeneracy 
factor, T is the temperature, is the frame four-vector, which in the observer rest frame 
of is 

u^^(l,^). (1.4) 

Then number of particles within the element of 3-dimensional freeze-out surface 
in the Minkovski space-time da^ — {(Px, 0) is defined by 

dN{a) = fda^ (1.5) 

The distribution of particles at this surface is given by the Cooper- Fryer formula [15]: 

In Boltzmann limit, which is good approximation for considered temperatures for most 
hadrons, we can omit '1' in determinator of distribution function: 

/(p, x) ^ 7Ae-"-^/^. (1.7) 

The important parameters of the SHM, which control the relative yields of parti- 
cles, are the particle specific fugacity factor A and space occupancy factor 7. The fugacity is 
related to chemical potential ^ — Tin A. The occupancy 7 is, nearly, the ratio of produced 
particles to the number of particle expected in chemical equilibrium (7 = 1 is chemical 
equilibrium) . 

The fugacity A is associated with a conserved quantum number, such as net-baryon 
number, net-strangcncss, heavy flavor. Thus antiparticlcs have inverse value of A, and A 
evolution during the reaction process is related to the changes in densities due to dynamics 
such as expansion. Here we always consider particle anti-particle symmetry or A = 1. This 
condition is almost satisfied for RHIC and LHC energies. 7 is the same for particles and 
antiparticlcs. Its value changes as a function of time even if the system does not expand, 
it describes buildup of the particular particle species. For this reason 7 is changes rapidly 
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during the reaction, while A is more constant. Thus it is 7 which carries the information 
about the time history of the reaction and the precise condition of particle production 
referred to as chemical freeze-out. 

These distribution functions with phase space occupancy 7 can describe quark 
in QGP and also hadrons multiplicity after hadronization. For phase space occupancy for 
hadron is a product of 7s for each quark in this hadron. For example for kaon K = Iqls, 
for pion 7^ = 7^, etc. 

Then we have, as expected for Boltzmann limit from Eq. (11.61) and (11. 7p : 

Integrating this equation over momentum, the number of particles of type 'i' with 
mass rrii per unit of rapidity is in our approach given by: 

= 7,<V. (1.9) 

Here dV is system the volume of a fireball, and n^'^ is a Boltzmann particle density in 
chemical equilibrium: 



n? = 9^j ^A,exp(-^p2 + ^2/y) ^ x,^^g.W{mjT), (1.10) 



and 



W{x) = x^K2{x) ^ 2 for X ^ 0. (1.11) 

Both, iTLiC^ — > mj, and kT —>■ T, are measured in energy units when h,c,k ^ 1. 
In non-relativistic Boltzmann limit Eq. Ol.lUp can be expand as: 



= 2^^^V2r^'"P(-"^^/^)('+8;;^ + T28UJ ^ ^ ^ 

Often we can use the first term alone for heavy flavor hadrons, since T/m « 1, however 
the asymptotic series in Eq. (I1.12p converges slowly (if at all) and one should proceed with 
caution. 

We use occupancy factors 7^*°^ and 7^^ for QGP and hadronic gas phase respectively, 
tracking every quark flavor {i = q, s, b, c) . We assume that in the QGP phase the light 
quarks and gluons are adjusting fast to the ambient conditions, and thus are in chemical 
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equilibrium with '-)^q 1. For heavy, and strange flavor, the value of 7*°^ at hadronization 
condition is given by the number of particles present, generated by prior kinetic processes, 
see Eq. (fLQjl . 

The yields of different quark flavors originate in different physical processes, such 
as production in initial collisions for c, 6, s, and for s also production in thermal plasma 
processes. In general we thus cannot expect that 7^^ will be near unity at hadronization. 
However, the thermal strangeness production process GG — > ss can nearly chemically 
equilibrate strangeness flavor in plasma formed at RHIC and/or LHC [26], and we will 
always consider, among other cases the limit 7^ ^ 1 prior to hadronization. 

Differing from other recent studies which assume that the hadron yields after 
hadronization are in chemical equilibrium at list in respect of light quarks or light and 
strange quarks together (7g(s) = 1) [27,28], to evaluate yields of final state hadrons we 
enforce conservation of entropy (determine 7,), and the flavor s, c, 6 quark pair number 
(determine corresponded 7^) during phase transition or transformation. 

The faster the transition, the less likely is that there is significant change in strange 
quark pair yield. Similarly, any entropy production is minimized when the entropy rich 
QGP breakup into the entropy poor HG occurs rapidly. The entropy conservation con- 
straint fixes the final light quark yield. We assume a fast transition between QGP and HG 
phases, such that all hadron yields are at the same physical conditions as in QGP breakup. 

Assuming that in the hadronization process the number of 6, c, s quark pairs does 
not change, the three unknown 7^, 7^, 7^^ can be determine from their values in the QGP 
phase, 7^, 7^, 7^ (or A^^*^) and the three flavor conservation equations, 

Nf = N^ = dNi, i = s,c, b. (1.13) 

In order to conserve entropy: 

_ _ 5, (1.14) 

a value 7^ > 1 is required for T > 180 MeV when in the QGP phase •j^q = 1. 7^ > 1 
is needed to compensate entropy lost because of decreasing number of degrees of freedom 
in hadronic gas and increasing mass of hadrons. This implies that yields of hadrons with 
light quark content are, in general, not in chemical equilibrium, unless there is some ex- 
traordinary circumstance allowing a prolonged period of time in which hadron reactions 
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can occur after hadronization. Chemical non-equilibrium thus will influence the yields of 
heavy flavored particles in final state as we shall discuss in this work. 

As noted at the beginning of this section, the use of the hadron phase space 
(denoted by H above) does not imply the presence of a real physical 'hadron matter' 
phase: the SHM particle yields will be attained solely on the basis of availability of this 
phase space as noted at the beginning of this section. Another way to argue this is to 
imagine a pot of quark matter with hadrons evaporating . Which kind of hadron emerges 
and at which momentum is entirely determined by the access to the phase space, and there 
are only free-streaming particles in the final state. 

Thinking in these terms, one can imagine that especially for heavy quark hadrons 
some particles are pre-formed in the deconfined plasma, and thus the heavy hadron yields 
may be based on a value of temperature which is higher than the global value expected 
for other hadrons. For this reason we will study in this work a range 140 < T < 260 MeV 
and also consider sensitivity to this type of two-temperature chemical freeze-out of certain 
heavy hadron yield ratios. 

We use this non-equilibrium hadronization model to predict heavy flavor hadrons 
yields for RHIC and LHC conditions. 

Here we also use grand canonical hadronization model for heavy flavor hadroniza- 
tion. The grand-canonical hadronization condition is conservation of symmetry of strange 
and antistrange hadrons in reactions: 

< Ui > - < nj >^ 0. (1.15) 

When number of flavor is small this condition must be replaced by sharper condition: 

KUi-n-i >= 0. (1-16) 

This more exact condition results to canonical suppression of given flavor yield. There 
is the large effect when the quarks multiplicity for this flavor is about unity and effect 
increase for multiflavor hadrons. This shows that the smaller multiplicity of given flavor (s) 
quarks is, it is more difficult for them to find each other to bind to hadron. Moreover here 
for hadrons with one heavy quark we normalize yield by quark multiplicity, and canonical 
suppression is the same for both yields and canceled. 
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For LHC even multiplicity of c quarks per rapidity is large dNc/dy ~ 10 and for 
hidden charm (cc) mesons we almost do not have canonical suppression. There is small 
canonical effect for hidden charm mesons at RIHC, which we did not consider. This effect 
may be large for B^, bb mesons. 

1.2.3 Strangeness enhancement 

Strangeness enhancement in relativistic heavy ions collisions, compared to pp collisions, was 
proposed to be an indication of presence of deconfined state by Rafelski and Muler [17-19]. 
In [16, 17] it was shown that strangeness is mostly produced in QGP phase in reactions 

q + q^s + s, (1.17) 
g + g^s + s. (1.18) 

Strangeness production in the reaction (11.181) is dominant in qgp. Strangeness can also be 
produced in hadronic gas in reactions 

TT + N K + Y; (1.19) 

TT + n ^ K + K, (1.20) 

where Y is strange baryon (hyperon). However the strangeness production rate in thermal 
QGP the strangeness production rate is an order of magnitude larger than in hadronic 
gas [16, 17]. 

Therefore strangeness is mostly produced in QGP phase. Before hadronization in 
QGP strangeness multiplicity is near equilibrium. Then in our model during fast hadroniza- 
tion number of strange quarks does not have time to change. As the result after hadroniza- 
tion the multiplicity of strange hadrons is higher than for models which just put '-/^'^ = 1 
the same as in QGP. The explanation is that in hadronic gas strange hadrons are more 
massive than free strange quark in QGP. Therefore distribution function with the same 7^ 
multiplicity of strangeness may be smaller than in QGP. This effect depends on hadroniza- 
tion temperature. In order to conserve strangeness the j^'^ > 1 may be needed. We will 
add quantative explanation in chapter [2. 3[ 

Strangeness content is defined by strangeness to entropy ratio s/S, which can be 
measured. s/S increases with energy of collision increase. The explanation is that the hot 
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state where the energy of gluons exceeds the threshold for strangeness production hves 
longer. For RHIC at central collisions s/S* ~ 0.03, for LHC it is expected to be about 0.04. 

s/S increases for higher energy of collision. Possible explanation for higher energies 
the hot state with temperature above threshold for strangeness production lives longer. s/S 
increases with centrality of collision or number of participants Npart [20] 

Also experimentally strangeness content can be found by measuring ratios of yields 
of hadrons with different strangeness content, for example K/tt or A/p which are oc 'js/lq- 
s/S increases with 'js/lq grow. The numerical results for connection between ratio s/S 
and 'js/lq are shown in chapter [3] 

In figure [LT] we show experimental results for ground states of strange baryons A, 
S, Q yields, normalized to the yields in pp or pBe collisions E{i) as a function of Npat for 
Au+Au and Pb+Pb and pp collisions. For each species, i, the yield enhancement, E{i), 
above that expected from Npart scaling was calculated using: 



Y^eld^JN^ 



The inclusive proton data illustrate the effects for non-strange baryons. The figure is taken 
from [21]. 

The difference between particles and antiparticles is due chemical potential, which 
is higher for lower energy collisions. We see increases in the yield of strangeness containing 
hadrons with Npart as predicted in [20]. For H yield (oc 7^) the increase is larger than for 
A (oc 7s) and the largest effect is for VL yield (oc 7^). 

1.2.4 J/\& suppression 

The observation of heavy hadrons containing more than one heavy quark, for example 
charmonium can give information about deconfined QGP phase (Matsui and Satz 1986). 
In deconfined phase charm and anticharm quarks from different nucleon-nucleon collisions 
diffuse, meet and can produce charmonium. The probability of charmonium production 
and therefore the enhancement of observed yield is expected as compared to single nucleon- 
nucleon collisions without deconfined phase [25,27,33]. Similar enhancement is possible 
for Br mesons. 
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Figure 1.1: (color online) E[i) as a function of Npart for A, A, S+, S^, Q~ and 
inclusive p. Boxes at unity show statistical and systematical uncertainties combined in 
the p+p (p+Be) data. Error bars on the data points represent those from the heavy-ions. 
The solid markers are for Au+Au at ^/s=200 GeV and the open symbols for Pb+Pb 
at y^=17.3 GeV. The arrows on the right axes mark the predictions from a GC (grand 
canonical) formalism model when varying T from 165 MeV (-E(S~)=10.7, ii^(A)=2.6) to 
170 MeV (i?(S~)=7.5, E(A)=2.2) [22]. The red arrows indicate the predictions for S and 
the black arrows those for A. 

However in SPS (NA50 experiment) and RHIC experiments otherwise the J/\E' 
suppression is reported [34,35]. This suppression can not be explained by standard nuclear 
absorption model. In NA50 experiment J/\E' yield is reconstructed through observation of 
decay J/"^ fifi only. Then nuclear absorption model with in medium mass modification 
was proposed. Small shift of J/\E' mass makes decay J/'^ —>■ DD also possible. However 
in medium mass modification has not been observed. 

Our present work suggests that it is important to account for the binding of heavy 
flavor with strangeness, an effect which depletes the eligible supply of heavy flavor quarks 
which could form ) J/\E'(cc) and Bc(6c, 6c) [2]. 

The strangeness enhancement and J/\& suppression are considered as the two most 
profound evidence of QGP state. There are also indications that after hadrons form there 
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is a kinetic phase where hadrons can rescatter. Heavier hadrons do not decay only but 
they can also be recreated by lighter particles. After kinetic freeze out hadrons decay only 
and their yields are reconstructed by observation of their decay products. 

1.3 Kinetic theory particles production and equilibration in QGP and thermal 
hadronic gas. 

1.3.1 Boltzman equation 

Relativistic Boltzmann equation is 



1 , 9 
—Pt: — 



where fi{p, x) are particles distribution functions, the index i denotes the type of particle, 
is its 4-momentum, F^^ is external force and rj^ and xt are the emission and the absorption 
coefficients for the production of a particle of type "i" via the physical process labeled by 
q. The collision term presentation, right hand side is taken from [38]. Here we will consider 
cases with F'^ = 0. For two-to-two particles reaction: 

1 + 2^3 + 4 (1.22) 

the collision integral for absorption of particle '3' can be written (similar to [38]): 

xVs = y cip4cipicip2W^p3,P4;pi,P2/4(P4,a;)/3(p3,a;), (1-23) 

where W^p3,p4;pi,p2 is transition function, connected to transition probability w;p3,p4;pi,p2 

V = 

W^P3,P4;pi,P2- Similar for particle '3' emission coefficient: 

'7*'' = y C?P4C?PlC?P2W^pi,p2;P3,P4/l(P4,a;)/2(p3,x), (1.24) 

Here we took dp = d^p6{pf — mf)9{p^^) = d^p/{2E). If we assume that in reaction fll.22p 
particles obey Fermi or Bose statistic the probability of transition 34 12 is 

«^P3,P4;Pi,P2 = 77T-^{P3P4 \M\pip2fS'^ {pi + P2 - p-i - Pa) (1 ± f{x,Pi)) (1 ± f{x,P2)) ■ 

(1.25) 
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Similar we can write a probability for particle '3' emission: 

Wpi,p2;p3,p4 = -^^^(^1^2 \M\pzPa)^5^ [pi + P2 - P3 - Pa) (1 ± f{x,P3)) (1 ± f{x,PA)) . 

(1.26) 

The factors (1 ± f{x,p)) are Bose enhancement or Fermi suppression factors, which shows 
that we can have more than one or only one particle in given final state. The transition 
probability is proportional to the matrix element {piP2 \M\ p^Pi)"^-, which we took as known 
from others studies. 

1.3.2 Equation for particle density evolution and reaction rates 

In this dissertation we investigate resonances evolution during kinetic phase (or ther- 
mal hadronic gas). The kinetic equations, we use, are similar to those used before for 
strangeness production and equilibration in QGP [16,17,39-41]. These equations can be 
obtained by integration of Boltzman equation 11.221 over momentum of studied particle i. 
We can do it if we know momentum dependence of distribution function. It is possible in 
case if thermal equilibrium establishes faster by scattering reactions than chemical equi- 
librium of particle i. Then to study particle '3' yield evolution in reaction 11.221 we can 
assume thermal distribution for all particles: 

Here historically we changed A7 to T, which we call fugacity. 

We assume that reaction (11.221) does not change momentum distribution, the T 
changes. Temperature can change with volume expansion. In our case from entropy 
conservation T^V ~ constant. 

For example, strangeness is produced mostly in thermal gluons fusion reac- 
tion (ll.lSp . Backward reaction also becomes important near equilibrium point. These 
reactions studied before in approach considered here are two-to-two particles reactions [11221 

In example of strangeness production these reactions are bosons fusion into 
fermion-antifermion pair and backward: 



6 + 6^/ + /. 



(1.28) 
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The others possible reactions, which we will consider here are 

f + f--f' + f'; (1.29) 

b + f{f)--b + f{f). (1.30) 

Under assumptions described above we integrate Boltzmann equation 11.221 for 
particle 3(or 4) over momentum phase space dpi = d^pi6{pf — '>TLj)6{p^) we obtain the 
evolution equation for current of produced particle 3 or 4 in reactions fll.22p : 

^^^12^34 ^^^"34 ^12 /i oi\ 

~ dVdt dVdt ' ^ ■ ^ 

The current of produced particle is defined by Eq.f ll.ll) . The covariant derivative is 

j[', = ^d,{^jn. (1.32) 

The Lorentz invariant particle production and annihilation rates are 
dWi2^34 T1T2 9192 f d^pi , , . f d^p2 f , . f d^p3 f d^ps 4 



dVdt l + /(27r)i2_/ 2Ei '^^'^ ' J 2E2 ' J 2E3 J 2E4 

J2\{piP2 \M\p3P,)f5' (p,+p,-p,-p,) (1 ± hj{T3,P3)) (1 ± /,j(T4,P4))(.1.33) 



spin 



dW34^i2 T3T4 5(35(4 f d^ps , , . [ d^pi f / ^ [ d?Vz f d^Ps ,^ .4 



dVdt 2 (27r)i2 7 2E3 ''^'^-"'J 2E4 ' J 2E3 J 2Ei 

5^ KP3P4 \M\piP2)f5' (p,+p2-p3-p,) (1 ± /6,/(Ti,Pi)) (1 ± /5j(T2,P2))(1.34) 

spin 

Using the relation 

1 T ff,t = Tr^e^P^/^hj, (1.35) 

and time reversal invariance of matrix element {{psPa \M\pip2)f = \{piP2 \M\p3P4)f , we 
obtain equation connecting fermions pair production and annihilation rates: 

1 rfl^l2^34 _ 1 dW3i-.i2 _ „ 

dVdt - T3T4 dVdt - ^^-^^^ 

From this equation we see that fermions production and annihilation rates are equal and 
reactions are in equilibrium when Ti, = T f = 1. Local equilibrium is also possible when 
T5 = Tj 7^ 1. However in the system with many reactions it may be impossible to satisfy 
equilibrium conditions for all reactions when Tj 7^ 1. Also in QGP the gluons are likely at 
their highest possible density with T„ = 1. 
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Using Eq. (11.361) . we can rewrite Eq. fll.3ip as 

J^'; = (T1T2 - T3T4)i?12^34. (1.37) 

For homogeneous expansion this equation can be written as 

(T1T2 - T3T4)i?l2^34. (1.38) 



1 (iiV3{4) 



V dt 

We can rewrite this equation as equation for T3(4): 

+ -^T + n,^,^- - (T1T2 - T3T4)i?i2^34. (1.39) 
We introduce reaction relaxation time 

ri2 34 = ^^3(4)/^T3(4) 

2a/(TiT2)-Ri2^34 

This time is on the order of magnitude of time needed to reach equihbrium condition for 
fugacities: 

r^eg^eg = T^-^Tf . (1.41) 

In simpler example V =const, T =const, Ti = T2 = T and T3 = T4 = T' Eq. (11.391) is 

1 dV ( T'2\ 1 

which has for T,(t = 0) = the simple analytical solution [17]: 

r = Ttanh(t/2ri2^34). (1.43) 

For t 00, near to chemical equilibrium, T'/T ^ 1 — e^*/'^i2«34^ while for t — > 0, at the 
onset of particle production with small T' we have T' = t/(2r'). Hence, near to chemical 
equilibrium it is appropriate to use factor 2 in definition of relaxation time Eq. (ll.40l) . 

For Boltzmann distribution, for example for strange particles production in 
hadronic gas (reactions (11.191) . (11.201) ). we may use reaction cross section in center mass 
frame a to estimate production rates in two body processes and the relation. Then Lorentz 
invariant reaction rate is [8,39]: 

-Ri2^34 = (crvrei) nin2, (1.44) 
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where Vrei is relative velocity of particle 1 in respect of particles 2. The cross section is 
connected to matrix element by 

(TV12E1E2 = j d^P3S{pl-m3)9{pl)d'^p^6{pl-m4)9{pl)6'^{pi+P2-p3-P4){piP2 iMlp^p^)^ 

(1.45) 

The rate can be evaluated, using 

V12E1E2 = X2{s) = {s- (mi + m2)^){s - (mi - ma)'), (1.46) 

as [8,39], where y/s = {E1 + E2) is total energy of interacting particles 1 and 2 in the center 
mass frame, mi and m2 are masses of the initial interacting particles 

327r4 1 + / J^^^ ^ 

(compared to reference [8] our definition is changed -Ri2-*34 -Ri2-*34/('^i'^2)) where gi, 
g2 and Ti, T2 are degeneracies and fugacities of the initial interacting particles. 

1.3.3 Resonance production during kinetic phase 

Resonances are very short lived hadrons, baryons and mesons, with width (inverse lifespan) 
on the order of 1 - 100 MeV. Because of their very short lifespan yields of resonances can 
not be observed directly. Their yields are reconstructed, using their decay invariant mass 
method. 

Some of resonance to similar non-resonance hadron yields ratios reported by RHIC 
and SPS experiments show deviations from those observed in pp collisions and from pre- 
dicted by statistical hadronization model along. In figure 11.3.31 we show experimental 
resonance to non-resonance ratios as a functions of centrality (number of participants) in 
Au-Au collisions reported in [43] (preliminary) on the left and in [44] on the right. 

From figure 11.3.31 we see that some of resonance to non-resonance ratios shows 
noticeable dependence on collision centrality (enhancement or suppression). The possible 
explanation is that reactions in kinetic phase of hadronic gas can influence. The idea that 
resonance can be regenerated in kinetic phase was pointed out within ultra relativistic 
quantum molecular dynamics (UrQMD) model [45]. 

In chapter 4 we present model, which explains the S(1385) and A(1230) enhance- 
ment and A(1520) suppression compared to statistical hadronization without kinetic phase 
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Figure 1.2: Left [43]: The K*^/K, p^/n, A++/p and (p/K ratios as a function of number of 
charged xhadrons in p+p (open symbols) and various centraHties in Au+Au (sohd symbols) 
collisions. Right [44]: Resonance to stable particle ratios for p + p and Au + Au collisions. The 
ratios are normalized to unity in p + p and compared to thermal and UrQMD model predictions 
for central Au + Au. 



and pp collisions in good agreement with experimental results in figure . In our approach 
we assume thermally equilibrated hadronic gas in kinetic phase. We assume that hadronic 
gas temperature changes together with volume {T^V ~ const from entropy conservation). 
We consider only relevant 

1 + 2^3; (1.48) 

reactions, estimating resonance yield. Their rates in most cases are much faster then for 
2-to-2 particles reactions. Threshold energy for resonance production is smaller. It is 
possible that 2-to-2 particles reactions also have influence in some cases. Non-equilibrium 
can accelerate these reactions (rate is proportional to corresponding Ts). This question 
we leave for future research. Then we use equations similar to equations from sections 11.31 
extended to the case of reactions fll.48p for detailed balance and particles evolution. We 
will derive these equations in chapter HI The new feature in our detailed balance for 
reactions 11.481 is that we connect Lorentz invariant rate for reactions in both direction in 
thermal medium with resonance decay time in vacuum, which is known. This way we do 
not need to know cross section or matrix element for resonance production. 
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1.3.4 Entropy in QGP fireball 

When heavy nuclear collide at high energy we expect that QGP is formed. The QGP 
consist of quarks (fermions) and gluons (massless bosons). Number of degrees of freedom 
in this phase is 

7 

9 = 9a + l9g, (1-49) 

where gluon degeneracy gg = 2s{N^ — 1) = 16 (2^ is spin degeneracy, = 3 is the 
number of colors) and quark degeneracy Qq = 2sNcnf = Qrif, nj is number of light flavor 
{rrif < T). If semi-massive strange quarks are present nj = 2.5. Factor 7/4 in Eq. fll.49p 
shows the presence of particles and antiparticles (factor 2) and the smaller fermion phase 
space compared to bosons, defined by exclusion principal. 

The entropy content is seen in the final state multiplicity of particles produced 
after hadronization. More specifically, there is a relation between entropy and particles 
multiplicities, once we note that the entropy per particle in a gas is: 

w w = ! = 

for massless Bose, classical (Boltzmann) and Fermi gases, respectively. Effectively, for 
QGP with u, s, d, G degrees of freedom, S^/N^ ~ 4 is applicable for large range of masses. 
Thus: 

dS^^AdN'^. (1.51) 

This in turn means that final state particle multiplicity provides us with information about 
the primary entropy content generated in the initial state of the QGP phase. 

It is today generally believed that there is entropy conserving hydrodynamic ex- 
pansion of the QGP liquid. Entropy is conserved in the fireball, and the conservation of 
entropy density a flow is expressed by: 

4^7^ =0. (1.52) 

where u'^ is local four velocity vector. A special case of interest is the Bj0rken scenario [23]. 
In this scenario, we assume that 

1. the energy of the colliding particles is so large that the flow of energy and matter after 
heavy ions collision remains unidirectional along the original collision axis; 
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2. the transverse extend of the system is so large that the existence of the edge of mater 
in a direction transverse to the coUisional axis is of a httle relevance. 

This scenario suggests that the natural variables for the dynamics of rapidly ex- 
panding in longitudinal direction flow are proper time r(t, z) and rapidity y{t, z): 

r = {e-z-YI\ ^/=^ln(^), (1.53) 

where z is coordinate in longitudinal direction. In this case Eq. (11.521) can be solved exactly 
assuming as initial condition scaling of the physical properties as a function of rapidity. This 
implies that there is no preferred frame of reference, a situation expected in very high energy 
collisions. Even if highly idealized, this simple reaction picture allows a good estimate of 
many physical features. Of relevance here is that the exact solution of hydrodynamics in 
(1+1) dimensions implies 

^ = Const. . (1.54) 
dy 

Thus entropy 5* is not only conserved globally in the hydrodynamic expansion, 
but also per unit of rapidity. Though we have (1+3) expansion, Eq. (11.541) holds as long 
as there is, in rapidity, a flat plateau of particles yields. Namely, each of the domains 
of rapidity is equivalent, excluding the projectile-target domains. However, at RHIC and 
LHC energies these are causally disconnected from the central rapidity bin, where we 
study the evolution of heavy flavor. The entropy we observe in the final hadron state 
has been to a large extent produced after the heavy flavor is produced, during the initial 
parton thermalization phase, but before strangeness has been produced. In order to model 
production of hadrons for different chemical freeze-out scenarios of the same reaction, 
we need to relate the entropy content, temperature and volume of the QGP domain. 
We consider for a m, d, G-chemically equilibrated QGP, and allowing for partial chemical 
equilibration of strangeness, the entropy content. 

The entropy density a can be obtained from the equation 

where the thermodynamic potential is: 

Fq(T,A„1^) = -TlnZ(QT,A„\/)Q. (1.56) 
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Inside QGP the partition function is a product of partition function of gluons Zg, light 
quarks Zq and strange quarks Z^, hence: 

lnZ = lnZg + lnZq + lnZ,; (1.57) 

where for massless particles with \q = 1 

InZ, = ^yT^ (1.58) 

lnZ„ = l^-^VT\ (1.59) 
^ 4 90 ^ ^ 

Here we take into account that the number of degrees of freedom of quarks and gluons is 
influenced by strongly interactions, characterized by strong coupling constant a^'- 

9g = 2,8, (^l-^a, + ..}j ; (1.60) 

gq = 2,3,2; + . (1.61) 

The case of strange quarks is somewhat more complicated, since we have to con- 
sider the mass, the degree of chemical equilibration, and guess-estimate the strength of 
QCD perturbative interactions. We have in Boltzmann approximation: 

InZ, = 2p/a^VT^ (1.62) 
gs = 2s3c7?0.5W{ms/T)(l-k^y (1.63) 

W{x) = x'^K2{x), where K2{x) is Bessel function of the 2nd order. We allow both for 
strange and antistrange quarks, factor 2p/a (which is for massless fermions 2-7/8 = 7/4). k 
at this point is a temperature dependent parameter. Even in the lowest order perturbation 
theory it has not been evaluated for massive quarks at finite temperature. We know that 
for massless quarks A; ~ 2. Considering expansion in m/T, for large masses the correction 
reverses sign [29], which result supports the reduction in value of k for m c:^ T. We will 
use here the value k = 1 [26]. 

The entropy density following from Eq. (11.551) is: 

47]-2 7 4 A 

^ = -^(9, + p,)T' + :^2,/.,g,T' + -. (1.64) 
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For strange quarks in the second term in Eq. fll.64p we set the entropy per strange quarks 
to 4 units. In choosing Sg/Ns = 4 irrespective of the effect of interaction and mass value 
iris/T we are minimizing the influence of unknown QCD interaction effect. 

The last term in Eq. fll.64p comes from differentiation of the strong coupling con- 
stant in the partition function with respect to T, see Eq. fll.55p . Up to two loops in the 
/3-function of the renormalization group the correction term is [30]: 



A = {boal + bial) 



27r^4 %57r ^4 
3 18 



:i.65) 



with Uf being the number of active fermions in the quark loop, "n,/ ^ 2.5, and 



For the strong couphng constant as we use 

°-<^'" l+CMr/r.) - C- = 0.760 ±0.002. (1.67) 

where Tc = 0.16 GeV [8]. This expression arises from running of the energy scale 

at = 2ttT, and the value as{Mz) = 0.118 (M^ =91.19 GeV). A much more sophisticated 
study of the entropy in the QGP phase is possible [31], what we use here is an effective 
model which agrees with the lattice data [32]. 

1.4 Electron-positron-photon plasma 

1.4.1 Electron-positron-photon generation and equilibration 

The other part of dissertation studies the possibility of e"'"e~7 plasma production and 
equilibration in strong laser field and also heavy (muon, pion) particles production and 
equilibration in this plasma. We also included here subsection about early universe. In 
the laboratory these heavy particles can be diagnostic tool for properties of e~^e~'~f plasma, 
similar as in case of heavy ions collisions. 

The electron-positron plasma can be found in many astrophysical objects as active 
galactic nuclear, pulsars, gamma ray bursts. Over last fifteen years the huge progress in 
laser intensity was achieved. The formation of the relativistic, electron-positron-photon 
e~,e"'",7 plasma (EP^, temperature T in MeV range) in the laboratory using ultra-short 
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pulse lasers is one of the current topics of interest and forthcoming experimental effort [48, 
49]. 

One of questions we study if this electron-positron plasma can be in thermal and 
chemical equilibrium with photons. Here we show that to have this equilibrated opaque 
plasma it is necessary to focus pulse energy in small size 10-1 nm. This follows from 
opaqueness condition. On the other hand in this case of small size less pulse energy is 
needed to create large electric field, close to breakdown Schwinger field Eg. 

The Schwinger field is a field necessary for virtual electron-positron pair to gain the 
energy 2meC^ during the time 6t defined by Heisenberg uncertainty principal 6t = h/meC^- 
The energy gain length is c6t = \b, where A;, is de Broil wave length. Then the Schwinger 
field is 

= = lO^^V/cm. (1.68) 

The laser field is connected to laser intensity Ii as E"^ = ZqIi, where Zq is vac- 
uum impedance, Zo = 377^2. We find that to create Eg intensity Is = lO^^W/cm"^ is 
necessary [48]. 

In [50] the time scale necessary to convert laser field energy into e'^e^ pairs energy 
was evaluated to be in oder of 1 — 10^^ fs for the corresponded field E ^ (0.4 1.0)Es. 
Pulse duration has to correspond to this field to plasma conversion time. This time and 
field ranges may be enough to produce desirable high density e~^e~ plasma. 

These physical conditions should become accessible in the foreseeable future upon 
the development of wavelength compression technology employing an optical wavelength 
laser beam reflected from a relativistic mirror, generated by a pulsed high intensity 
laser [51]. In [51] the thin plasma slap plays mirror role accelerating in the radiation 
pressure dominant regime. The flying mirror reflects counter-propagating radiation caus- 
ing its frequency multiplication by squared Lorentz factor of mirror because of double 
dopier effect. The scheme of radiation reflection from accelerated double-surface mirror is 
shown in figure 11.31 

1.4.2 Pion and muon production in e+e 7 plasma 

We also study the production of heavy particles in e"''e~7 plasma. The purpose of this 
research is to use observation of heavy particles (pion, muon) yields as a tool for study of 
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Driver 



reflected 




Source 



Figure 1.3: [51] The scheme of the double-surface mirror. The uhra-intense driver going 
from the left accelerates the mirror in the radiation pressure dominant regime. In its turn, 
the mirror reflects the intense source sent from the right. 



properties, similar to what we do in the case of heavy ions collisions. Also it may be useful 
for study of the heavy particles production reactions itself. 

The 7r°, vr^, /i^ can be produced in e"'"e~7 plasma. For T << myro (starting 
from temperature about few MeV) the neutral pions are most effectively produced in two 
photons fusion: 



7 + 7 ^ vr . 



can be produced in ttotto charge exchange scattering: 



7r° + 7r° ^ 7r+ + vr" 



as well as in two photon, and in electron-positron fusion processes 



:i.69) 



;i.70) 



7 -|- 7 — > vr 



e + e TT + vr" 



:i.7i) 
:i.72) 



We find that for vr^ production, the last two processes are much slower compared to the 
first, in case that ttq density is near chemical equilibrium. Similarly, the two photon fusion 
to two 7r°: 

7 + 7 7r° + 7r°, (1.73) 
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as expected, has rate much smaller than rate of one vr" production. It is a reaction of 
higher order in a and the energy is shared between two final particles. 

In the plasma under consideration, muons can be directly produced in the reac- 
tions: 



We will show in section 16.31 that already at temperatures > 5 MeV the large yield 
of pion and muon can be observed from e"'"e~7 plasma. 

These pion and muon production reactions take place in early universe and "pi^ 
production does not freeze-out with universe expansion, as we will show here in section [71 
This results that density of 7r° is relatively large, comparable to proton p and neutron n 
densities up to the temperatures of few MeV. vr^ and muon production freeze out about 
few MeV. However up to this temperatures their density is also comparable to n and p 
densities. 



7 + 7 /i+ + /i , 
e"^ e" ^ /i"^ + 



(1.74) 
(1.75) 
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CHAPTER 2 

STATISTICAL HADRONIZATION AND ESTIMATION OF PHASE 
SPACE OCCUPANCY FACTORS 



2.1 Introduction 

As we discussed in section 11.2.21 we assume fast hadronizaton. Then physical parameters 
of fireball as temperature, volume, entropy, strangeness and heavier flavor multiplicities do 
not change during hadronization. 

In this chapter we introduce the notion of conservation of entropy in section [2l2] and 
strangeness in section 12.31 expected to be valid in the fast hadronization process at LHC, 
and discuss how this impacts the SHM statistical parameters. We consider the entropy 
in a system with evolving strangeness in subsection 12.2.11 and show that the number of 
active degrees of freedom in a QGP is nearly constant. Another highlight is the discussion 
of sudden hadronization of strangeness and the associated values of hadron phase space 
parameters in subsection I2.3.4[ Throughout this paper we will use explicitly and implicitly 
the properties of QGP fireball and hadron phase space regarding entropy and strangeness 
content developed in these two sections 12.21 and 12.31 

2.2 Entropy conservation at hadronization 
2.2.1 Number of degrees of freedom in QGP 

The entropy density following from Eq. (11.551) is: 

Ivr^ 7 4 A 

^ = -^i9g + -^g,)T' + -,'^^/.gsT' + -. (2.76) 

For strange quarks in the second term in Eq. (I2.76P we set the entropy per strange quarks 
to 4 units. In choosing Ss/Ng = 4 irrespective of the effect of interaction and mass value 
rris/T we are minimizing the influence of unknown QCD interaction effect. 

The last term in Eq. (I2.76P comes from differentiation of the strong coupling con- 
stant as in the partition function with respect to T, see Eq. (ll.55p . Up to two loops in the 
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/9-function of the renormalization group the correction term is [30]: 



A = {boa^ + bia 



3 18 



(2.77) 



with Uf being the number of active fermions in the quark loop, n/ ~ 2.5, and 



For the strong couphng constant as we use 

^ YT^h^rTT^' ^ = 0-760 ±0.002, (2.79) 

where = 0.16 GeV. This expression arises from renormahzation group running of asif^), 
the energy scale at /i = 27rT, and the value as{Mz) = 0.118. A much more sophisticated 
study of the entropy in the QGP phase is possible [31], what we use here is an effective 
model which agrees with the lattice data [32]. 

Eq. fl2.76p suggests that we introduce an effective degeneracy of the QGP based on 
the expression we use for entropy: 

&T) = g,{T) + \g,{T) + 2^.^ + (2.80) 

Which allows us to write: 

^ = %9?,T^ (2.81) 

and 

dS Att^ o ^dV ^ 



We show the QGP degeneracy in figure[231 as a function of T G [140, 260] MeV, top 
frame for fixed s/S = 0, 0.03, 0.04 (from bottom to top), and in the bottom frame for the 
strangeness chemical equilibrium, 7^ = 1 (dashed) and approach to chemical equilibrium 
cases (solid). When we fix the specific strangeness content s/S in the plasma comparing 
different temperatures we find that in all cases g^g increases with T. For s/S = we have a 
2-fiavor system (dotted line, red) and the effective number of degrees of freedom g^^ varies 
between 22 and 26. The solid line with dots (green) is for s/S = 0.03, and the dot-dashed 
line (blue) gives the result for s/S = 0.04. 
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Figure 2.4: (color on-line) The Stefan-Boltzman degrees of freedom g^g based on entropy 
content of QGP, as function of temperature T. Upper frame: fixed s/S, the solid line with dots 
(green) is for a system with fixed strangeness per entropy s/S = 0.03, while dot-dashed (blue) 
line is for s/S = 0.04. The dotted (red) line is for 2-flavor QCD s/S = {u,d,G only); The 
bottom frame shows dashed (black) line 2-|-l-flavor QCD with nig = 125 it 35 MeV (chemically 
equilibrated u,d,s,G system). The (thick, thin) solid lines are for QGP in which strangeness 
contents is increasing as a function of temperature, see text. 
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In the bottom panel of figure [2^ we note that hke for 2 fiavors, case (s=0), for the 
2+1-fiavor system (7^ = 1) g^^ increases with T (dashed fine, black), g^^ varies between 
30 and 35.5. The thin dashed lines indicate the range of uncertainty due to mass of the 
strange quark, which in this calculation is fixed with upper curve corresponding to = 90 
MeV, and lower one rus = 160 MeV. The expected decrease in value of with T will thus 
have the effect to steepen the rise in the degrees of freedom with T. 

We now explore in a QGP phase the effect of an increasing strangeness fugacity 
with decreasing temperature. This study is a bit different from the rest of this paper, 
where we consider for comparison purposes hadronization for a range of temperatures but 
at a fixed value of s/S. A variable 7P(r) implies a more sophisticated, and thus more 
model dependent picture of plasma evolution. However, this offers us an important insight 
about g%. 

We consider the function: 

This is consistent with the kinetic computation of strangeness production [26]. At T = 140 
MeV we have chemical equilibrium in the QGP phase, while and at the temperature 
T = 260 MeV we have 7^ = 0.25. The result for is shown as a thick (black) line 
in figure 12. 4[ with the range showing strange quark mass range = 125 ± 35 MeV. We 
see that in a wide range of temperatures we have 29.5 < g^g < 30.5. 

The lesson is that with the growth of 7^ with decreasing T the entropy of the 
QGP is well described by a constant value g^^ = 30 ± 0.5. Since the entropy is (nearly) 
conserved and g^^ is (nearly) constant, Eq. (12.82^ implies that we can scale the system 
properties using the constraint T^dV/dy =Const. We stress again that these results arises 
in a realistic QGP with 2 + 7^-fiavors, but are model dependent and of course rely on the 
lattice motivated description of the behavior of QGP properties. On the other hand it is 
not surprising that the rise of strangeness chemical saturation with decreasing temperature 
compensates the 'freezing' of the g, G-degrees of freedom with decreasing temperature. 



50 



Table 2.1: Reference values of volume, temperature, entropy, particle multiplicity 



dV/dy [fm3] 


T[MeV] 


dS'^/dy 


dN'^/dy 


dN^/dy 


800 


200 


10,970 


2,700 


5,000 


2300 


140 


10,890 


2,700 


4,500 



2.2.2 Entropy content and chemical (non-) equilibrium 

We use as a reference a QGP state with dV/dy = 800 fm^ at T = 200 MeV, see table O 
We find from Eq. fll.5ip the Q and H phase particle multiplicity. The hadron multiplicity 
stated is what results after secondary resonance decays. The total hadron multiplicity 
after hadronization and resonance decays was calculated using on-line SHARE 2.1 [57]. If 
a greater (smaller) yield of final state hadrons is observed at LHC, the value of dV/dy need 
to be revised up (down). In general expansion before hadronization will not alter dS/dy. 
We can expect that as T decreases, V^^^ increases. Stretching the validity of Eq. (12.821) to 
low temperature T = 140 MeV, we see the result in the second line of table 16.61 

For QGP, in general the entropy content is higher than in a comparable volume 
of chemically equilibrated hadron matter, because of the liberation of color degrees of 
freedom in the color-deconfined phase. The total entropy has to be conserved during tran- 
sition between QGP and HG phases, and thus after hadronization, the excess of entropy is 
observed in excess particle multiplicity, which can be interpreted as a signature of decon- 
finement [59,60]. The dynamics of the transformation of QGP into HG determines how 
this additional entropy manifests itself. 

The comparison of entropy in both phases is temperature dependent but in the 
domain of interest i.e. 140 < T < 180 MeV the entropy density follows: 

> 3a". (2.84) 

Since the total entropy S is conserved or slightly increases, in the hadronization process 
some key parameter must grow in the hadronization process. There are two options: 
a) either the volume changes: 

3yH > -^Q. ^2.85) 

or 
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(2.86) 



In a slow, on hadronic time scale, transition, such as is the case in the early Universe, we 
can expect that case a) prevails. In high energy heavy ion collisions, there is no evidence 
in the experimental results for the long coexistence of hadron and quark phases which is 



subsection 12.3.41 

2.3 Strangeness in Hadronization 

2.3.1 Strangeness abundance in QGP and HG 

The efficiency of strangeness production depends on energy and collision centrality of heavy 
ions collisions. The increase, with value of centrality (participant number), of per-baryon 
specific strangeness yield indicates presents of strangeness production mechanism acting 
beyond the first collision dynamics. The thermal gluon fusion to strangeness can explain 
this behavior [16], and a model of the fiow dynamics at RHIC and LHC suggests that the 
QGP approaches chemical equilibrium but also can exceed it at time of hadronization [26]. 

The strangeness yield in chemically equilibrated QGP is usually described as an 
ideal Boltzman gas. However, a significant correction is expected due to perturbative QCD 
effects. We implement this correction based on comments below Eq. (11.631) . We use here 
the expression: 



with the Boltzman limit density Eq. (11.101) . and mass rus = 125 MeV, Qs = 6, = 1. The 
QCD correction corresponds to discussion of entropy in subsection 11.3.41 

We obtain strange quarks phase space occupancy 7^ as a function of temperature 
from condition of equality of the number of strange quark and antiquark pairs in QGP and 



required for volume growth. Consequently, we have ~ and a large value of 7^ is 
required to conserve entropy. The value of 7^ is restricted by 




(2.87) 
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HG. Specifically, in the sudden QGP hadronization, quarks recombine and we expect that 
the strangeness content does not significantly change. For heavier fiavors across the phase 
boundary this condition Eq. fll.lSp is very well satisfied, for strangeness the fragmentation 
effect adds somewhat to the yield, 

dN? ^ dN9 



> 



(2.89) 



dy ~ dy 

Using the equality of yields we underestimate slightly the value of strangeness occupancy 
that results. We recall that we also conserver entropy Eq. (I1.14p which like strangeness can 
in principle grow in hadronization. 



> 1 

CI 

H o 



Q 



(2.90) 



using Eq. fll.l4p we underestimate the value of 7^^. 

Counting all strange particles, the number of pairs is: 

dN'? dV 



H / H eq I H 2 eq\ 

dy - dy^ {%^K+% V) 

+ 7^(27?%^ + n; + P.< 
+3 if'rQ], (2.91) 



where rf^ are densities of strange hadrons (mesons and baryons) calculated using Eq. (11. 9p 
in chemical equilibrium. Pg is the strangeness content of the rj. The way we count hadrons 
is to follow strangeness content, for example = + n^o = + '^^o- We impose 
in our calculations s = s. The pattern of this calculation follows an established approach, 
SHARE 2.1 [57] was used in detailed evaluation. 

2.3.2 Strangeness per entropy s/S 

Considering that both strangeness, and entropy, are conserved in the hadronization process, 
a convenient variable to consider as fixed in the hadronization process, is the ratio of these 
conserved quantities s/S. In chemical equilibrium we expect that in general such a ratio 
must be different for different phases of matter from which particles are produced [26,61,62]. 

We compare QGP and HG specific per entropy strangeness content in figure 12.51 
We show as function of temperature T the s/S ratios for chemically equilibrated QGP and 
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Figure 2.5: (Color on line) Strangeness to entropy ratio s/S as function of temperature T, for 
tlie QGP (green, solid line for nis = 160 MeV, blue dash-dot line for mg = 90 MeV) with k = 1, 
see Eq. (j2.88p : and for HG (light blue, dashed line) phases for 7^ = 7^ = Ag = = 1 in both 
phases. 

HG phase. For the QGP the entropy S in QGP is calculated as described in section [2^2} and 
we use A; = 1 in Eq. (12.88^ . The shaded area shows the range of masses of strange quarks, 
considered, results for nis = 90 MeV (upper (blue) dash-dotted line) and rris = 160 MeV 
((green) solid line) form the boundaries. The central QGP value is at about s/S = 0.032. 

The short-dashed (light blue) line shows the hadron phase s/S value found using 
SHARE 2.1 . For HG near to usual range of hadronization temperature T ~ 160 MeV we 
find s/S* ~ 0.025. In general formation of QGP implies and increase by 30% in s/S. Both 
HG and QGP phases have a similar specific strangeness content at T = 240-260 MeV, 
however it is not believed that a HG at such high temperature would be a stable form of 
matter. This HG to QGP dissociation, or QGP hadronization depends on the degree of 
strangeness equilibration in plasma [66], and other dynamical factors. 

In the QGP the value of s/S for the range of realistic hadronization temperature 
140 < T < 180 MeV is in general larger than in HG. This implies that generally, the 
abundance of strange hadrons produced in hadronization over saturates the strange hadron 
phase space, if QGP state had reached (near) chemical equilibrium. Moreover, since we 
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are considering the ratio s/S and find in QGP a value greater than in HG, for chemical 
equilibrium in QGP the hadronization process will lead to 7^/7^ > 1. 

One can wonder if we have not overlooked some dynamical or microscopic effect 
which could adjust the value of s/S implied by QGP to the value expected in HG. First we 
note that the fast growth of the volume V cannot change s/S. Moreover, any additional 
strangeness production in hadronization would enhance the over-abundance recorded in 
the resulting HG. Only a highly significant entropy production at fixed strangeness yield 
in the hadronization process could bring the QGP s/S ratio down, masking strangeness 
over-saturation. A mechanism for such entropy production in hadronization is unknown, 
and moreover, this would further entail an unexpected and high hadron multiplicity excess. 

One could of course argue that the perturbative QCD properties in the QGP are 
meaningless, the entropy in QGP is much higher at given temperature. However, the 
properties of QGP have been checked against the lattice results, and the use of lowest 
order expressions is justified in these terms [32]. Moreover, the value of s/S" is established 
way before hadronization. 

2.3.3 Wroblewski ratio Ws 

At this point it is appropriate to look at another observable proposed to study strangeness 
yield, the Wroblewski ratio [63]: 

Ws = , (2.92) 

{uu) + {dd) ^ ' 

Ws compares the number of newly produced strange quarks to the produced number of light 
quarks. In an equilibrated deconfined phase Ws compares the number of active strange 
quark degrees of freedom to the number of light quark degrees of freedom. 

The ratio s/S compares the strange quark degrees of freedom to all degrees of 
freedom available in QGP. Therefore as function of T the ratios s/S and Ws can behave 
differently: Considering the limit T ^ a. constant s/S indicates that the reduction of 
s-degrees of freedom goes hand in hand with the 'freezing' of gluon degrees of freedom, 
which precedes the 'freezing' of hght quarks. This also imphes that for T — > Tc in general 
Ws diminishes. The magnitude of mg, the strange quark mass decisively enters the limit 
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For T » T(. the ratio Wg can be evaluated comparing the rates of production of 
hght and strange quarks, using the fluctuation-dissipation theorem [64] , which allows to 
relate rate of quark production to quark susceptibilities Xi (see Eqs. (11) and (12) in [64]): 

An evaluation of R^^ as function of temperature in lattice QCD has been achieved [65]. 
For T ^ 2.5Tc the result obtained, Wg —>■ R^ — 0.8, is in agreement with the expectation 
for equilibrium QGP with nearly free quarks, with mass of strangeness having a small but 
noticeable significance. With decreasing T, the ratio R^ 0.3 for = T^. However, this 
value of TUs is too large, the physical value should be nearly half as large, which would 
result in a greater R^. Moreover, for T ^ Tc the relationship of Ws to R^, Eq. fl2.93p 
is in question in that the greatly reduced rate of production of strangeness may not be 
satisfying the conditions required in Ref. [64]. 

Comparing the observables s/S and Ws we note that the experimental measure- 
ment requires in both detailed analysis of all particles produced. At lower reaction 
energies there is additional complication in evaluation of Wg due to the need to subtract the 
effect of quarks brought into the reaction region. Turning to the theoretical computation 
of s/S and Wg we note that the thermal lattice QCD evaluation of s/S is possible without 
any approximation, even if the actual computation of entropy near the phase boundary 
is a challenging task. On the other hand, the lattice computation of Wg relies on pro- 
duction rate of strangeness being sufficiently fast, which cannot be expected near to the 
phase boundary. Moreover, the variable s/S probes all QGP degrees of freedom, while Wg 
probes only quark degrees of freedom. We thus conclude that s/S is both more accessible 
theoretically and experimentally, and perhaps more QGP related observable, as compared 
to Wg, since it comprises the gluon degrees of freedom. 

2.3.4 Strangeness chemical non-equilibrium 

In order that in fast hadronization there is continuity of strangeness Eq. (12. 89 p . and entropy, 
Eq. f l4.130p the hadron phase 7^ 7^ 1 and 7^ 7^ 1. We have to solve for 7^ and 7^ 
simultaneously Eqs. ^MSdM). 

In figure 12.61 we show as a function of T the strange phase space occupancy 7^, 
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- - 1^, s/S=0.045 

- ■ ^1^, s/S=0.04 

- - i^,\/S=0.035 
— . ^7^, s/S=0.03 

y", s/S=0.025 



. . . Y"=exp(m /2T) 

H q " 




140 160 180 200 220 240 260 
T [MeV] 



Figure 2.6: (color on-line) Phase space occupancy as a function of T: 7^ (blue, solid line), 7^ 
(dash-dotted lines, from top to bottom) for s/S = 0.045, for s/5 = 0.04 (thick line), s/S = 0.035, 
s/S = 0.03 (thick line), s/S = 0.025; 7^'' (red, dashed line). 

obtained for several values of s/S ratio (from top to bottom 0.045, 0.04, 0.035, 0.03, 0.025) 
evaluated for S'^ = S^. The solid line shows 7^ for = and S'^ = S^. The maximum 
allowed value Eq. fl2.87p is shown dashed (red). 

In figure [2771 we show results for 7^/7^ (where 7^ = 1 we show 7^ ) . We consider 
the three cases: 7^ = 1 , 7^ = ■y^^, and entropy conservation = for s/S = 0.045, 
s/S = 0.04, s/S = 0.035, s/S = 0.03, s/S = 0.025 (dash-dot lines) (lines from top to 
bottom). We see that except in case that strangeness were to remain well below chemical 
equilibrium in QGP {s/S ~ 0.03), the abundance of heavy flavor hadrons we turn to 
momentarily will be marked by an overabundance of strangeness, since practically in all 
realistic conditions we find 7^ > 7^- 

In figure [23] we show s/S ratio as function of 7^/7^. The solid line is for T = 200 
MeV, 7j^ = 0.83, S^ = S^ , dashed line for T = 170 MeV, -ff = 1.15 S^ = S^ and 
dash-dot line for T = 140 MeV, 7^ = 1.6 MeV, S^ = S^ . We also consider 7^ = 1 case 
for T = 170 MeV (dot marked (purple) solid line). In this case strangeness content 7s/7g 
is higher than for S^ = S^ with the same T and s/S. In the limit 7^ = 7^*^, Eq. 12.871 
(T = 200 MeV, solid, thin line; T = 170 MeV, dashed line) the strangeness content 7^/7^ 
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Figure 2.7: (color on line) 7^/7^ (=7? at 7^ = 1) as a function of hadronization temperature 
T. Top frame: 7^ = 1, middle frame: 7^ = 7^*^, and bottom frame: = S^. Lines, from 
top to bottom: s/S = 0.04 (blue, solid line), s/S = 0.03 (green, dashed line), s/S = 0.025 (red, 
dash-dot line) 
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Figure 2.8: (color on line) Strangeness to entropy ratio, s/S, as a function of ^s/jq for T = 200 
MeV, 7g = 0.083, = 5«(solid line), T = 170 MeV, 7g = 1.15, = (dashed line), T = 140 
MeV, 7g = 1.6, = S"3 (dash-dotted line); 7^ = 1 (dot marked sohd); 7, = j^'': T = 200 MeV 
(thin solid line), T = 170 MeV (thin dashed line). 



is minimal for given T and s/S. 

These results suggest that it is possible to measure the value of s/S irrespective 
of what the hadronization temperature may be, as long as the main yield dependence is 
on the ratio 7l^/7?- Indeed, we find that the ratio 4>/K~^; 



K+ 



(2.94) 



is less sensitive to hadronization temperature compared to its strong dependence on the 
value of s/S. In figure [2^9] we show the total hadron phase space ratio ip/K^ as function 
of T for several s/S ratios, and for 7^^ = 1 (chemical equilibrium, dashed (red) line). The 
yield contains the contribution from the decay of (f) into kaons which is a noticeable 
correction. 

We record in table [2^ for given s/S and volume dV/dy the corresponding total 
yields of strangeness, which may be a useful guide in consideration of the consistency of 
experimental results with what we find exploring heavy flavor hadron abundance. 
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Figure 2.9: (Color on line) The ratio as a function of T. Dashed line (red) is for 

chemical equilibrium. Solid line with dots (green) s/S = 0.03, solid line (blue) s/S = 0.04, 
dash-dot line (per) is for s/S = 0.022. 



Table 2.2: Specific and absolute strangeness yield for different reaction volumes at T = 200 
MeV. 



s/S 


ds/dy 


dV/dy [fm-3] 


T [MeV] 


0.045 


550 


1000 


200 


0.04 


360 


800 


200 


0.035 


250 


700 


200 


0.03 


165 


600 


200 


0.025 


106 


500 


200 


0.022 


83 


500 


200 
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2.3.5 Phase space occupancy 7^ and 7^ 

The first step in order to determine the yields of heavy flavor hadronic particles is the 
determination of the phase space occupancy 7^ and 7^. 7^ is obtained from equality of 
number of these quarks (i.e. of quark and anti quark pairs) in QGP and HG. The yield 
constraint is: 

^ = ^ K"<P + ifinlt^ + 27?nS, + ^ifnlD] ; (2.95) 
where open 'op' charm yield is: 

= ifnS + ifnL + l^n^^ + ^f^^nl^ + 7?^ n^^,- (2.96) 

Here and n^^ are densities of D and Dg mesons, respectively, in chemical equilibrium, 
n^^^ is equilibrium density of baryons with one charm and two light quarks, n^^^ is density 
of baryons with one charm (or later on one bottom quark) and two strange quarks (f2°, 
f2°) in chemical equilibrium and ?t,^?j is equilibrium particle density with both, a charm 
(or bottom) and an anticharm (or antibottom) quark (C=0, B=0, S=0). The equilibrium 
densities can be calculated using Eq. fll.9p . 7^ can now be obtained from Eq. (12.951) . 

Similar calculations can be done for 7^. The only difference is that we need to 
add number of Be mesons to the right hand side of Eq.f l2.95i) . 

^ = ^H^Hn-L^ — (2 97) 

dy ^^^^ ^^rfy- ^ ^ 

is density in chemical equilibrium of Be- In the calculation of 7^ the contribution of 

term with ubc is very small and we did not consider it above. 

The value of 7c is in essence controlled by the open single charm mesons and 

baryons. For this reason we do not consider the effect of exact charm conservation. The 

relatively small effects due to canonical phase space of charm are leading to a slight up- 

renormalization of the value of 7c so that the primary dNd dy yield is preserved. This effect 

enters into the yields of multi-charmed and hidden charm hadrons, where the compensation 

is not exact and there remains slight change in these yields. However, the error made 

considering the high yield of charm is not important. On the other hand for multi-bottom 

and hidden bottom hadrons the canonic effect can be large, depending on actual bottom 

yield, and thus we will not discuss in this paper yields of these hadrons, pending extension 

of the methods here developed to include canonical phase space effect. 
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We consider in figure 12.101 the temperature dependence of both 7^ (top) and 7^ 
(bottom) for the heavy flavor yield given in Eqs. ( ]3.99|3.100p . In the non-equilibrium 
case (solid lines) the space occupancy 7^ is obtained from Eq. (12.911) and 7^ is chosen to 
keep = S^. 7^^^ depend on 7^ and 7^: the value of A''^ in Eq. (l2.9ip is chosen to have 
s/S = 0.04 after hadronization, the corresponding 7^ and 7^ are shown in figures [2l6l and 
12.71 . Since applicable 7^ may depend on hadronization dynamics and/or details of equation 
of state of QGP, we show charm quark phase space occupancies also for maximum possible 
value of 7^ — > 7g^), also considered at s/S = 0.04 for all hadronization temperatures. We 
can compare our results with the chemical equilibrium (dashed lines) setting 7^ = 7^ = 1 
in Eq. (12.951) . At hadronization condition T = 160 ± 20 MeV temperatures we see in 
figure [2?T0] a significant (considering the fast changing logarithmic scale) difference between 
the chemical equilibrium, and non-equilibrium (s/S=0.04) results. 

In figure [2TTT] we show the ratio 7^{,)/7^fe)eg ^ function of hadronization temper- 
ature T. This helps us understand when the presence of chemical nonequilibrium is most 
noticeable. This is especially the case should heavy flavor hadronization occur at the same 
temperature T = 140-170 MeV as is obtained for non-heavy hadrons, and/or when the 
entropy content of light hadrons is maximized with 7^ 7^"^. When no additional entropy 
is formed in hadronization, that is = S'^, 'lf(^h)/l4b)eq exceeds unity for T > 200 MeV, 
at which point the heavy flavor hadron yields exceed the chemical equilibrium expecta- 
tions. In general we find that heavy hadron yields if produced at normal hadronization 
temperature would be effectively suppressed, compared to statistical equilibrium results, 
by the high strangeness yield. This happens since the phase space is bigger at 7^^ > 1, 
and thus a smaller 7^^ is required to reach a given heavy flavor yield. 

7^ and 7^ are nearly proportional to dNi,^c/dy, respectively. The deviation from 
the proportionality is due to the abundance of multi-heavy hadrons and it is small. To 
estimate this effect more quantitatively we first evaluate: 




(2.98) 



i. e. the value expected in absence of multi heavy hadrons. Next we compare with the result 
when we take into account the last three terms in Eq. (I2.95p . The influence of these therms 
depend not only on dNc/dy but also on dNc/dy /dV/dy. For fixed dV/dy = 800 fm~^ in the 
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Figure 2.10: (Color on line) 7^(6 = 1) (upper panel), and 7^ (c=10) (lower panel), as functions 
of temperature of hadronization T. The solid lines are non-equilibrium for s/S = 0.04 with 
SQ = S^, dashed hues are equilibrium case 7s=7g=l and dot-dash lines are ioi s/S = 0.04 with 
maximal value of 7g(7g = 7^*") {dV/dy = 800 fm^). 
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Figure 2.11: (Color on line) J^/"f]^eq 7?/7ceg' ^ functions of temperature of heavy flavor 
hadronization T. The solid with dot marks line is for "y^/j^gg with s/S = 0.04, solid line is for 7^ 
with s/S = 0.04, dot-dash and dashed lines are for s/S = 0.04 with maximal value of 7^ 7^'' 
for 7^/7^eq and for 7^/7^eg) respectively. 

range of dN^/ dy = (5, 30), we find that 'j^/Nc (and therefore yields of open charm hadrons) 
changes at temperature T = 140 MeV by ~ 6% for the s/S = 0.04. For the chemical 
equilibrium case 7^ = 7^ = 1, 7^/Ai'c changes up to 15% at the same conditions. For the 
particles with hidden charm or 2 charm quarks the yields are proportional 7?, therefore 
changes in their yields will be about twice larger. For RHIC iVc < 3 and dV/dy = 600 
fm~'^ the dependence of yields on A^^^ is much smaller. 

The multiplicity dN^/dy can also influence 7^, since as we noted it also includes 
a term proportional to 7^^^^- range of N^. = (5, 30), 'jb/Nt changes at temperature 

T = 0.14 MeV by ~ 0.5% for s/S = 0.04. Since the mass of 6-quark is much larger than 
that of c-quark, the effect due to multi-bottom states is negligible. 
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CHAPTER 3 

HEAVY FLAVOR HADRONS IN STATISTICAL HADRONIZATION OF 
STRANGENESS AND ENTROPY RICH QGP 



3.1 Introduction 

A relatively large number of hadrons containing charmed and bottom quarks are expected 
to be produced in heavy ion ( AA) coUisions at the Large Hadrons CoUider (LHC) . Because 
of their large mass c, c, b, b quarks are produced predominantly in primary parton-parton 
collisions [52], at RHIC [53], and thus even more so at LHC. These heavy flavor quarks 
participate in the evolution of the dense QCD matter from the beginning. In view of 
the recent RHIC results it can be hoped that their momentum distribution could reach 
approximate thermalization within the dense QGP phase [54]. 

In the calculations in this chapter we assume the same evolution stages as present 
in the beginning of dissertation introduction, except the kinetic. Here we assume that this 
phase does not have influence on charm and bottom hadrons. 

It is important to observe that in the presence of deconflned QGP phase heavy 
hadrons containing more than one heavy quark are made from heavy quarks created in dif- 
ferent initial NN collisions. Therefore yields of these hadrons are expected to be enhanced 
as compared to yields seen in single NN collisions [25,27]. We note that the Bc(bc,bc) and 
J/'^{cc) and more generally all bound cc states yields were calculated before in the kinetic 
formation and dissociation models [24,25]. Our present work suggests that it is important 
to account for the binding of heavy flavor with strangeness, an effect which depletes the 
ehgible supply of heavy flavor quarks which could form Bc(6c, fee) and J/'^{cc) [2]. 

Enhanced production yield of multi-heavy hadrons can be considered to be an 
indicator of the presence of deconfined QGP phase for reasons which are analogue to those 
of multi-strange (anti) baryons [16]. Considering that we have little doubt that QGP is 
the state of matter formed in the very high energy AA interactions, the study of yields 
of multi-heavy hadrons is primarily explored in this work in order to falsify, or justify. 
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features of the statistical hadronization model (SHM) employed or the model itself in the 
context of formation of the heavy flavor hadrons. 

For example, differing from others recent studies which assume that the hadron 
yields after hadronization are in chemical equilibrium [27,28], we form the yields based on 
abundance of u. d, s quark pairs as these are available at the chemical freeze-out (particle 
formation) conditions in the quark-gluon phase. This approach is justified by the expecta- 
tion that in a fast break-up of the QGP formed at RHIC and LHC the phase entropy and 
strangeness will be nearly conserved during the process of hadronization. We will investi- 
gate in quantitative terms how such chemical non-equilibrium yields, in the conditions we 
explore well above the chemical equilibrium abundance, influence the expected yields of 
single, and multi-heavy flavor hadrons. 

In the order to evaluate the yields of final state hadrons we enforce conservation of 
entropy, and the flavor s, c, b quark pair number during phase transition or transformation. 
The faster the transition, the less likely it is that there is signiflcant change in strange quark 
pair yield. Similarly, any entropy production is minimized when the entropy rich QGP 
breakup into the entropy poor HG occurs rapidly. The entropy conservation constraint 
flxes the flnal hght quark yield. We assume a fast transition between QGP and HG phases, 
such that all hadron yields are at the same physical conditions as in QGP breakup. 

In the evaluation of heavy particle yields we form ratios involving as normalizer 
the total heavy flavor yield, and for yields of particles with two heavy quarks we use as 
normalizer the product of total yields of corresponding heavy flavors such that the results 
we consider is as little as possible dependent on the unknown total yield of charm and 
bottom at RHIC and LHC. The order of magnitude of the remaining dependence on heavy 
flavor yield is set by the ratio of yield of all particles with two heavy quarks to yield 
of particles with one heavy quark. This ratio depends on the density of heavy flavor at 
hadronization, [dNc/dy) /{dV /dy). The results we present for LHC are obtained for an 
assumed charm and bottom quark multiplicity: 

dN 

^ = c = 10, (3.99) 
dy 

^ = 6 = 1. (3.100) 
dy 

and dV/dy — 800 fm^ at T = 200 MeV. Theoretical cross sections of c and b quarks 
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production for RHIC and LHC can be found in [55,56]. In certain situations we will 
explore how variation of the baseline yields Eq. ( ]3.99p and Eq. ( ]3.100p impact the results. 
In particular among the yields of multi-heavy hadrons, this influence can be noticeable, 
see discussion in the end of section I2.3.5[ We note that the number of b quarks can not 
change during expansion, because of large mass rrih » T. It is nearly certain that all 
charm in QGP at RHIC is produced in the first parton collisions, for further discussion 
of LHC see Ref. [26] - it appears that for all practical purposes also in the more extreme 
thermal conditions at LHC charm is produced in the initial parton interactions. 

In order to form physical intuition about the prevailing conditions in the QGP 
phase at time of hadronization, we also evaluate the heavy quark chemical reference density, 
that is the magnitude of the chemical occupancy factor in QGP, considering the pre- 
established initial yields of c and b from parton collision. For this purpose we use in the 
deconfined QGP phase: 

TUc = 1.2 GeV, 
TJib = 4.2. GeV 

We also take Aj = l,i = u,d,s for all light flavors, since the deviation from particle- 
antiparticle yield symmetry is rather small and immaterial in the present discussion. 

When computing the yields of charmed (and bottom) mesons we will distinguish 
only strange and non-strange abundances, but not charged with non-charged (e.g. D~{cd) 
with D^{cu)). We assume that the experimental groups reporting results, depending on 
which types of D-meson were observed, can infer the total yield (charged+non-charged) 
which we present. We treat in similar way other heavy hadrons, always focusing on the 
heavy and the strange flavor content and not distinguishing the light flavor content. 

This chapter is organized as follows: we use the elements of the SHM model 
introduced in section [l.2.2l to evaluate heavy flavor hadron yields. This allows us to discuss 
the relative yields of strange and non-strange heavy mesons in section 13.21 and we show 
how this result relates the value of the strangeness chemical (non-) equilibrium parameters. 
In this context, we also propose a multi-particle ratio as a measure of the hadronization 
temperature, and explore how a multi-temperature, staged, freeze-out would impact the 
relevant results. 
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We turn to discuss the heavy flavor hadron yields for given bulk QGP constraints 
in section [373I where we also compare when appropriate to the strangeness and light quarks 
chemical equilibrium results. We use the charm and bottom quark phase space occupancy 
parameters (subsection 12.3. 5p and turn in subsection 13.3.11 to discussion of the yields of 
single heavy mesons, which we follow with discussion of yields of single heavy baryons in 
subsection 13.3.21 In last subsection 13.3.31 we present the expected yields of the multi-heavy 
hadrons, in so far these can be considered in the grand canonical approach. We conclude 
our work with a brief summary in section 13. 4[ 

3.2 Relative Charmed Hadron Yeilds 
3.2.1 Determination of 7s/7g 

We have seen considering s/S and also s and S individually across the phase limit that 
in general one would expect chemical non-equilibrium in hadronization of chemically equi- 
librated QGP. We first show that this result matters for the relative charm meson yield 
ratio D/Ds, where Ds{cs) comprises all mesons of type (cs) which are listed in the bottom 
section of table 13.31 and D{cq) comprise yields of all (eg) states listed in the top section of 
table 13. 3[ This ratio is formed based on the assumption that on the time scale of strong 
interactions the family of strange-charmed mesons can be distinguished from the family 
non-strange charmed mesons. 

The yield ratio D/Dg calculated using Eq. (11.91) and Eq. (ll.lOl) is shown in fig- 
ure 13.121 Using Eq. (I1.12p we see that this ratio is inverse proportional to 7^/7^ and 
weakly dependent on T: 



A deviation of 'js/lg from unity in the range we will see in section [2.3.41 leads to a noticeable 
difference in the ratio D/Ds- We show in figure [3?T2] results for T = 140, 160, 180 MeV. In 



D_ _ 7g T.igDsi'mj^^iexp{-mDsi/T) 
Ds 7? T,igDim^Df exp{-mDi/T) 



fiT) 




(3.101) 



this temperature range the effect due to 7^ / Yq^^ 7^ 1 is the dominant contribution to the 
variation of this relative yield. 
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Table 3.3: Open charm, and bottom, hadron states we considered. States in parenthe- 
sis either need confirmation or have not been observed experimentally, in which case we 
follow the values of Refs. [67, 68] . We have charm-bottom symmetry required for certain 
observables. 
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Figure 3.12: (Color on line) D/Dg ratio as a function oi'jf/j^ for T = 140 MeV (blue, dashed 
line), T = 160 MeV (green, solid line) and T = 180 MeV (red, dash-dot line). 
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Figure 3.13: (Color on line) J/^ cfi/DsDs ratio as a function of hadronization temperature T. 
3.2.2 Check of statistical hadronization model 

We next construct a heavy flavor particle ratio that depends on hadronization tempera- 
ture only. To cancel the fugacities and the volume we consider the ratio J/^cj)/ DgDs in 
figure 13.131 Here J/^ yield includes the yield of (cc) mesons decaying into the J/^ . All 
phase space occupancies cancel since J/^ oc 7^^, (p oc 7^^, Dg oc 7^7^ and similarly 
Ds oc 7^7^. When using here the particle Ds{cs) and antiparticle Ds{cs) any chemical po- 
tentials present are canceled as well. However, for the LHC and even RHIC environments 
this refinement is immaterial. 

This ratio J/^(f)/ DgDs, turns out to be practically constant, within a rather wide 
range of hadronization temperature T, see figure 13.131 The temperature range we study 
140 < T < 280 MeV allows us to consider an early freeze-out of different hadrons. To be 
sure of the temperature independence of J/^(f)/DsDs we next consider the possibility that 
hadronization temperature T of charmed hadrons is higher than hadronization temperature 
To of (p. We study this question by exploring the sensitivity of the ratio J/^cj)/ DsDg to 
the two temperature freeze-out in figure [3TT^ see bottom three lines for Tq = 180, 160, 140 
MeV with 7^ from condition S**^ = 5*^, see figure [231 If charmed hadrons hadronize later, 
T > To, and T - Tq < 60 the change in J (p / D ratio is small (about 20%). If this 
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Figure 3.14: (Color on line) J/^ 4>{Tq)/DsDs ratio is evaluated at two temperatures, T for 
heavy flavor hadrons, and Tq for (p as a, function of T — Tq, with three values of Tq = 140, 160, 180 
MeV is considered with = . 

were to be measured as experimental result, 

^^ 0.35, (3.102) 

DsD, 

one could not but conclude that all particles involved are formed by mechanism of statistical 
hadronization. 

This interesting result can be understood, considering the behavior of the 
7^(To)/7^(T) ratio, which increases rapidly with increasing T — Tq (see the top three 
lines in figure [XH]). This ratio almost compensates the change in 0-yield, an effect we al- 
ready encountered in the context of the results we show below in figure [231 For large T — Tq 
the ratio J/^(f)/ DgDg begins to decrease more rapidly because 7^ increases for 5*^ = S'^, 
see figure (12. 6p . 
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3.3 Yields of heavy flavored hadrons 
3.3.1 D, Ds, B, Bs meson yields 

In next sections we will mostly consider particles yields after hadronization and we will omit 
superscript H in 7s. Considering Eq. ( 11. 9p . we first obtain 7c as a function of "is/lq ratio 
and T. Substituting this 7c and appropriate equilibrium hadron densities into Eq. (11 .Qp we 
obtain yields of D{B) and Ds{Bs), as functions of the 7s/7g ratio, at fixed temperature, 
which are shown on figure I3.15[ In the upper panel we show the fractional yields of charmed 
D/Nc and Dg/Nc mesons, and in the lower panel B/Nb and Bg/Nh for T = 200 MeV (solid 
line), T = 170 MeV (dashed line), T = 140 MeV (dash-dot line). Fractional yield means 
that these yields are normalized by the total number of charm quarks and, respectively 
bottom quarks A*";,, and thus tell us how big a fraction of available heavy fiavor quarks binds 
to non-strange and strange heavy mesons, respectively. Using figure 12^81 the ratio 7^/75 
can be related to the s/S ratio. 7^ was chosen to conserve entropy during hadronization 
process, see figure [231 In general the heavy non-strange mesons yield decreases and strange 
heavy meson yield increases with js/lq- The yields D,B and Ds,Bs are sum over exited 
states of D,B and Ds,Bs respectively, see table 13^ for the 'vertical tower' of resonances 
we have included. 

Using 7c, 7s, 7q at a given T (see figures [2^ 12.71 12.101) we have now all the inputs 
required to compute absolute and relative particle yields of all heavy hadrons which we can 
consider within the grand canonical phase space. When we consider chemical equilibrium 
case, we use naturally 7^ = 7^ = 1. 

In figure [3A6] we consider the yields shown in figure [3TT5] as a functions of hadroniza- 
tion temperature. The dashed blue and green lines were obtained for chemical equilibrium 
yields of D and Dg respectively. The extreme upper and lower lines are for fractional D 
and Dg yields with s/S = 0.03 (dot marked, blue and green lines, respectively), while 
the central lines are for s/S = 0.04 (solid, blue and green lines). Also we show fractional 
yields for maximal possible value 7^ — > 7^^, see figure [2T6] for 'j^'^(T)) (dash-dot lines) and 

We note that there is considerable symmetry at fixed T between the fractional 
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Figure 3.15: (Color on line) Upper panel, fractional charm meson yield, and lower panel, 
fractional bottom meson yields as a function of "fsMq ratio for fixed hadronization temperature 
T. Upper lines in each panel are for D(B) mesons, solid line is for T = 200 MeU, 7g = 1.1, dashed 
line is for T = 170 MeV, 7^ = 1.15 and dash-dot line is for T = 140 MeV, 7^ = 0.83 {S^ = S^). 
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Figure 3.16: (Color on line) Upper panel, fractional charm meson yield, and lower panel, 
fractional bottom meson yields. Equilibrium (dashed lines) and non-equilibrium for s/S = 0.03 
(point marked solid line) and s/S = 0.04 (solid line) for D/Nc (blue lines, upper panel); Dg/Nc 
(green lines, upper panel); for D/Nc and Dg/Nc with s/S = 0.03 and s/S = 0.04 for 7g = 
(dash-dotted lines); B/N^ (solid line, lower panel); and Bg/Ni, (point marked solid line, lower 
panel), for B/N^ and Bg/Nh with s/S = 0.03 and s/S = 0.04 for 7^ = 7^'' (dash-dot lines); as a 
function of T. 
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Figure 3.17: (Color on line) Ratios D /Ds (upper panel) and B/Bs (lower panel) are shown as 
a function of T for different s/S ratios and in chemical equilibrium. Solid line is for s/S = 0.04, 
dash-dot line is for s/S = 0.03, dashed line is for 7^ = 7^ = 1. 
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yields of charmed, and bottom mesons, for the same condition of s/S. The chemical equi- 
librium results show significant difference between strange and non-strange heavy mesons. 
In the case of chemical equilibrium, for the considered very wide range of hadronization 
temperatures Dg/Nc — Bg/Nf, ~ 0.2 are nearly constant. A significant deviation from this 
result would suggest the presence of chemical non-equilibrium mechanisms of heavy flavor 
meson production. 

The yields of Dg/Nc and Bg/Nh are very similar, and similarly so for D/Nc and 
B /Niy. Thus the relative yield of either of these mesons measures the relative yield of charm 
to bottom participating in the statistical hadronization process: 

This is a very precise result, which somewhat depends on the tower of resonances included, 
and thus in particular on the symmetry in the heavy quark spectra between charmed and 
bottom states which we imposed. 

It is useful to reconsider here the ratio D/Dg (B/Bg) which is proportional to 
'jq/'js, see figure [3TT2] for Dg/D presented as a function of 7s/7g. We consider this ratio 
now as a function of T, the upper panel in figure 13.171 is for charm, the lower for bottom. 
We see that there is considerable symmetry in the relative yields between charmed and 
bottom mesons with upper and lower panels looking quasi-identical. Except for accidental 
values of T where the equilibrium results (blue, dashed lines) cross the fixed s/S results, 
there is considerable deviation in these ratios expected from chemical equilibrium. For 
LHC with s/ S = 0.04 this ratio is always noticeably smaller than in chemical equilibrium 
(solid purple line is for = and purple, dash-dot line is for 7^ = 7^''). Even for 
RHIC-like conditions with s/S = 0.03 this ratio is smaller than in chemical equilibrium for 
all temperatures when entropy conservation in hadronization is assumed, S^ = S^ (dot 
marked solid, green line). 

3.3.2 Heavy baryon yields 

As was the case comparing charm to bottom mesons we also establish a symmetric set of 
charmed and bottom baryons, shown in the table 13. 4[ Many of the bottom baryons are 
result of theoretical studies and we include that many states to be sure that both charm 
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Table 3.4: Charm and bottom baryon states considered. States in parenthesis are not 
known experimentally and have been adopted from references [69]. 



hadron 




M[GeV] 


hadron 




M[GeV] 


g 


A+(l/2+) 


udc 


2.285 


Ab0(l/2+) 


udb 


5.624 


2 


A+(l/2-) 


udc 


2.593 


A60(l/2-) 


udb 


(6.000) 


2 


A+(3/2-) 


udc 


2.6266 


A60(l/2-) 


udb 


(6.000) 


2 


E+(l/2+) 


qqc 


2.452 


S0(l/2+) 


qqb 


(5.770) 


6 


E*(3/2+) 


qqc 


2.519 


Eg*(3/2+) 


qqb 


(5.780) 


12 


Sc(l/2+) 


qsc 


2.470 


S.(l/2+) 


qsb 


(5.760) 


4 


Sl(l/2+) 


qsc 


2.5741 


s;(i/2+) 


qsb 


(5.900) 


4 


S.(3/2+) 


qsc 


2.645 


S;(3/2+) 


qsb 


(5.900) 


8 


0^(1/2+) 


ssc 


2.700 


f^b(l/2+) 


ssb 


(6.000) 


2 


^^c(3/2+) 


ssc 


(2.700) 


^^6(3/2+) 


ssb 


(6.000) 


4 



and bottom are consider in perfect symmetry to each other. In figure 13.181 (upper panel) 
we show hadronization temperature dependencies of yields of baryons with one charm 
quark normalized to charm multiplicity Nc- We show separately yields of baryons without 
strange quark (Ac + T,c)/Nc, and with one strange quark S=l (Ec/Nc). We show two cases 
for s/S = 0.04 with conserved entropy at hadronization = (solid lines) and with 
maximum possible entropy value 7^ = 7^'' (dash-dot lines). The chemical equilibrium case 
7g = 7s = 1 is also shown (dashed lines). The upper lines of each type are for {Ac + T,c)/Nc, 
the lower lines are for Sc/Nc- A similar result is presented for bottom baryons in the lower 
panel of figure 13.181 We note that the result for bottom baryons is more uncertain since 
most baryon masses entering are not experimentally verified. 

We note that the results shown figure 13.181 imply that under LHC conditions at 
least 15% of heavy flavor can be bound in heavy baryons, but possibly 30%. For large 
7g = 7^^ > 1 we see increase in (Ac + T,c)/Nc yields compared to chemical equilibrium 
and especially compared to entropy conserved hadronization = . This is so since 
yields are proportional to 7g, 7s7g. This results to relative suppression the DjNc (see 
figure 13.161) . 

In figure 13.191 we show ratio cqq/cqs = (Ac + Ec)/Hc as a function of 'Js/lq for 
T = 200 MeV (dash-dot line), T = 170 MeV (solid line) and T = 140 MeV (dashed line). 
This dependence is linear, the slope depends only on hadronization temperature T. The 
■js/lq ratio can be converted to s/S ratio using figure [2^ 
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Figure 3.18: (Color on line) Equilibrium (dashed lines), sjS = 0.04, S"^ = (solid lines), 
s/S = 0.04, 7g = 7^*^, (the upper panel) upper lines for each type are for ratio (Ac + T,c)/Nc 
and lower lines are for Sc/Nc (upper panel) and (lower panel) upper lines of each type are for 
(Afe + T,b)/Nc and lower lines are for 'Eh/Nf, as functions of T. 
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Figure 3.19: (Color on line) The ratios cqs/cqq = Hc/(Ac + Sc) (upper lines) and 
css/cqq = l^c/(Ac + S^) (lower lines) for T = 200 MeV (dash-dot hne), T = 170 MeV (sohd 
line) and T = 140 MeV (dashed line) as functions of 'Js/jq ■ 

The yield of multi-strange charmed baryon, Qc{css) is, similar to the light multi- 
strange hadrons, much more sensitive to chemical non-equilibrium. In figure 13.201 we see 
a large increase in fractional yield of Qc{css)/Nc for s/S = 0.04 and = (solid line) 
compared to the chemical equilibrium (dashed line) expectation for the entire considered 
range of hadronization temperature. As expected, this yields increase with T. This also 
means that higher formation temperature can be invoked to explain an unusually high 
yield. We expect that at LHC more than one percent of total charm yield will be found in 
the flc{css) state. 

3.3.3 Yields of hadrons with two heavy quarks 

We consider multi-heavy hadrons listed in the table 13.51 The yields we will compute are 
now more model dependent since we cannot completely reduce the result, it either remains 
dependent on the reaction volume dV/dy, or on the total charm(bottom) yields dN/dy. 
For example the yields of hadrons with two heavy quarks are approximately proportional 
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Figure 3.20: (Color on line) il,c{css)/Nc as function of T: dashed line for chemical equihbrium; 
solid lines are for = S^, dashed dotted lines are for 7^ = 7^*": both for s/S = 0.03 and 
s/S = 0.04 (upper lines). 

to l/{dV/dy) because 7^^ for heavy quarks is proportional to 1/dV/dy, see Eq. (12.951) : 



dN,i 



dy 
dNsc 



hid H2 



dV 



1 



oc 



dy dV I dy ' 

H udV 1 

« 7c 76 — oc 



(3.104) 
(3.105) 



dy dy dV/dy 

Moreover, unlike it is the case for single heavy hadrons, the canonical correction to grand- 
canonical phase space does not cancel out in these states, adding to the uncertainty. 

Thus the result we present must seen as a guiding the eye and demonstrating a 
principle. In figure 13.211 we show the yield of hidden charm cc mesons (see table 13.51) nor- 
malized by the square of charm multiplicity A*"^ as a function of hadronization temperature 
T. We consider again cases with s/S = 0.03 (upper panel) and s/S = 0.04 (lower panel), 
solid line is for S^ = S'^, dot-dash line is for 7^ = 7^^, and dot-dash line is for 7^ = 7^*^. 
The chemical equilibrium cc mesons yields are shown (dashed lines on both panels) for two 
different values of dV/ dy = 600 fm^ for T = 200 MeV (upper panel) and dV/ dy = 800 fm'^ 
for T = 200 MeV (lower panel). 

The yield of cc mesons is much smaller for s/S = 0.04 than in equilibrium for the 
same dV/dy for large range of hadronization temperatures. For s/S = 0.03 the effect is 
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Figure 3.21: (Color on line) ccjNl yields as a function of hadronization temperature T, at 
dVjdy = 600 fm-3 for T = 200 MeV, s/S = 0.03 (upper panel), dV/dy = 800 frn'^ for T = 200 
MeV, s/S = 0.04 (lower panel). Results shown are for S'^ = (solid lines), for 7^ = 7^'' 
(dash-dot lines), and for chemical equilibrium case (dashed lines, s/S is not fixed). 
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Table 3.5: Hidden charm and multi heavy hadron states considered. States in parenthesis 
are not known experimentally 



hadron 




mass(GeV) 


g 


VcilS) 


CC 


2.9779 


1 


J/^{1S) 


CC 


3.0970 


3 


Xco(lP) 


CC 


3.4152 


1 


Xci(lP) 


CC 


3.5106 


3 


hcilP) 


CC 


3.526 


3 


Xc2(lP) 


CC 


3.5563 


5 


Vci2S) 


CC 


3.638 


1 


ij{2S) 


CC 


3.686 


3 




CC 


3.770 


3 




CC 


3.929 


5 




CC 


4.040 


3 




CC 


4.159 


3 




CC 


4.415 


3 


Be 


he 


6.27 


1 


' — 'CC 


ccq 


3.527 


4 




CCS 


(3.660) 


2 



similar, but suppression is not as pronounced. For 7^ = 7^'' suppression the yield of hidden 
charm particles is always smaller than equilibrium value. This suppression occurs due to 
competition with the yield of strange-heavy mesons, and also, when 7^ > 1, with heavy 
baryons with two light quarks. The enhanced yield of D,Ds and heavy baryons in effect 
depletes the pool of available charmed quark pairs, and fewer hidden charm cc mesons are 
formed. For particles with two heavy quarks the effect is larger than for hadrons with one 
heavy quark and light quark(s). 

In figure [3^22] we compare the J/^ yield to the chemical equilibrium yield "^/J/^eq, 
as a function of 'y'^/'y^, each line is at a fixed value 7^. This ratio is: 



J/* 



hid 



N, 



hid eg 



/ceq 



(3.106) 



J/^'/J/U'oq always decreases when 75/7^ increases. For 7^ = 7cr J/^/J/^cq is smaller 
than unity even when 7^ — > 0, because of large phase space occupancy of light quarks. 

cq > 1 for small 7^ and small 7^/ 7^ . This ratio decreases with 7^/ 7^ grow. This 
J/^ I J/^eq ratio behaviour is similar to experimental results for SPS enerdies shown in fig- 
ure ?? [36]. At SPS energies the chemica potential the chemical potential has influence on 
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Figure 3.22: (Color on line) Ratio J/^ / J/'^eq = ll/lleq ^ function of 7^/7^ at fixed value 
of 7^ and if required, entropy conservation. Shown are: T = 200 MeV at 7^ = 0.83 (dot-dash 

T = 170 MeV at 7^ = l(upper dashed 



line) and at 7^ = 7^** = 1.42 (lower solid line (purple) ); 
line) , at 7g = 1.15, (upper solid line (red)), and at 7^ 
T = 140 MeV, 7<j = 1.6 



1.51, (lower dashed line); and 



result. However it have not to change the effect qualitively. The reason of J/U' suppression 
in this case can be the same as in our model. 

Considering the product of Jf^ and yields normalized by A^"^ we eliminate nearly 
all the uncertainty about the yield of charm and/or hadronization volume. However, we 
tacitly assume that both Jf^ and hadronize at the same temperature. In figure [3l23] we 
show J /'^(p/N'^ as function of 'js/'^q- There is considerable difference to the ratio considered 
in figure [3T3l We see mainly dependence on 7^/75. As before, see section 13.2.21 JM/ is the 
sum of all states cc from table [23] that can decay to J/U'. We show results for T = 200 MeV 
(solid lines), T = 170 MeV (dashed line) and T = 140 MeV (dash-dot line). The 7,, for 
each T, is fixed by entropy conservation condition during hadronization (figure 12.61) (thick 
lines) or by 7g = 7^^ (thin lines). For T = 140 MeV these lines coincide. T = 170 MeV, 
7g = 1 case is also shown (solid line with dot markers). The s/S values, which correspond 
to given 7s/7g ratio can be found in figure [2781 Figure [3723] shows that despite the yield 
(f) / {dV / dy) increasing as {■js/lqY, J/^4>/N^ is increasing as 7s/7g considering compensation 
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's 'q 

Figure 3.23: (Color on line) J/^(f)/N^ states yields as a function of js/lq ratio for T = 200 
MeV, 5*3 = (solid line) and 7, = 7^^ (solid thin line), T = 170 MeV: 5^ = (dashed hne), 
7g = 7^^ (thin dashed line) and 7^ = 1 (solid line with dot marker); for T = 140 MeV, = 
(dash-dot line) 

effects. 

A similar situation, as in figure [3^221 for hidden charm, arises for the Be meson yield, 
see figure 13211 where Bc/NcNi, ratio is shown as a function of hadronization temperature 
T, for the same strangeness yield cases as discussed for the hidden charm meson yield. 
Despite suppression in strangeness rich environment, the meson yield continues to be 
larger than the yield of B^ produced in single NN collisions, where the scale yield is at the 
level of ~ 10~^, see cross sections for bb and Be production in [56] and in [70], respectively. 

In figure 13.251 we show scaled yields of ccq and ccs baryons as a function of 
temperature. Upper panel shows aside of the equilibrium case (dashed lines) the yields 
for s/S = 0.04 with = S'^ (sohd lines) and with 7^ = 7^*" (dash-dot line) for dV/dy = 
800fm-3 for T = 200 MeV. Lower panel is for dV/dy = 600 fm'^ and s/S = 0.03. For the 
ccq baryons the chemical nonequilibrium suppression effect is similar to what we saw for cc 
and Be mesons. Equilibrium yield is much larger than non-equilibrium for T < 230 MeV 
when s/S = 0.04 and = S'^, and for T < 190 MeV when s/S = 0.03 and S" = S^. In 
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Figure 3.24: (Color on line) Be mesons yields as function of T for chemical equilibrium case 
with dV/dy = 600 fm^^ for T = 200 MeV (the upper panel, dashed line), for s/S = 0.03 with 
dV/dy = 600 fm^^ for T = 200 MeV (the upper panel, solid line), for chemical equilibrium 
case with dV/dy = 800 fm^^ for T = 200 MeV (the lower panel, dashed line) for s/S = 0.04 
dV/dy = 800 for T = 200 MeV (lower panel, sohd line) 
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s/S=0.03, dV/dy=600 fm'^ for T=20q MeV" 

140 160 180 200 220 240 260 

T [MeV] 

Figure 3.25: (Color on line) ccq/N"^ (upper lines in each panel) and ccs/N^ (lower lines in each 
panel) baryon yields as a function of T. Upper panel: chemical equilibrium case with dV/dy = 800 
fm-3 for T = 200 MeV(dashed line), s/S = 0.04 with dV/dy = 800 /m"^ for T = 200 MeV: 
= S'^ (solid line) and = 7^'' (dash-dot line); and lower panel: chemical equilibrium case 
with dV/dy = 600 fm'^ for T = 200 MeV (dashed hue), and s/S = 0.03 S"3 = (solid line) 
and 7"^ = 7g'' (dash-dot line). 

case 7g = 7^^, the yield of ccq is always smaller than equilibrium. The yield of ccs baryons 
has similar suppression, but it becomes larger than equilibrium for smaller temperatures 
and yield enhancement for higher T is larger for = then in case of ccq because of 
large number of strange quarks. 

In the figure [3^26] we show ratios ccq/Jf^ (upper panel) and ccs/J/^ (lower panel) 
as a function of hadronization temperature. These ratios do not depend on dV/dy. 
ccq/J/^ oc 7g does not depend on s/S. For ccq/ J/^ ratio we show three cases: chemical 
equilibrium 7s = 7g = 1 (dashed line), = S"^ (solid line) and 7^ = 7^'' (dash-dot line). 
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Figure 3.26: (Color on line) ccqj Jj^ (upper panel) and ccs/ J/^ (lower panel) ratios as a function 
of T. Upper panel: chemical equilibrium case (dashed line), = 5*^ (solid line) and 7^ = 7^'' 
(dashed-dot line); and lower panel: chemical equilibrium case with (dashed line), s/S = 0.04: 
5*^ = (solid line with dot marker) and 7^ = 7^*" (thin dash-dot line); s/S = 0.03 (solid line) 
and 7g = 7^'' (thin dash-dot line). 

For ccs/J/^ {ccq/J/"^ oc 7^) we show chemical equilibrium case (dashed line), s/S = 0.04: 
= S*^ (solid line with point marker) and 7^ = 7^*^ (thin dash-dot line); s/S = 0.03: 
S^ = S'^ (solid line) and 7^ = 7^'' (thin dash-dot line). The overall all yields of double 
charmed (strange and non-strange) baryons and anti-baryons is clearly larger than the 
yield of J/^. 
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3.4 Conclusions 

We have considered here in some detail the abundances of heavy flavor hadrons within the 
statistical hadronization model. While we compare the yields to the expectations based 
on chemical equilibrium yields of light and strange quark pairs, we present results based 
on the hypothesis that the QGP entropy and QGP flavor yields determine the values of 

phase space occupancy i = q,s, c, b, which are of direct interest in study of the heavy 
hadron yields. 

For highest energy heavy ion collisions the range of values discussed in literature 
is 1 < 7^ < 1.65 and 0.7 < jf/jf < 1.5. However 7^ and 7^ values which are much 
larger than unity arise. This is due to the need to describe the large primary parton based 
production, and considering that the chemical equilibrium yields are suppressed by the 
factor exp(— m/T). 

Our work is based on the grand canonical treatment of phase space. This approach 
is valid for charm hadron production at LHC, since the canonical corrections, as we have 
discussed, are not material. On the other hand, even at LHC the much smaller yields of 
bottom heavy hadrons are subject to canonical suppression. The value of the parameter 
7^ obtained at a fixed bottom yield N^, using either the canonical, or the grand canonical 
methods, are different, see e.g. Eq. (15) in [71]. Namely, to obtain a given yield Nf, in 
canonical approach, a greater value of 7^ is needed in order to compensate the canonical 
suppression effect. However, for any individual single-6 hadron, the relative yields, i.g. 
B/Bg do not depend on 7^ and thus such ratios are not influenced by canonical suppression. 
Moreover, as long as the yield of single-6 hadrons dominates the total bottom yield: Nf, ~ 
B + Bs + Ab + . . ., also the Ni, scaled yields of hadrons comprising one 6-quark i.e. ratios 
such as B/Nb, Bg/Nb, Bc/Nb, etc, are not sensitive to the value of 7^ and can be obtained 
within either the canonical, or grand canonical method. On the other hand for bb mesons 
and multi-6 baryons the canonical effects should be considered. Study of the yields of these 
particles is thus postponed. 

We address here in particular how the yields of heavy hadrons are influenced by 
if/l^ 7^ 1 and 7, 7^ 1. The actual values of 7^/7^ we use are related to the strangeness 
per entropy yield s/S established in the QGP phase. Because the final value s/S is 
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established well before hadronization, and the properties of the hadron phase space are 
well understood, the resulting 7^/7^ are well defined and turn out to be quite different 
from unity in the range of temperatures in which we expect particle freeze-out to occur. 
We consider in some detail the effect of QGP hadronization on the values of 7^ and 7^. 

One of first results we present (figure [3TT3|l allows a test of the statistical hadroniza- 
tion model for heavy flavor: we show that the yield ratio cc ss/ {cs cs) is nearly independent 
of temperature and it is also nearly constant when the is allowed to freeze-out later (fig- 
ure [3]T1]), provided that the condition of production is at the same value of strangeness 
per entropy s/S. 

We studied in depth how the (relative) yields of strange and non-strange charmed 
mesons vary with strangeness content. For a chemically equilibrated QGP source, there 
is considerable shift of the yield from non-strange D to the strange Dg for s/S = 0.04 
expected at LHC The expected fractional yield Dg/Nc — Bg/Ni, ^ 0.2 when one assumes 
1^ — Iq^ — '^1 expected enhancement of the strange heavy mesons is at the level of 
30% when s/S = 0.04, and greater when greater strangeness yield is available. 

As the result we find a relative suppression of the multi-heavy hadrons, except 
when they contain strangeness. This suppression depends on both factors 7^ and gammaq. 
When phase space occupancy of light and strange quark is relatively high the probability 
for charm quarks to make hadrons with strange quarks increases and probability to find 
the second charm quark among light and strange quarks decreases. Therefor the cc yield 
suppression increases when js/lq ratio increases for constant 7^. This result is qualitatively 
in agreement with experimental results obtained for SPS energies [36]. 

On the other hand, the yield of cc/N^ ^ 210~^ is found to be almost independent 
on hadronization temperature when entropy at hadronization is conserved. That is be- 
cause for larger T 7, decreases. The suppression effect decreasees, compared to SHM and 
become even negative for T > 200 MeV, resulting to the cc yield almost independent on 
temperature. We don't know exactly equation of state in QGP and so the value of 7^ which 
is needed to conserve the entropy may be different. If 7^ is larger for higher temperatures, 
suppression of cc is larger for a fixed s/S. The same result is found for i?c ~ 5 — 6 lO^^A'^cA'";,, 
that yield remains considerably larger (by a factor 10 — 100) compared to the scaled yield 
in single nucleon nucleon collisions. 
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We have shown that the study of heavy flavor hadrons will provide important 
information about the nature and properties of the QGP hadronization. The yield of Bc(6c) 
mesons remains enhanced while the hidden charm cc states encounter another suppression 
mechanism, compensating for the greatly enhanced production due to large charm yield at 
LHC. 
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CHAPTER 4 

CHEMICAL EQUILIBRATION INVOLVING DECAYING PARTICLES AT 

FINITE TEMPERATURES 



4.1 Introduction, motivation, overview 

In this chapter we consider relativistic master population equation and equations for reac- 
tion rates similar to considered in section ll.3[ but extended to the case of two — body ^ 
one — body reaction ll.48[ To our surprise, we reahzed that such two — body ^ one — body 
reaction has so far NOT been addressed in the relativistic context in literature. This study 
begun with the question at which temperature in the expanding early Universe the reaction 
7 + 7 ^ vr° 'freezes' out, that is the vr^ decay overwhelms the production rate and the yield 
falls out from chemical equilibrium yield. This reaction has the lowest threshold, one pion 
mass, with two thermal particles available to reach it. Thus this reaction should be still 
operational at a relatively low temperature when all other hadron production reactions 
cease to be effective. 

In chapter [6] we will also consider this reaction in e"'"e~7 plasma, created by high 
intensity laser pulse. Aside of cosmological implications, insights from this study are clearly 
of relevance to the general understanding of quark gluon plasma and hadron gas evolution 
in relativistic heavy ion collision. For example this study allows us to consider the chemical 
yields arising in reactions such as p ^ tttt, tt^ ^ 77, A ^ Ntt, and so on [4]., chapter [5l 

We recall here that the reaction 11.481 considered in the rest frame of the decaying 
particle implies the constraint mi + 1712 < since 

■"^3 = (Pi + PiY = {"mi + 7722)^ + 2{EiE2 - mim2 - V\ ■ V2) > ("^1 + "^■2)^, (4.107) 

considering that the condition E\E\ > {17111712 + Pi ■ 'P2Y implies (mi|p2| — ""^21^1)^ > 
2mim2f)i ■ P2 — 2mim2|pi| \p2\ which is always true. The equality sign corresponds to 
the case that mi + 1712 = m-^ for which the reaction rate vanishes by virtue of vanishing 
phase space. This text book exercise shows that the raction Eq. (??) is possible, has a 



91 



'good' phase space size, and it invites to evaluate the rates of the processes of interest in 
the rest frame of particle '3', boosting, as appropriate, from/to laboratory frame. To do 
this effectively we will need to formulate the master population equations in explicitely 
covariant fashion. 

This constraint Eq. fl4.107p forbids many reactions. For example, the hydrogen 
formation p + e — >H is forbidden since for a bound state mn < mp + rrie. Thus there 
must be a second particle in the final state, the electron capture involves either a radiative 
emission, p + e — >H+7 or a surface/third atom, which picks the recoil momentum. The 
situation would be different if there were 'resonant' intermediate state of relative long 
lifespan with energy above ionization threshold. Such 'doorway' resonances are available 
in many important physical processes, including e.g. the d + t ^ a + n fusion. 

The general kinetic master equation approach to reactions of type ll.22l for the yield 
(chemical) equilibration in nuclear and particle physics has been studied frequently in the 
context of heavy ion reactions [16]. However, the simpler situation was not considered 
in this framework, and the adaptation is not trivial given novel quantum and relativistic 
effects involving particle decay. 

At temperatures T~m(^ = c = A; = l) the particle number present follows 
rapidly the relativistic statistical phase space. Due to the conservation of energy and 
momentum, reaction Eq. (11.481) is subject to a particular kinematical constraint. 

In the present work we present for the first time the dynamics of reaction Eq. 
(I1.48P in a microscopic description of particle production and the associated decay reac- 
tions within the frame work of kinematical master equation obtained from the Boltzmann- 
Uehling-Uhlenbeck (BUU) equation under thermal bath [72]. 

We will ask questions such as "Will decaying particles reach chemical equilibrium, 
and if so, how fast? Does the presence of background thermal particles stimulate or slow 
down reaction rate?" . In the foillowing section 14.21 we write down the kinetic equations for 
time evolution of number density n of decaying particle and equations for invariant rates. 
We show that the time variation of density of particle 3 is 



where dW3^i2/dVdt is the decay rate of particle 3 and Tj is fugacity for the particle i. 




(4.108) 
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Here the number density of particle i in thermal (kinetic) but not necessarily in the 
chemical equilibrium is given by: 



fb/f, defind by Eq. \1.27\ is the covariant form of the usual Bose or Fermi distribution 
function defined in the rest frame of the thermal bath, and describes the corrresponding 
quantity in a general reference frame where the thermal bath has the relative velocity 
defined by u'^, see Eq J1.4[ Note that the distribution function / is a Lorentz scalar but 
spatial density rii is not. 

The particle 3 attains the chemical equilibrium when the following condition among 
fugacities is satisfied: 



This, as expected, is equivalent to the Gibbs condition for the chemical equilibrium. In 
section H73] we evaluate invariant rate using decay time in vacuum in rest frame of decaying 
particle, and discuss the behavior of the average decay rate of an unstable particle in the 
presence of thermal bath. In section 14.41 we apply our formalism to three examples: 

a) relaxation time of formation p meson through vr + vr ^ p in a baryon-free hot hadronic 
gas, where mesons are considered in thermal and chemical equilibrium; 

b) the decay and production relaxation time of resonance S(1385) in reaction E(1385) ^ 
Att in hot hadronic gas; 

c) 7r° equilibration the reaction 7 + 7 ^ tt*^ in thermal e+e~7 plasma or early universe 
However, one should note that at sufficiently low temperatures the local density of 7r° is 
too low to apply the methods of statistical physics. 

4.2 Kinetic equations for decaying particles 

4.2.1 Decaying particle density evolution equation 

Consider an unstable particle, say 3, which decays uniquely into other two particles. 




(4.109) 



T1T2 = 



T3. 



(4.110) 



3^1 + 2 



(4.111) 
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in the vacuum. In a dense and high temperature thermal ambient phase particles 2 and 3 
are present, and the inverse reaction: 

1 + 2^3 (4.112) 

can occur to produce the particle 3. If we assume that the abundance of particle 3 changes 
solely by thermal production of particle 3 by particles 1 and 23 and its decay via Eqs. ( 14.112t 
14.1111) . then we write the time variation of the number density as: 

dt ~ dVdt dVdt ' I • J 

where dWi2-^^/ dVdt is the production rate per unit volume of particle 3 via Eq. fl4.112p 
and dWs-^12 / dV dt is the decay rate of particle 3 per unit volume. 

In a normal situation, the abundance of particles 1 and 2 are determined by the 
other processes which produces these particles. For example, consider the reaction p ^ tttt 
in dense hot matter formed the heavy ion collisions. Then pions can be easily created 
by inelastic collisions of other hadrons and thus in principle we have to deal with multi- 
component systems. However, there exists special situation where the total abundances of 
particles 1 and 2 are initially determined, and in the following the time variation of number 
densities of particles 1 or 1 is established by the above reactions In such cases we have: 

dni^2 ^ dWs^u _ dWi2^3 ^ 
dt ~ dVdt dVdt ' I • J 

In the following, we assume that the system is spatially homogeneous and all of the 
particles are in thermal equilibrium. Furthermore, we consider that the interaction time 
among particles is short enough so that all the dynamical information can be obtained from 
the single particle distribution function / (p) for each particle. In a thermal equilibrium, 
this function is specified completely by 2 parameters, T the temperature and T the fugacity. 
In this paper, we assume that the thermal back-ground is inert, so that we keep T constant, 
but the fugacity T changes in time through the chemical reactions so does the density of 
each component of the gas. 

The thermal production rate dWi2-^3/dVdt and the decay rate of the particle 3 
under the thermal background dW^^u/dVdt can then be expressed using these distribution 
functions for each of particles involved in the reaction. 
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4.2.2 Decay and production rates 

According to the boson or fermion nature of the particle 1, we have to consider different 
cases. If the particle 1 is boson, then there are two different cases of the decay and 
production mode: 

bosons < — > bosoni + boson2, (4.115) 



bosons < y fermioui + fermion2. (4.116) 

On the other hand, if the particle 3 is fermion it should decay into a boson and a fermion: 

fermious < > bosoni + fermion2. (4.117) 



The Lorentz invariant transition probability per unit time and unit volume corresponding 
to the process fl4.112p can be expressed as 

X (27r)^ 6' (pi +p,-p,)^Yl I |M| ps)f (1 ± /,,/(T3, P3))(,4. 118) 

where / = 1 for the reaction of two indistinguishable particles 1 and 2, and / = if 1 and 
2 are distinguishable. The factor l/{gig2) and the summation are due to averaging over 
all initial spin states. The last factor accounts for the enhancement or hindrance of the 
final state phase due to the quantum statistical effect, as is introduced first by Uehling and 
Uhlenbeck [72]. The sign is for the case when the particle 3 is boson and ' when it 
isa fermion. It is clear that Eq. f l4.118p is manifestly Lorentz invariant and it can be used 
in any frame of reference. 

Now we write in the same way the decay rate of the process (14.21) . per unit volume 

we have: 

dW3^l2 93 f d^P3 X /or „ N /" ^^Pl f d^P2 /o 



dVdt (2vr)3 j 2eM'^'^^'^ j 2EA!^^ j 2^(2^^'") ^ bi +P2 -Pa) x 

E 1(^3 \M\pm)? (1 ± A/(Ti,Pi)) (1 ± A,/(T2,P2)) (4.119) 



^31 + 

spm 



Here, the one particle state is normalized as 

{p'\p) = 2/ (27r)^ 53 (p' _ p-) . (4.120) 
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The quantum statistical effects on the final state affects the decay rate, compared to the 
free-space case. The decay process is simulated by the presence of thermal background 
particles 1 and 2 if they are both mesons, and hindered if they are both fermions. 

Note that the pure thermal production rate dWi2^3/dVdt of the particle 3 is 
related to its decay rate dW3^i2/dVdt through the time-reversal relation of the transition 
matrix element which can be shown in the following. Using Eq. fll.35p . we can rewrite the 
Eq. flCTSl) as 

^ydr = TTl^' j 2E;{^^ j 2E^^ j 

^|(p3|M|pip2)|'A(Ti,Pi)A/(T2,P2)A/(T3,P3)exp(n-p3/T); (4.121) 

spin 

and for the Eq. fl4.119p . using energy-momentum conservation Pi+P2= Ps we obtain 

dW3^U 1 f d^pi f d^p2 f (fps ^4 .4. , ^^ 



dVdt l + I ' 'J 2Ei{2nyJ 2E2{2nyJ 2E3 (27r)^ 
J2\{PiP2\M\p3)fMT,,p,)f,jiT2,P2)hAr3,P3)expiu-p3/T); (4.122) 

spin 

Using the time reversal symmetry of the transition matrix element, 

I (P3 \M\ P1P2) f = I (P1P2 \M\ P3) f , (4.123) 

we find: 

dWi2^3 1 ^1^3^12 ' 



dVdt T1T2 dVdt T3 " ^^-^^^^ 

which is the detailed balance relation for the process of formation and decay of unstable 
particle. Therefore chemical equilibrium T1T2 = T3 corresponds to equal decay and 
production rates as we expected. Using this relation, Eq. fl4.113p can be written in the 
form of Eq. (14.1081) . 

Given a thermal bath with a fixed temperature T, we wish that the change of 
number density is related directly to the change of fugacity. This is achieved by defining 
the decay time by 

ra - (4.125) 



where: 
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Therefore, from Eg. (14. 1081) the time derivative of the fugacity of the particle 3 is: 

t3 = (TiT2-T3)-. (4.127) 

For the case where the abundances of 1 and 2 are determined only from the reac- 
tions (3 ^ 1 + 2), then analogous expressions for particles 1, 2 are obtained by introducing 
Ti for each particles as 

n = (4.128) 

4.3 Calculations of invariant decay (production) rate 
4.3.1 General case 

The vacuum decay width of particle 3 in its own rest frame is found in textbooks. In our 
notation: 



J Ztj2 277) (73 ~~ 



To 2m3l + /y 2E^{2t,Y J 2E2{2r, 



2m3^3 4(/+l)(2 
1 1 p 

^msgs 4(/ + 1) nrris 



' ^ ' spm 



J2\{P,-p\M\m3)f (4.129) 



spm 



Here p = Pi = P2 and E12 = yp^ + mfg are the magnitude of the momentum and, 
respectively, the energy, of particles 1 and 2 in the rest frame of the particle 3. From 
energy conservation: 

Ei,2 = , P =E,^,-m,^2 = ^ 2^ + 4mi " ^^'^^^^ 

We denote by the decay rate of the particle 3 in the rest frame of the thermal 
bath, E^ and p^ are the corresponding energy and the momentum. The thermal decay rate 
per unit volume dW3^i+2/dVdt should then be the average (over the inverse of this life 
time) in the thermal bath frame: 

dW3^l+2 _ 93 

dVdt ~ {27ff J -3 '3 

where E^r^m^ is the decay time of the particle 3 with moment ^3. 



[ C^'P3A/(T3,P3)^^, (4.131) 
J E3 
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Comparing this expression Eg. (14. 1311) with the complete Eq. fl4.122p . we conclude 
that the in medium, at finite temperature T, decay rate r!^ is given by: 

-El(!'iP2l''''IP!>>l'A./(Ti.Pi)A./(^2'K)Tr'T2-'exp(t,.ps/r), (4.132) 

spin 

which is a Lorentz invariant form. We note that u ■ = denotes the energy of the 
particle 3 in the rest frame of the bath. 

Using the vacuum rest-frame decay time, Eg. (14. 129 p . we find that Eq. (l4.132p takes 
the form: 

1 1 e^3/T 

-7 = -^^iPs)- (4.133) 

To 2 



The function $(^3) is: 



dC -1 , , , . (4.134) 



with 



E\ E^ E2 E^ PPs 

ai = — , 02 = —, b= — — and C = cos 6^ = cos(p2 A pi). (4.135) 

m^I m^I m^I 

The integral ^(ps) can be evaluated analitically. The integrant in this equation 
{( = x) is even, therefore 

$(»^) = 2 / dx 1 , , , 1 , , , = / dx —-r Tzr^ 

(4.136) 

Introducing y = e^^ this integral can be written as 



b A '^(T7Vi±y)(y±T2e--2) 



b{e-i+-^ - T1T2) ii ^ V (Tr'e'^i ±1/) {y± T2e-''^) 



M i;-y,"::2g- 1 . (4.137, 
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The result is 



We note that in the non-relativistic hmit (ma ^ T,p^), this quantity tends to 



Finally, the average particle 3 decay rate per unit volume in a thermally equilibrated system 
is given by 

dVdt (2x2) T„ i„ Tjiea/i-il ^ ' 

4.3.2 Decay and production rates in Boltzmann limit 

The equations become much simpler in case of Boltzmann limit when we can omit unity 
in distributions Eq. fll.27l) . This is possible when 

that is, when Tj ^ 1 or T ^ mi/2. The condition T <^ mi/2 comes from fact that the 
minimal energy of lighter particles is mi/2 in the particle 1 rest frame. In this limit the 
decay time in the particle 1 rest frame from Eq. (14. 1321) r' tq so that from Eq. (14. 1251) we 
have for the average decay rate r in the reference frame (the rest frame of the bath) as 

^3 ^ ^O rJ\; ,^ (4.142) 

Jo p dpe^-^'^rris/hs 

- -ro^^P^. (4.143) 



(1/7) "i^i(mi/T)' 

As we see in Eq. (l4.13T|) . the average decay time in lab frame is proportional to the 
(inverse) average of Lorentz factor of particle 3. We will discuss this effect next in quan- 
titative manner, and note that the ratio of Tg to Tq is shown in figure 14.271 as dash-dot 
line. For T ^ ms this ratio goes to unity because the Lorentz factor becomes 1. For large 
T, the rate increases because of the larger average energy of particle 3 or equivalently the 
larger average Lorentz factor. Therefore, if we have small number of all particles (Tj <^ 1) 
so that Eq. (14. 1411) is yet valid for T > m^, the average particle life time increases with T 
due to relativistic effects. 



99 



4.4 Examples 

4.4.1 Production of p mesons via p ^ tttt process 

Here we consider example of p meson thermal decay and production: 

p°^7r+ + 7r-, (4.144) 
p^^7r^ + 7r°. (4.145) 

We consider the pions to be in chemical equilibrium with chemical potential /i^r = 0. In this 
case all particles are bosons and mi = m2 and in integral (14.1341) we have Ei = E2 = mp/2 
in p rest frame. Integrant in $(p) function is symmetric function. Then we can write 







where 



2T 



2T 

The integral (I4.146P can be evaluated in this case as 



(4.147) 
(4.148) 



Then we substitute $ into Eg. (14. 131 1) and using Eg. (14. 1241) we can calculate p decay and 
production rates. To calculate r we use definition (14.1251) . 

In figure 14.271 we presents p decay time in lab frame as a function of temperature 
T for Tp = = 1, solid line, for Tp = T^r = 1.5, solid with dot marker, Tp = T^r = 0.1, 
dashed line and dash-dot line is for Boltzmann limit Eg. ( 14.1431) . We consider range of 
temperatures between 50 and 300 MeV which includes guark gluon plasma hadronization 
temperature (~ 140 -180 MeV). We show case Tp = = 0.1 to check transition to 
Boltzmann limit. We can see that for this case result is close to Boltzmann approximation 
for our range of T as it is expected. In case Tp = = 1 we have chemical eguilibrium. In 
this case and for Tp = T^r = 1.5 for small T <^ mp/2 we have ratio t/tq near Boltzmann 
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T [MeV] 

Figure 4.27: The ratio t^/tq as a function of temperature T in the reaction p ^ vrvr. The dashed 
hne is for Boltzmann hmit showing just time dilation. Nearly this limit arises (dot-dashed line) 
for Tp = = 0.1. Solid lines are for Tp = T^r = 1 (top, red) and Tp = T^r = 1.5 (bottom). 

limit, near unity. For such small T, when Boltzmann limit is applied, decay time r doesn't 
depend on T. When T increases quantum effects take place then r begin to decrease with 
increase T. The larger T the faster r decreases with temperature. 

The case with T^r = Tp ^ 1.5 can take place after quark qluon plasma hadroniza- 
tion. Light hadrons multiplicities has to be above chemical equilibrium for hadronization 
temperature smaller than 180 MeV to conserve entropy during hadronization. 

4.4.2 Baryon resonance (S(1385)) lifespan calculations in dense hadronic gas 

In this subsection we consider the effect of oversuturated pion component in hadronic gas 
and the effect of the motion of the decaying resonance with respect to the thermal rest 
frame on its lifespan and then also on resonance production relaxation time, considering 
example S(1385) ^ Att. 

For the temperatures of interest (hadronization of QGP and below) m/^ and my; >> 
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Figure 4.28: (color on line) The Bose enhancement factor 1 + fn{Ei) in 11(1385) rest frame as 
a function of light quark fugacity for the reaction S(1385) Avr at T = 140 MeV (blue, solid 
line), at 160 MeV (green, dash-dot line) and 180 MeV (red, dashed line). The dots show the 
initial value of fugacities for the three possible hadronization cases. 



T. In this case with sufficient accuracy we can rewrite function $(^3) as 



1 , (e'^^+^-T^) 
In ■ 



(4.150) 



Here fugacities for A and vr correspond to those for particles 1 and 2, respectively. There 
are no significant medium effects upon decay rate of S(1385) and A resonances. However 
the pions have energy E2 = 250 MeV (Eq. fl4.130p ) in the S rest frame and the Bose 
enhancement effect is possible in the oversaturated hadronic gas after QGP hadronization. 

For the low temperatures considered here we can assume that S resonances almost 
do not move. Thus the enhancement effect in the thermal bath frame is close to the 
enhancement in the S(1385) rest frame. The decay rate increases by Bose enhancement 
factor 1 + /,r (here /,r = ^(-£^2) ^))- figure 1^1281 we show Bose enhancement factor 
as a function of light quark fugacity 7^ for temperature Tq = 140 MeV (blue, solid line), 
To = 160 MeV (green, dash-dot line), Tq = 180 MeV (red, dashed line). The large dots show 
Bose enhancement factor for our initial 7^ determined from entropy conservation in fast 
hadronization. The fugacity 7^ = 1.6 is close to maximum expected value at Tq = 140 MeV. 
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80 100 120 140 160 180 
T [MeV] 

Figure 4.29: (color on line) The ratio of the in medium lifespan T3 with the vacuum lifespan 
To as a function of temperature T for the reaction 5](1385) <-> Avr. The dashed (red) line is for 
hadronization at To = 180 MeV, jq = 1.0; the dot-dashed line (green) for hadronization at 160 
MeV, 7q = 1.27; solid line (blue) is for hadronization at 140 MeV and 7^ = 1.6. 

The maximum fugacities for each temperature correspond to Bose - Einstein singularity. 
The Bose enhancement effect is largest for maximum 7^ and it diminishes for small 7^. At 
fixed entropy the greatest enhancement is for smallest ambient temperature, see the dot 
on solid line in figure I4.28[ 

In figure H^29] we show the corresponding decrease in the lifespan, the ratio t^/tq 
as a function of temperature T in the reaction S(1385) ^ Avr. We consider temperature 
range from corresponding hadronization temperature until T = 70 MeV. We assumed, that 
T,r is a constant. Fugacities of heavy resonances do not influence the result. The lowest 
t-^/tq ratio is for 7^ = 1.6 at Tq = 140 MeV when we have maximum value of 7^ for given 
temperature. If we compare this value of t-^/tq = 0.65 with inverse Bose enhancement 
factor 1/(1 + /^(£'2,T)) = 0.54 for this T and 7^ (see figure 1^^281) we see that these values 
are near to each other (difference is about 20% ) as expected for » T. For smaller 
T, 7g decay time goes to its vacuum value. 

The same calculations are applicable for heavier E*. When the difference of mass 
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of the initial and final state resonance decreases, the Bose enhancement effect increases, 
since it involves small momenta. The largest effect is for reaction E(1670) ^ A(1520) + vr. 
On the other hand, for the reactions which satisfy condition ms — (mi + 7712) > 300 MeV 
the enhancement effect becomes very small. 

4.4.3 Thermal Production of 7r° 

As mentioned in the Introduction, it is interesting to examine the mean life time of vr" in 
the end of hadronic gas stage of the universe where the temperature becomes several tens of 
MeV. Then the 7r° production in two photons fusion. Eg. (11.691) determines the abundance 
of 7r°. 

The difference with previous example is that the photons are massless and they 
are in chemical equilibrium (Ti = T2 = 1). Then we can rewrite function (14.1381) as 

with 



a = (4.152) 
b = (4.153) 

Again we use Eq. (l4.13ip and (I4.124p we can calculate ttq decay and production rates. To 
calculate r we use definition (14.1251) . 

In figure 16.2641 we show ratio of vr'' decay time in the presence of thermal particles 
to the decay time in vacuum in ttq rest frame: r^o/r^o- In this figure the large range 
of temperature is shown 10 — 10^ MeV. For T^o = 1 the ratio r^o/r°o the temperature 
dependence is similar to that for p decay, considered in previous chapter. It increases at 
first until relativistic effects become noticeable. Then, after T goes slowly down with 
temperature, when quantum in-medium effect becomes important. Range of change of r is 
not large for this large range of temperature. The smallest Tj^o/t^^ is about 0.6 at T = 10'^ 
and the maximal value of this ratio is about 1.2. 

The cases when T^^ 7^ 1 are different from those considered for p decay because 
photons are stay in chemical equilibrium. When T^^ < 1, T^^ = 0.5 and Tttq = 0.1 (purple 
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Figure 4.30: The ratio t/tq for vr*^ decay/production as a function of temperature T. Dotted, 
blue line is for T^o = 3, dash-dot, green line is for T^o = 2, dashed, red line is for T^o = 1.5, 
solid, turquoise line is for T^o = 1, purple solid line with triangle markers is for T^o = 0.5, brown 
solid line with dot markers is for T^o = 0.1. 

line with dots and brown line with tringles, respectively), the r also begins to decrease 
slowly after T ^ 20 MeV because of quantum effect from photons distribution. Then for 
T > 300 it is slowly increases because relativistic effects becomes slightly dominant. 

When Tttq > 1, T^ro = 1-5 (red, dashed line), T^^ = 2.0 (green, dash-dot line) and 
T^g = 3.0 (blue, dotted hne), there is Bose - Eistein critical point when 

T = m^o/log(A). (4.154) 

1 d\ is increasing faster near this critical point than 7r° production rate. Decay relax- 
ation time r^o goes sharply up, diverges, near critical point. 

In figure I473T] we show r^o for more realistic temperatures. This temperature range 
can be interesting for early universe evolution and for e"'"e~7 plasma created by laser pulse. 
It turns out that there are both relativistic and quantum effects which contribute and 
they (nearly) cancel at this range of temperature. The relativistic effect arises because 
r^o in Eq. (16.2641) is in lab frame while the known r°o is in the pion rest frame. In the 
relativistic Boltzmann limit the correction is obtained considering the related time dilation 
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Figure 4.31: The ratios t^o/t^^ as functions of temperature T for relativistic Boltzmann limit 
(blue, dashed line) and for quantum distribution in chemical equilibrium, T^r = = 1 (green, 
solid line). 

effect Eq. fl4.143p . We find that this effect imphes that r^o in the lab frame increases with 
temperature. This effect is shown by dashed (blue) line in figure 14.311 Furthermore, with 
increasing temperature quantum distribution functions for photons and for the produced 
particle need to be considered. This leads to the result shown as solid line (green) in figure 
14.311 Thus in general r^o > r°o, by up to 14%. 

4.5 Conclusions and Discussion 

In this chapter, we examined in detail the kinetic master equation for the process involving 
formation of an unstable particle through the reaction Eq. (11.481) in a relativistically co- 
variant fashion. Assuming that all the particles in the process are in thermal equilibrium, 
we calculate the thermal averaged decay and formation rate of the unstable particle based 
on the BUU equation. Using the time reversal symmetry, we show that the time evolution 
of the density of the unstable particle as Eq. (14. 1081) . Therefore in chemical equilibrium 
particles fugacities are connected by Eq. (I4.110p as expected. We have explicit the thermal 
decay rate of unstable particle, obtaining Eq. (l4.140p . which is our principal result. 
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Using the formalism developed above, we examined the general properties of the 
thermal particle decay/production rate. We see that for T <^mi where the Boltzmann limit 
can be applied, the decay width is reduced to Ti/tq and production width is T2T3/ro. For 
larger values of T but Tj ^ 1 so that the Boltzmann approximation is valid, then decay 
width and production width tend simply to Ti/r and T2T3/T, respectively, where r is 
essentially proportional to average Lorentz factor and doesn't depend on Tj. When some 
of rrii/T and Tj are about unity or larger we see dependence of r on Tj. 

We applied our formalism to 3 examples, p ^ vr + tt, S(1385) ^ and 7r° ^ 7 + 7. 
The first and second processes can take place both in a hot hadronic gas created by the 
heavy ion collisions and in the expanding early Universe. In particular for the heavy 
ion reaction case, our analysis, coupled to the hydro dynamical expansion of the system 
will furnish additional information of the dynamics of the system. We will study baryon 
resonances evolution in heavy ions collisions in next chapter. The relaxation time for 7r° 
decay remains close (within 50%) to relaxation time in vacuum for large temperature range. 
In chapter [6] we will apply this for 7r° evolution in e+e~7 plasma, created by the intensive 
laser pulse and in early universe. 
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CHAPTER 5 

RESONANCE PRODUCTION IN HEAVY IONS COLLISIONS 
5.1 Introduction 

Hadron resonances are produced copiously in the quark-gluon plasma (QGP) fireball break 
up into hadrons (hadronization, chemical freeze-out) e.g. at RHIC [76-78, 80] . Within the 
statistical hadronization model (SHM) approach [57,58], the initial yields are described by 
chemical fugacities T, and hadronization temperature T. The production of heavy reso- 
nances is suppressed exponentially in m/T. Once formed, resonances decay. If this occurs 
inside matter, detailed balance requires also production of resonances, called 'regeneration' 
and / or 'back-reaction' . 

If the chemical freeze-out occurs much earlier than thermal, the initially produced 
resonances are practically invisible due to rescattering of decay products [81] . The observed 
yield of resonances is fixed by the physical conditions prevailing at the final breakup of the 
fireball, at which time last scattering occurs, this is the 'kinetic freeze-out'. The present 
work addresses two questions: 

a) how observable resonance yield depends on the difference between chemical freeze-out 
temperature (e.g. point of hadronization of QGP) and the kinetic freeze-out temperature; 

b) how this yield depends on the degree of initial chemical non-equilibrium at hadroniza- 
tion. 

One can see this work as an effort to improve on the concept of chemical freeze-out for 
the case of resonances: given the relatively fast reactions their yield remains sensitive to 
the conditions prevailing between chemical and thermal freeze-out, even if this time is just 
1 fm/c. 

Hadron resonances are observed in a surprisingly large yield when a quark-gluon 
plasma (QGP) fireball breaks up into hadrons [44,76-80]. This is unexpected, since the 
invariant mass signature formed from decay products could be erased by rescattering of the 
strongly interacting decay products [81]. In order to describe evolution of the resonance 
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abundance one can perform a microscopic transport simulation of the expanding system. 
In this approach the regeneration of resonances was previously studied by Bleicher and 
collaborators [45-47] . There are many detailed features of particle interactions to resolve in 
a microscopic model description and thus it seems appropriate to simplify the situation. We 
study resonance decay and regeneration using the momentum integrated population master 
equations, and assuming hydrodynamic expansion inspired model of fireball dynamics with 
conserved entropy content. In all our considerations we presume that the yield of pions tt 
is so large that we can assume it to be unaffected by any of the reactions we consider, thus 
we fix pion yield in terms of fugacity and temperature values. As a result, the final short 
hfespan resonance yield can be considerably different from statistical hadron gas (SHG) 
benchmark expectation. 

The other result, we obtain here, is that the long lived resonances, such as A(1520), 
can be considerably suppressed in their yield. This effect is amplified for the case when 
the initial hadron fugacities, and thus particle yields, are above chemical equilibrium. 
This situation is expected for a hadronizing QGP phase. The low A(1520) yield has been 
reported both in RHIC and SPS experiments [76, 79] . 

Here we consider two models. The first simplified model we apply for A(1232) /Ntot 
and E(1385)/Atot ratios calculation. In this model we consider only one dominant (fastest) 
reaction for each resonance: 



Then, for the 'fast' baryon resonances considered here we keep the sum of yields constant: 



The baryon annihilation, strangeness exchange such asA^ + ir<->A + 7r reactions, and 
population exchanges with higher resonances are assumed not to have a material impact 
within the time scale during which the temperature drops from chemical to kinetic freeze- 
out condition. 



A(1232) ^ AT + tt; 
S(1385) ^ A + TT. 



(5.155) 

(5.156) 



A + AT = Ao + A"o = A^o°* = Const., 
E(1385) + A = Eo(1385) + Aq = A[,°* = Const.. 



(5.157) 
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The experimentally observable hyperon yield appearing in our final result is 

Atot = S(1385) + A + E°(1193) + Y* (5.158) 

due to experimentally inseparable S°(1193) 7+A decay and the decay of further hyperon 
resonances Y* . Similarly, when we refer to A^tot we include baryon resonances in the count. 

In this model we also do not take into account medium effect on the reaction rates. 
It will be added in the second model. These effect is described in section 14.4.21 

The second model we apply for A(1520)/Atot and S(1385)/Atot ratios. In this 
model we include many reactions, which we will describe in section 15.3.11 In case of 
A(1520) (which is suppressed) it is necessary to consider a few reactions when for S(1385) 
the resut does not change much compared to first model. The resonance suppression, or 
enhancement, mechanism works as follows. In thermal hadronic gas the reaction (11.481) . 
can occur in both directions: the resonance decay 3 1 + 2, and the back-reaction 
(regeneration) resonance formation 1 + 2 — 3. When the reaction goes with the same rate 
in both directions, we have chemical detailed balance, e.g. particles yields do not change 
in this period of temporal evolution of the system. This does not necessarily mean that 
we have a chemical equilibrium. Instead it may be a transient condition for which none of 
the three particles is equilibrated chemically - we will show when this can happen. 

In the study of resonance decay and regeneration we are using the momentum 
integrated population master equations. We assume a fireball expansion model governed 
by hydrodynamic inspired flow with conserved entropy content. In our considerations we 
presume that the yield of pions vr is so large that we can assume it not to be materially 
affected by any of the reactions we consider. Thus we fix pion yield in terms of an ambient 
fugacity and temperature value, and in essence the total (per unit rapidity at RHIC) yield 
is fixed since we conserve entropy. 

An important assumption implied below is that the rapidly expanding hadron 
system maintains for the relevant particles a fully thermal (Boltzmann) momentum distri- 
bution. To describe the evolution of hadron abundances in the kinetic phase we track in 
time the yields of single strange hadrons after their initial formation. This is implemented 
in terms of time dependence of the chemical fugacities T(t), and the time dependence of 
the hadronization temperature T[t). 
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We look in detail at three potential evolution scenarios: 

a) a high temperature breakup at Tq ~ 180 MeV where the entropy content of the equili- 
brated QGP and HG-phase are similar; 

b) the To ~ 160 MeV case where chemical non-equilibrium among produced hadrons is 
already required; and 

c) at To ~ 140 MeV which is favored by descriptions of stable hadron production, and in 
which case a strong chemical non-equilibrium situation arises. 

For the late stage of the expansion, at relatively low density the assumption of 
thermal momentum distribution may not be anymore fully satisfied. In particular pions 
of high momentum could be escaping from the fireball. For this reason we will consider 
here a second scenario, which we call "dead channel". In this scenario we assume that 
the reaction fl4.124p goes mainly in the direction of resonance 3 decay and the resonance 
formation is switched off for 

ms - (mi + ma) > 300 MeV. (5.159) 

Without a complete kinetic model including equilibration and particle emission 
we do not know the exact energy in condition (15.1591) and timescale (during expansion) 
for which Boltzmann distribution is violated and dead channels appear. It is possible 
that reality lies between the two cases (kinetic Boltzmann distribution and dead-channels) 
considered here which, in our opinion, are the two most extreme limits. 

In section 15.21 we calculate A(1232)/A^tot and T,{1385) / Atot considering reac- 
tions (15.1551) and (I5.156p . However, in section 15.3.11 we investigate many further reac- 
tions in which resonances A(1520) an S(1385) participate. Thus we are obliged to develop 
a completely numerical evolution, for which the analytical study of S(1385)/A° provides 
a benchmark check of our approach. We discuss the temporal evolution of HQ particle 
fugacities T(t) in section 15.4.11 In section 15.4.21 we present results for the evolution of 
particle S(1385), A(1520) multiplicities during kinetic phase. In section [5.4.31 we obtain 
the observable 'ob' ratios A(1520)ob/Atot and S(1385)ob/Atot- We discuss our results in 
section 15.51 
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5.2 Short lived resonances A(1232) and S(1385) (simplified model) 



5.2.1 A multiplicity evolution equation 

In the following we will be referring explicitly to the A yield governed by cta = 1 /Fa = 1-67 
fm. All equations apply equally to S(1385) yield (partial decay width F^^a — 35 MeV) 
and we will compare our results with experiment for this case. We note that even though 
the S(1385) decay width is much smaller than Fa, the number of reaction channels and 
particle densities available lead to a significant effect for E(1385), comparable to our finding 
for A. 

The evolution in time of the A (or S(1385)) resonance yield is described by the 
process of resonance formation in scattering and decay, population equation (I4.114p . where 
particle 3 is A(1323) or S(1385), particle 1 is the ground state or A, and particle 2 is 
a pion. Allowing for Fermi-blocking and Bose enhancement in the final state, the two in- 
matter rates dW^-K^A/dydt and dWA^Nn/dVdt are described by Eq. fl4.118p and (14.1191) . 
The distribution functions for S, A, A^, A are Fermi and Bose for pions, Eq. (11.271) . 

Using detailed balance equation (14.1241) the master equation, Eq. (14. 1141) . can now 
be cast into the form: 

1 dN^ f T^Tat ^ \ rfPFA^iVTr -|^gQ^J 



V dt V Ta / dVdt 
This is a rather intuitive and simple result, yet only recently the 1^2 population master 

equations have been considered [3]. Equation (15.1601) implies for dN^/dt = the chemical 

equilibrium condition: 

This equation is solved by the global chemical equilibrium T^^ = = = 1. However, 
there are also other, transient, equilibrium states possible, given a prescribed value of e.g. 
the background pion abundance, T"^ ^ 1. When the initial state is formed away from 
transient equilibrium condition, we recognize that for Ta < T^T^r the A production is 
dominant, and conversely, for Ta > '^■k'^n the A decay dominates. 

We now introduce into the population master equation (I5.160p the effective lifes- 
pan, Ta aiming to find an equation similar to classic radioactive decay population equation. 
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We define the in medium A-lifespan to be: 

""^ = 'VdWA^j.^/dVdf ^^-^^^^ 

We recognize that in the Boltzmann hmit this corresponds to the ratio of equihbrium yield 

to the rate per unit time at which the equilibrium is approached. We obtain for Eq. (15. 1601) : 

dNA „ „ . dNA 1 . . 

— — = (T^T;v - Ta) -tt; . (5.163) 

dt dl Ta 

In case that the ambient temperature does not vary with time, and thus only 

populations evolve due to change in fugacities, we have dN/dt = dN/dT dT/dt and the 

following dynamical equation for the fugacity arises: 

TA^ = (T^T^ - Ta) . (5.164) 

This is 'classical' population equation form where the fugacity plays the role of the clas- 
sical densities. When the dynamical values of Tj(t) are used in the quantum Bose/Fermi 
distributions, the effects of blocking, and stimulated emission are explicit. 

If we instead were to introduce the lifespan by ta = {NAlV)/{dWA~*NTTldVdt), 
this implies for all particles (Bose, Fermi, Boltzmann) the classical population equation, 
e.g. dNA/dt = (TttTjv/Ta — l)iV'A/TA, and the quantum effects are now hidden in the 
definition of f. Both definitions coincide for the case of a dilute system, and differ most 
for dense systems. In the limit of very dilute, vacuum system, the relaxation time is the 
same as the lifespan of the particles. The computed yields of particles as function of time 
are not dependent on the finesse of the relaxation time definition. 

We now set up for semi-analytical solution of master equation ( 15.163^ . For mul- 
tiplicities A and considering the small yield and m ^ T we will use the Boltzmann 
distribution: 



Na _ ^ 



^ - TN-z-^gNX%K2{xN), (5.166) 



'^Ai-^9Ax\K2{xA), (5.165) 

Nn ^ T' 
T - ^^2^2 

where xa,n = ^a,n/T, K2{x) is Bessel function. Considering that fugacities, temperature 
and volume vary in time, we rewrite the left hand side of Eq. ( 15. 1631) : 

dNA _ dTA ^ dlnjxlK^jxA)) ^ ^ djVT') ^^^^ 



NAdr TAdT dT VT^ dr 
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We changed from t to r to make explicit the fact that we work in fluid-element co-moving 
frame and thus do not consider the effect of flow on the volume time dependence. 
Combining Eg. (15. 1631) with Eq. (l5.167|) we obtain 

^ = (T.T^ - Ta) — + Ta- + Ta-, (5.168) 
dr Ta Tt ts 

1 d\n{xA'^K2{xA)) 



Tt dT 
1 dln{VT^) 



T. (5.169) 
t. (5.170) 



Ts dT 

The last term is negligible, ^ ta, tt since pions dominate and we have near conservation 
of entropy which for massless particles would in fact imply VT^ =Const. 

Since entropy must be (slightly) increasing, while T is decreasing with time, ts > 0. 
Similarly, tt > 0, since the temperature decreases with time, and x'^K2{x),x = m/T 
increases with T : 

x^K2{x) ^ V0.57ra;^/^exp(-a;); (5.171) 



1 mA / 3 T \ T , . 

1 - o— ••• 7^- (5.172) 



Therefore: 

Tt T \ 2mA J T 
We now evaluate the magnitude of tt invoking a model of matter expansion of the 
type used e.g. in [26], where the longitudinal and transverse expansion is considered to be 
independent. In this model 

^ = nRliT)T, (5.173) 
dy 

where r is the proper time in the local volume element, this is exact for a 1-d ideal hydro 
flow. The growth of the transverse dimension can be generically described by 

R±{r) =Rq+ [ v{T')dT\ (5.174) 

J TO 

where we take velocity 

2 

v{t) = Vmax- arctan(4(r - To)/tc), (5.175) 

71 

where Vmax ~ 0.5 — 0.8c (c is speed of light and we take c = 1), relaxation time Tc ~ 0.5 
fm. To ~ 0.1 — 1 fm. 
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In the proper rest frame of the outflowing matter, 

^ oc = nRl{T)T^^ ^ Const.. (5.176) 

dy dy dy 

We will use dz/dy ~ r. 

The growth of the transverse dimension can be generically described by Eq. 

(imi). 

From Eq. (15.1 74p and Eq. (15. 1761) by elementary evaluation we obtain: 

T ^ _1 ^2WM±i) . (,177) 

Equation (15.1771) evaluated near hadronization condition is yielding the magnitude 
of see Eq. (l5.172|) . If the maximum expansion velocity is practically instantly achieved, 
vt/R_[_ ~ 1. This leads to maximum value of T/T ~ — 1/t. However if a more realistic 
profiles are assumed, T/T is diminished in magnitude as much as 30%. We thus conclude 
that 

0.5 rriA 1 0.7mA 
Th T Tt Th T 

which for hadronization time Th < 10 fm can compete with the width of the A-resonance, 
l/^A — 120 MeV. As this shows, the details of the expansion model are not critical for the 
results we obtain. In actual calculations we employ v{t) described in [26], where we assume 
that the expansion is already at maximum velocity at the time of chemical freeze-out. The 
resulting dependence T(r) after chemical freeze-out is shown in figure 15.321 We note that 
the time between chemical and thermal freeze-out At is not longer than about 2.5fm/c, 
and can be as short as Ifm/c. However, even such a short scattering period is enough to 
alter the visible yields of strong resonances, in fact most pronounced effect we find in the 
latter case, since the longer time allows a greater degree of chemical equilibration. 
We now can solve Eq. (l5.197p . Employing Eq. (l5.157p we have: 

^ + f(r)TA = g(r), (5.178) 



T(t) 



1 1 

ta Tt' 



(5.179) 



A^tot 1 
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Figure 5.32: Temperature T as function of 5t, the proper time interval between chemical and 
thermal freeze-out or chemical freeze-out temperature (from top to bottom) T = 180, 160, 140 
MeV and thermal freeze-out T > 120 MeV. 

where and are densities of A and N resonances with Ta = T^r = 1. The solution 
of Eg. (15.2011) is elementary: 



TA(r) 



■Y-O 



qe-'^h dr e -^^^ 



Th 



(5.181) 



where Th is initial expansion time at hadronization, and tu < t < Tmax, upper time hmit 
chosen to yield Tmax = 120 MeV, i.e. Tmax — 8 fm. 

5.2.2 Results for A(1232) and S(1386) resonance multiplicities 

In order to evaluate the final A multiplicity we need also to know initial particles den- 
sities right after hadronization which we consider for RHIC head-on Au-Au collisions at 
y^SNN = 200 GeV. We introduce the initial hadron yields inspired by a picture of a rapid 
hadronization of QGP with all hadrons produced with yields governed by entropy and 
strangeness content of QGP by quark recombination. In this model the yields of mesons 
and baryons are controlled by the constituent quark fugacity 7^: 



^0 



7ji 



i A,7V = 7a- 



(5.182) 
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Thus for 7g > 1 we always have the initial condition 

^ =7j>l (5.183) 

and the yield of A will increase in the time evolution. 

For each entropy content of the QGP fireball, the corresponding fixed background 
value of 7g can be found once hadronization temperature is known, see section 12. 3.41 
For T = 140 MeV pions form a nearly fully degenerate Bose gas with 7^ ~ 1.6. In 
the following discussion, aside of this initial condition, we also consider the value pairs 
T = 150 MeV, 7q = 1.42, T = 160 MeV, 7, = 1.27, T = 170 MeV, 7^ = 1.12 and T = 180 
MeV with 7, = 1. 

We assume in this section that m ^ T the density A is relatively low, thus there 
is no significant dependence of I/ta (the same for E(1385)) on T and Ta; in essence 
ta = h/T/^ takes the free space value ta — h/120MeV. Although from section |4:.4.2I we 
know that noticeabale but not very large large effect on rs(i385) from dense pion gas exits. 
We will take this effect into account in section , where more detailed model is presented. 

As already noted, we do not need to follow the evolution in time for the pion yield, 
which is fixed by conservation of entropy per unit rapidity, as incorporated in Eq. (15.1771) . 
Thus it is (approximately) a constant of motion. This can be seen recalling that the entropy 
per pion is nearly 4 within the domain of temperatures considered. Thus the conservation 
of entropy implies that pion number is conserved. With VT^ ~ Const., this further implies 
that during the expansion 

= 7J = Const., 

which we keep at the initial value. 

In figure 157331 we present results for ratios A/Aq (solid lines) and N/Nq (dashed 
lines) as functions of temperature T, beginning from the presumed initial hadronization 
temperature T through T^ax = 120 MeV. Aq and Nq are the initial yields obtained at each 
hadronization temperature. For T < 180 MeV, initially Ta < T^rT^, thus based on our 
prior discussion, we expect that the master equation leads to an initial increase in the yield 
of resonances. However, as temperature drops, due to the dynamics of the expansion the 
increasing yield of A turns over, and a final net increase of resonance yield is observed for 
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Figure 5.33: The ratio A/ (solid lines) and N/N^ (dashed lines) as functions of temperature 
T for select given pairs of values T, 7^, see text and figure box for details. 

T < 160 MeV. We note that for T > 180 MeV there is a continuous depletion of resonance 
yield. The nucleon yields move in opposite direction to the A-resonance. 

This behavior can be understood in qualitative manner as follows: The total num- 
ber A + is conserved therefore A multiplicity increases and multiplicity decreases 
until they reach transient chemical equilibrium {dN^^/dr = 0), corresponding to the maxi- 
mum point seen for A in figure 15.331 There is also influence of expansion: even if for some 
temperature the transient equilibrium condition fl5.16ip is reached, the system cannot stay 
in this equilibrium, Ta and T^v are increasing to conserve total number of particles. Ta 
increases faster because of larger A's mass. After Ta becomes larger than Tt^Tn A decay 
begins to dominate and their multiplicity is decreasing. The special case at hadronization 
temperature T = 180 MeV where, Tj = 1 and equilibrium condition is satisfied initially. 
As expansion sets in, A is decreasing because Ta > Tn (recall that here T^ = 1). In the 
SHM evaluation of yields one assumes that all ratios seen in figure 15.331 are unity. 

The initial hadronization yields which we used as reference in figure 15.331 are not 
accessible to measurement. Therefore, we consider in figure [5341 the fractional yield A/A'tot 
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(top frame), again as a function of temperature T. The results for hadronization tempera- 
tures To = 140 (solid blue line), To = 160 (dash-dot green line) and To = 180 MeV (dashed 
brown line) are shown. A'tot is fixed by hadronization condition and is not a function of 
time, as discussed. Thus the observable final rapidity nucleon yield corresponds to the ini- 
tial value at hadronization. Note that up to strange and multi strange baryon contribution, 
A'tot is the total baryon (rapidity) yield. 

Since in this study we have considered a subset of all relevant baryon resonances our 
chemical equilibrium reference yield (line for Ta = T^t) is not the same as the correspond- 
ing reference line for the full statistical hadronization model (SHM) evaluation, obtained 
using SHARE2, and presented as 2A+/p (upper frame) and 3S+/A'^ (lower frame). The 
SHARE2-SHM value A++/j9 ~ 0.2 at T ~ 160 MeV is consistent with the STAR d-Au 
results [44]. Also, comparing our with the SHARE2 result we note that SHARE2 yield 
is larger at chemical freeze-out. The magnitude of the difference in the yields at time of 
chemical freeze-out provides a measure of the magnitude of the corrections we can expect 
to arise in the full treatment at thermal freeze-out and/or systematic error for these yields. 

The nature of these effects is different for the two yield cases considered: the 
presence of heavier resonances which cascade by way of A leads to an increase of the thermal 
freeze-out yield. The correction is thus nearly as much as we see the SHARE2 yield higher 
at chemical freeze-out. For S(1385) the difference with SHARE2 arises from a difference of 
contributions of partial decays producing Atot, thus the correction is multiplicative factor 
which does not change, but is uncertain in magnitude due to lack of knowledge about the 
branching ratios. 

We believe that the A and E(1385) yields are underestimated by about 15% - 35%. 
(bigger effect for hadronization at higher T) . This implies that depending on hadronization 
temperature a relative yield range 0.16 < A+^/p = 0.5A/Atot < 0.26 arises, and similarly 
(see lower frame in figure 153^ 0.35 < S(1385)/Atot < 0.43 with the higher relative yield 
corresponding to the lower hadronization temperature. One of the key results of this 
work is the narrow range for S(1385)/Atot, and the fact that the initial chemical non- 
equilibrium effect leads to a reversal of the SHM model situation: the relative yields of 
massive resonances decreases with decreasing hadronization temperature. 
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Figure 5.34: Relative resonance yield, for (top) A/A'tot and (bottom) S(1385)/Atot as a functions 
of freeze-out temperature, for hadronization temperatures Tq = 140, 160, 180 MeV, see box and 
text for details. The dotted brown line gives the expected SHM chemical equilibrium result. 
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In order to compare with the experimental resuhs we note that the data pre- 
sented [76,77] are for charged S(1385), particle and antiparticle channels, (S=''(1385) + 
S='=(1385)/(Atot+Atot) — 0.29. This result needs to be multiplied with 3/2 to be comparable 
to results presented here which include S°(1385). Multiplying the value for hadronization 
at T = 140 MeV with thermal freeze-out at T = 120 MeV, and allowing for contribution by 
heavier resonances as indicated by SHARE2 our result is in perfect agreement with [76,77] 
However, given the narrow range of results we find, it seems that the high yield of S(1385), 
seen the error 0(20%) is nearly compatible with the entire range of chemical freeze-out 
temperatures here considered - the low T chemical freeze-out is favored by 1.5 s.d. over 
high T. 

The reader should take note that the 'thermal' model result presented in Ref. [76] 
corresponds to initial high temperature freeze-out in chemical equilibrium which is un- 
observable, since the high T hadronization resonance decay products have no chance to 
escape into free space. Thus this comparison of this model with experiment is flawed. The 
evolved yield is shown as (red) dashed line in figure [5311 and is found 25% below the value 
measured. The reason this happens is that the high T chemical freeze-out happens near 
chemical equilibrium and the yields follow closely the chemical equilibrium yield described 
by temperature, thus it is the thermal freeze out temperature which in this case controls 
the final observable resonance yield. 

5.3 Suppression of A(1520) and Enhancement of S(1385) 
5.3.1 Reactions scheme for A(1520) and S(1385) 

In figure [5735] we show the scheme of reactions which all have a noticeable effect on A(1520) 
yield after the chemical freeze-out kinetic phase. The format of this presentation is inspired 
by nuclear reactions schemes. On the vertical axis the energy scale is shown in MeV. There 
are three classes of particle states, which we denote from left to right as " N" (S=0 baryon), 
"S" {S = —1, 1 = 1 hyperon) and "A" {S = —1, / = hyperon). Near each particle bar we 
state (on-line in blue) its mass, and/or angular momentum and/or total width in MeV. The 
states A(1520) and S(1385) are shown along with the location in energy of A(1520) + tt 
and S(1385) + n respectively, both entries are connected by the curly bracket, and are 
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Figure 5.35: (color on line) Reactions scheme for A(1520) and S(1385) population evolutions. 
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highlighted (on-hne in red). The inclusion of the vr-mass is helping to see the kinetic 
threshold energy of a reaction. The lines connecting the A^, E,A columns are indicating 
the reactions we consider in the numerical computations. All reactions shown in figure [5T35] 
can go in both directions, as shown by the double arrows placed next to the numerical value 
of the partial decay width Fj in MeV. 

A(1520) decays with a total decay width of about 15.6 MeV, with two main chan- 
nels: 

S+ TT ^ A(1520), 6.5 MeV; (5.184) 

N + K ^A{1520), r^7 MeV. 

However, A(1520) reacts with several heavier S*-resonances, (S* = S(1670), S(1750), 
S(1775), S(1940), S(2030)): 

A(1520) + 7r ^ S*, (5.185) 

and these reactions have a larger reaction strength shown in figure [5351 A(1520) nearly 
behaves like a 'stable' hadronic particle since: 

a) it is dominantly coupled to heavier resonances; 

b) its natural lifespan is larger than the hadronic reaction rate. 

Hereto we note that (several) S* involved in Eq. (15.1851) participate in further 
reactions: 



A(1115) + 7r ^ S*; 


(5.186) 


S(1190) + 7r ^ S*; 


(5.187) 


N + K ^J:*; 


(5.188) 


S(1385) + 7r ^ S*; 


(5.189) 


A + K ^ S(1940,2030); 


(5.190) 


N + K{892) ^ S(1940); 


(5.191) 


E + T]^ S(1750). 


(5.192) 



All reactions shown above can excite S* resonances. Since the mass of A(1520) is near to 
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the S* mass, the yield of A(1520) is effectively depleted by the reaction chain 



A(1520) + ^ N + K, etc. 



(5.193) 



The balancing two step back-reaction can also occur, especially once A(1520) has been 
depopulated. Thus a dynamical reduced detailed balance yield of A(1520) would result if 
the system were at fixed volume rather than expanding. 

The multiplicity of S(1385) is mostly determined by its dominant decay and pro- 
duction in the reaction (I5.156P and to a lesser extent by the reaction 



The resonance S(1385) participates further in reactions with heavier S*; see reaction 
(15.1891) . but strength of these interactions is smaller than for similar reactions with A(1520) 
and smaller than the decay width of S(1385). Thus we find that the influence of these 
reactions on S(1385) yield is small. Another reason for a reduced effective depletion rate 
of S(1385) is that a lesser fraction of this resonance is needed to excite S*. Thus in such 
a reaction the depopulation effect decreases because of a larger mass difference between 
S(1385) and S* in comparison with A(1520) and S*. 

The reactions scheme for A(1520) reactions with dead channels is shown in fig- 
ure I5.36[ The difference between figure 15.351 and figure 15.361 is that some of the reaction 
lines have single-directional arrows, as is stipulated by the condition Eq. (I5.159p . 

5.3.2 Resonances densities, time evolution equations 

The evolution in time of the resonance yield is described by a master equation, similar 
to (14.1141) . where in general multichannel case the all processes of resonance formation in 
scattering is balanced by all natural resonance decay channels: 



where subscripts i, j denote different reactions channels when available. We further allow 
different subscripts i, j for the case where there are dead channels. Thus dWl^2^3/d'Vdt 



S(1190) + 7r ^ E(1385). 



(5.194) 




(5.195) 



3 



and dWl 



3^1+2 



/ dVdt are invariant rates (per unit volume and time) for particle 3 production 
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Figure 5.36: (color on line) Reactions scheme for A(1520) and S(1385) interactions in the "dead 
channel" model. 
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and decay respectively. In case all reactions occur in both directions the total number of 
fusion channels is the same as the total number of decay channels. 

Allowing for Fermi-blocking and Bose enhancement in the final state, where by 
designation particles 1 and 3 are fermions (heavy baryons) and particle 2 is a boson (often 
light pion) we have Eq.( I4.118P for resonance production and Eq. (14.1191) for resonanse 
decay rate. 

Using detailed balance Eq. fl4.124p we obtain for fugacity T3 the evolution equation: 



dr 



J2 r\T^L + f - + - - E 4 I ' (5.196) 




___5^TlT^,l+T3(l + l-$:i,), (5.197) 



Using detailed balance Eq. fl4.124p we obtain for fugacity T3 the evolution equa- 
tion [3, 7]: 

where characteristic time constants of temperature T and entropy S evolution, tt and ts, 
are from Eq. (15169|) and (ISTTOj) . 

The entropy term is negligible, ts ^ t-^, tt since we implement near conservation 
of entropy. We implement this in the way which would be exact for massless particles 
taking VT^ =Const.. Thus there is some entropy growth in HG evolution to consider, but 
it is not significant. In order to evaluate the magnitude of tt we use the relation between 
Bessel functions of order 1 and 2 (not to be mixed up with particles 1,2) d {z'^K2{z)) /dz = 
—z'^Ki{z). We obtain 

Tt a 2(3:3) T 

Tt > 0. We invoke a model of matter expansion. For a static system with ^ we see 
that Eq. (15.1970 has transient stable population points whenever 

J2r\Tl^-T,J2\ = ^- (5-199) 

'8 ■ Tn 

Next we address the functional dependence on time of Ti, T2. In the equation for 
Ti we have terms which compensate what is lost/gained in T3 see Eq. (I5.197p . Further we 
have to allow that particle '1' itself plays the role of particle 3 (for example this is clearly 
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the case for A(1520)). That allows a chain of populations relations as follows: 

(l' + 2'^ l) + 2^ 3, (5.200) 

Then we obtain: 

dTi ^ ^ 1 dNi/dTl ^ 1 dN^ldTl 

+ (5.201) 

The ratios of derivative of seen in the first line are due to the definition of relaxation 
time Eq. (??). The system of equations for baryons closes with the equation for Ty 

rfTi, ^ ^ 1 dN^/dTl ^ l_dNl/dTl 
dr ^ ^ r'^ dNv/dTv ^ ^' '^'r?dNv/dTv 

k n ' 

+ Tr (— + —]. (5.202) 

V^T TsJ 

In the present setting T2=7r =Const. by virtue of entropy conservation (see discussion 
below) and the same apphes to the case 2' = tt. However, if either particle 2 or 2' is a kaon, 
we need to follow the equation for T2,2'=k which is analogous to equation for particle 1 or 
1'. 

The evolution equations can be integrated once we determine the initial values of 
particle densities (fugacities) estabhshed at hadronization/chemical freeze-out. We deter- 
mine these for RHIC head-on Au-Au collisions at ^snn = 200 GeV. We introduce the 
initial hadron yields inspired by a picture of a rapid hadronization of QGP in which quarks 
combine into final state hadrons. For simplicity we assume here that the net baryon yield 
at central rapidity is negligible. Thus the baryon-chemical and strangeness potentials van- 
ish. The initial yields of mesons {qq, sq) and baryons {qqq, qqs) are controlled aside of the 
ambient temperature T by the constituent light quark fugacity 7^ and the strange quark 
fugacity 7^. 

The strangeness pair- yield in QGP is maintained in transition to HG. This fixes 
the initial value of 7^. In fact, since we investigate here relative chemical equilibrium reac- 
tions our results do not depend significantly on the exact initial value 7^ and/or strangeness 
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content. The entropy conservation at hadronization fixes 7^. For hadronization tempera- 
ture T{t = 0) = To = 180 MeV, 7^ = 1. However, when Tq < 180 MeV, 7g > 1 in order to 
have entropy conserved at chemical freeze-out. At Tq = 140 MeV 7^ = 1.6 that is close to 
maximum possible value of 7^, defined by Bose- Einstein condensation condition [1]. 

For reactions, such as shown in Eq. (11.481) . we have (lower index defines particle 
considered, where F = S, A is a hyperon) 

T?i=y) = 7J7., T°2=K) = 7,7s; (5.203) 

^('i=Ar) ^('2=710) defined by Eq. 15.1821 where the particle 1 in reaction (I1.48P is a 
baryon and particle 2 is a meson. The particle 3 is always a strange baryon: 

T?3=y) = 7?%, (5.204) 

As a consequence initially the pair of particles 1,2 reacts into 3. We have again condi- 
tion (15.183P satisfied. 

5.4 Numerical results 

5.4.1 Evolution of fugacities 

In order to evaluate the A(1520) and S(1385) multiplicities we must integrate Eq. (15.1971) . 
or Eq. (I5.20ip . or Eq. (15.2021) for each particle involved in figure 15.351 and perform simi- 
lar operations for reactions with dead channels in figure 15.361 This system of equations 
includes equations for A(1520), S(1385), five equations for E*s, equations for K(892) and 
A and equations for ground states A(1115), S(1190), N, K. All reactions in figures 15351 
are included. We solve this system of equations numerically, using classical fourth order 
Runge-Kutta method. 

Particle fugacities, except T^r, change rather rapidly. Figure 157371 shows the com- 
puted T(t) as a function of temperature T{t). We present here the scenario in which all 
reactions evolve in both directions, for the initial condition 7^ = 7^. The time, correspond- 
ing to the temperature shown at the bottom, is shown at the top of figure 15.371 in each 
frame. On the left we have hadronization at 140 MeV, in the middle at 160 and to the 
right at 180 MeV. Each frame has the same scale size for temperature unit, not time. For 
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Figure 5.37: (color on line) The fugacities T for selected particles are shown as a function of 
temperature T{t), for Tq = 140 MeV on the left, for Tq = 160 MeV in the middle and for Tq = 180 
MeV, on the right. See text for further details. 

Ts* we show two possible evolution examples, for S1750 (dash-dot dark line) and S(1775) 
(dashed line). These resonances have significant influence on the A(1520) yield. The solid 
hues are for Ta(i52o) (upper, red line) and Ts{i385) (lower, light blue line). The dash-dot 
and dashed light lines are for T2(ii9o) and Tao, respectively. The upper dotted line is for 
Tjy and lower dotted line is for T^. 

An important feature is that the Ts of massive hadron (resonances) increase very 
fast when T decreases. This is so since in absence of a rapid re-equilibration reactions, 
multiplicity of given resonance must be conserved. Then, according to Eq. (11.101) T.^ oc 
1/K2{mi/T), and thus for large rrii Tj oc exp(mj/T). We would expect Tj > T^, when 
rrii > rrij, and T decreases. This behavior is just like we found for the case of large charm 
fugacity [1]. However, because of the decay and regeneration reactions there are some 
deviations from this expectation in flgure 15.371 

For To = 180 MeV in most cases T3 > T1T2 {t > 0). Massive resonances decay 
to lower mass particles. The result is deflned by resonance mass, its decay width and 
decay products. For example T2(i775) is smaller than Ts{i75o) and Ta(i52o), because of its 
large decay width. Therefore excitation of S(1775) by A slightly dominates over S(1775) 
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T [MeV] T [MeV] 

Figure 5.38: The ratio S(1385)/S(1385)o on left and A(1520)/A(1520)o on rigiit as a functions 
of temperature T(t) for different initial hadronization temperatures To = 140, 160 and 180 MeV 
(blue/bottom, black/middle and red/top lines, respectively). Solid lines are for calculations with 
dead channels, dashed lines are for calculations without dead channels. 
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Figure 5.39: The rates for main channels of 11(1775) (on the left) and 11(1750) (on the right) 
decay and production as a functions of temperature T in the case when all reactions go in both 
directions and To = 140 MeV. Solid lines are for reaction S* ^ A(1520) + vr; dash-dot lines 
are for reaction S* N + K; dashed lines are for reaction S(1775) ^ A" + vr on the left and 
S(1750) <-> S + r/ on the right; blue and red lines are for decay and backward fusion reaction, 
respectively. 
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decay to A(1520) even in this case, when for most resonances the decay is dominant. For 
smaller initial hadronization temperatures Ta(i52o) becomes smaller than Te(i775), and even 
smaller than T2(i385) in some range of temperatures. This suppression occurs because of 
S(1775), and others S* regeneration. Because of large T^,, T2(i775) < TA(i52o)'^7r, the 
S(1775) production by A(1520) is dominant in the full range of T considered here. 

5.4.2 Final A(1520) and S(1385) multiplicities 

In this section we consider the evolution of the multiplicity of resonances A(1520), S(1385), 
S(1775) during the kinetic phase. We use the Boltzmann yield limit. By the symbol X{T) 
we refer to a particular resonance, and Xq is the initial multiplicity for that resonance. 
The dynamic yield of this resonance may be expressed as 

XjT) ^ Tx{t)T{tfK^{mx/T{t)) 
Xo TxoTlK,{mx/T,) ^^'^"^^ 

Figure 15.381 shows this yield as a function of T(t) for X = S(1385) (left) and 
X = A(1520) (right). We consider three initial conditions, temperature Tq = 140, 160, 180 
MeV, with corresponding 7^ = 1.6, 1.27, 1.0, respectively. The solid lines correspond for the 
model with dead channels and dashed one are for case when all reactions are symmetric in 
both directions. The thin dotted vertical line at T = 120 MeV marks the kinetic freeze-out 
temperature, assumed before in [4]. The main result is that the resulting relative yields for 
A(1520) and S(1385) behave qualitatively different from each other. In particular, as the 
temperature decreases, for the case Tq = 140 MeV we observe a strong yield suppression 
for A(1520), and a strong enhancement for S(1385) (as compared to initial SHM yields). 

To better understand the mechanism of A(1520) suppression, we analyze in some 
detail the case of S(1775) and S(1750) decay and production rates dW/dVdt. We assume 
here that these reactions can go in both directions. In figure 15.391 we show the reactions 
rates for the principal channels of decay and production as a functions of temperature T for 
S(1775) (left) and S(1750) (right), for the case of initial temperature Tq = 140 MeV which 
provides the largest A(1520) suppression. Solid lines are for the reaction S ^ A(1520) + 71, 
dash-dot lines are for reaction J] ^ N + K, dashed lines are for reaction S ^ A° + tt. Two 
set of lines are presented for the decay (on-line blue) and backward fusion reaction (on-line 
red), respectively. 
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As temperature decreases, all rates dW/ dtdV are increasing rapidly. This is mainly 
because fugacities T increase nearly exponentially when number of particles is conserved, 
see figure 15.371 We see that at the beginning of the kinetic phase all reactions go in the 
direction of S(1775) production, since S(1775) production rate is larger than its decay rate 
for all channels. Then at first S(1775) ^ A^+vr decay rate becomes dominant over S(1775) 
production rate in this channel, followed by the same for S(1775) ^ N + K channel. 

For the reaction S(1775) ^ A(1520) + 7r backward reaction is always dominant. As 
result, during the kinetic phase always more A(1520) resonances are excited into S(1775) 
than they are produced by S(1775) decay. The reason for this is the decay of S(1775) 
to the other channels, as long as Ts(i775) < TA(i52o)T'"7r- The lighter is the total mass of 
decay products, the earlier the decay reaction becomes dominant. This is due to the fact 
that the fugacity of T for heavier particles increases faster with expansion. Therefore, the 
decay rate becomes dominant earlier, when the difference between initial and final mass is 
larger. The net result is A(1520) suppression by S(1775) excitation. 

In figure 15.401 we show the yield of S(1775) normalized by its initial yield at 
hadronization: S(1775)/S(1775)o as a function of T{t). Like in the other figures above, 
solid lines are for the dead channels and dashed lines are for case when reactions go in both 
directions, solid (blue) lines are for To = 140 MeV, solid (black) lines for To = 160 MeV, 
and solid (red) lines are for Tq = 180 MeV. Each of the lines can be identified by their initial 
T-value. We see that when all reactions go in both direction the ratio S(1775)/S(1775)o 
increases at first similar to S(1385)/S(1385)o and A(1230)/A(1230)o ratios [4]. 

Compared to these ratios, S(1775)/S(1775)o ratio reaches its maximum value ear- 
lier, and after the maximum, the yield of S(1775) decreases faster. The reason for this 
behavior is that the mass of S(1775) is larger. The phase space occupancy Ts(i775), and 
therefore its decay rates, increase faster than the fugacity and decay rates for S(1385) and 
A(1230). Therefore decays S(1775) to some channels and its total decay rate become dom- 
inant earlier (see figure [5l39|) . Although the total decay width of S(1775) is approximately 
the same as for A(1230), the maximum value of this ratio is smaller. 

Said differently, the maximum yield of S(1775) does not have time to reach the 
value as high as that for A(1230). We thus learn that the time evolution of the yield of 
resonances with large decay width depends not only on their decay width, but also on mass 
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difference between initial and final states. Similar time evolution occurs for the other E*, 
which quantitatively depends on their mass, decay products masses and decay width. 

For most S*s, the decay products in the channel A(1520) + vr are heavier than 
the decay products in others channels, which are thus favored by phase space. For most 
resonances in our range of temperature, the decay into A(1520) + vr remains weak. The 
exception is S(1750) which decays also to S + r/, see figure 157391 (ms + m^ > mA(i52o) 
S(1750) begins to decay dominantly to A'^(1520) at relatively low temperature T = 116 
MeV, and continues to be produced by S + r/ fusion. 

As a result, allowing all reactions to go in both directions, the ratio 
A(1520)/A(1520)o has a minimum. This is specifically due to S(1750)) decay back to 
A(1520) at small temperatures as described above. However, when we satisfy Eg. (15. 1591) 
for dead channels the only decay occurs in the beginning of kinetic the dead-channel model 
phase. In that case the Ts*s are smaller, and the rate of reaction A(1520) + 7r S* exceeds 
the rate for backward reaction by larger amount, compared to the scenario without dead 
channels. This amplifies the effect of A(1520) suppression. In this case, S* decay to lighter 
hadrons right after they are produced by A(1520). We can see that for To = 140 MeV 
and To = 160 MeV A(1520) yield is always decreasing in the here considered temperature 
range. 

For S(1385) multiplicity we find a result quite different from A(1520) behavior 
discussed here, but similar to what we obtained in [4] by a very different method in a 
smaller basis set of states. In particular, the S(1385) yield is enhanced, but the maximum 
value of S(1385)/S(1385)o we find is a few percent higher, since we took into account 
the Bose enhancement of interaction rates, reaction (15.1941) . and S* production. Zl(1385) 
contribution to S* production is small, compared to the influence of the flrst two effects. 
The time (i.e. temperature) evolution of S(1385) practically does not depend on the 
presence of dead channels, and the maximum enhancement of S(1385) is even less sensitive. 
This in fact indirectly conflrms that S* has a small influence on E(1385) multiplicity. Thus 
we conflrm that: 

a) for To = 180 MeV S(1385) evolves with the system following the ambient temperature; 

b) for Tq = 160 MeV S(1385) shows some increase in yield; 

c) for To = 140 MeV there is a strong yield increase of S(1385). 
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Figure 5.40: The ratio S(1775)/S(1775)o as a functions of temperature T{t) for different initial 
hadronization temperatures To = 140, 160 and 180 MeV (blue/bottom, black/middle and red/top 
lines), respectively. Solid lines are for calculations with dead channels, dashed lines are for 
calculations without dead channels. 

While there is little sensitivity in the yield of S(1385) to issue of particle momen- 
tum distribution (little difference between the two models considered, dashed and sohd 
lines), the S(1385) yield is highly sensitive to initial hadronization condition. While for 
E(1385) the yield increases with decreased hadronization temperature, for A(1520) the op- 
posite is true, and in particular the smallest final A(1520) yield corresponds to the smallest 
hadronization temperature for both models. 

5.4.3 Experimentally measurable resonance ratios 

The initial hadronization yields, which we used as a reference in figure 15.381 in order to 
understand the physical behavior, are not measurable. What is commonly used as a 
reference for the yields of single strange hyperon resonances is the overall yield of the 
stable A°(1115), without the weak decay feed from S. Aside of the initially produced 
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Figure 5.41: The ratios i;(1385)/Atot (on left) and A(1520)/Atot (on right) as a function of 
temperature T of final kinetic freeze-out, for different initial hadronization temperatures Tq = 140, 
160 and 180 MeV (blue, black and red lines, respectively). Dashed lines are for calculations 
without dead channels, solid lines are for calculations with dead channels. The dotted purple line 
gives the expected SHM chemical equilibrium result. The dash-dot line is relative yield result 
result for SHM with Tq = T. 
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particles, the experimental yield of A°(1115) also includes resonances decaying during the 
free expansion after kinetic freeze-out, in particular (nearly) all decays of S(1385), and the 
experimentally inseparable yield of S°(1193) — 7 + A*^ decay and the decay of any further 
hyperon resonances Y*. 

Thus we normalize our final result with the experimentally observable final A^^^ 
hyperon yield: 

Atot = E°(1193) + 0.91S(1385) +A + Y*. (5.206) 

The factor 0.91 shows that 91% of end-state S(1385) decays to A. We also included in Atot 
calculations decays of H* A + K, which makes the result slightly dependent on 7s/7g 
ratio. We use 7s /7g = 1, since this ratio value is expected at top RHIC energy [1]. 

As noted, A(1520) and S(1385) experimentally observable yields also include any 
decays which occur in the free-streaming post-kinetic period. Thus we have: 

S(1385)ob = S(1385) + Y^dsss), (5-207) 
A(1520)ob = A(1520) + ¥^^,,20). (5-208) 

where ^^(lass) ^a(1520) hyperon multiplicities at kinetic freeze-out temperature, and 
which decay to S(1385) and A(1520), respectively. The multiplicities S(1385) and A(1520) 
are taken at the moment of kinetic freeze - out. 

In figure IS^in we present the fractional yields S(1385)/Atot (left), and A(1520)/Atot 
(right) as a function of temperature of final kinetic freeze-out T. The results for the 
hadronization temperatures Tg = 140 (blue lines), Tq = 160 (black lines) and Tq = 180 
MeV (red lines) are shown. Solid lines are for the case with dead channels and dashed lines 
are for the case when all reactions are going in both directions. 

In figure 15.411 the green dash-dotted line is the result when the kinetic freeze-out 
temperature T coincides with the hadronization temperature Tq. There is no kinetic phase 
in this case, only resonances decay after hadronization. This result is similar to SHARE 
result (purple, dotted line). The small difference is mainly due to us taking into account 
the decays 

S(1670,1750) ^ A(1520) + 7r, (5.209) 
which are expected/predicted in [82]. Similarly, for S(1385) our results for Tq = T are 
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different from SHARE results because we include the decay: 

E(1670) ^ E(1385) + vr, (5.210) 

expected/predicted in [83]. These additional resonances are part of current particle data 
set [84]. 

For all initial hadronization temperatures, as the freeze-out temperature decreases, 
the suppression for A(1520)ob/Atot ratio is larger than for A(1520)/A(1520)o (at the same 
temperature T of final kinetic freeze-out). This is particularly evident for dead channels 
and hadronization temperatures Tq = 160, 180 MeV (see figure 15381) . The effect is due to 
S(1775) suppression, as shown in figure [530] (and similar for other S*). For Tq = 140 MeV 
the additional suppression of A(1520), described above, is relatively small. 

For To = 140 MeV in the case without dead channels at final kinetic freeze-out 
T > 120 MeV, the final observed A(1520) suppression is even smaller, compared to its 
suppression in the kinetic phase at the same temperature (see figure [5738!) . The reason is 
that yield of S(1775) (and of the other S*s) is much enhanced for this range of temperatures 
see figure (5301 This additional S(1775) decays back to A(1520). That results in a smaller 
suppression at these temperatures. 

The above suppression effect increases in magnitude for higher hadronization tem- 
peratures, since the suppression of S(1775) and the sensitivity of A(1520)ob multiplicity 
to S* decays increase with temperature. However, when we consider dead channels (see 
figure 15.411) . the former effect of A(1520) suppression during evolution of kinetic phase 
increases for decreasing hadronization temperatures. Thus in the combined effect, the ob- 
servable relative suppression of A(1520)ob/Atot, is approximately of the same magnitude for 
all hadronization temperatures Tq. However, the initial hadronization yield of A(1520) is 
sensitive to temperature, and decreases rapidly with T. Therefore only for Tq = 140 MeV, 
a kinetic freeze-out temperatures ~ 95 — 105 MeV, and allowing for dead channels, the ra- 
tio Aob(1520)/Atot reaches the experimental domain Aob(1520)/Atot < 0.042 ±0.01 [76,79] 
shown in figure 15.411 by dashed lines. 

For the same initial conditions, that is for Tq = 140 MeV, we find the ratio 
E(1385)/Atot ~ 0.45 at T ^ 100 MeV (and for the entire range 95 - 135 MeV, in good 
agreement with experimental data [76,77]). In [4] this value of S(1385)/Atot is found at 
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T = 120 MeV, which was in the reference the presumed lowest possible temperature of 
the final kinetic freeze-out. Here we find that at T = 120 MeV the ratio S(1385)/Atot can 
be even higher (about 0.47), which is due to the Bose enhancement of in-medium E(1385) 
production rate (see discussion following figure [5l38|) . 

5.5 Summary and conclusions 

The resonant hadron states, considering their very large decay and reaction rates, can in- 
teract beyond the chemical and thermal freeze-out of stable particles. Thus the observed 
yield of resonances is fixed by the physical conditions prevailing at a later breakup of the 
fireball matter rather than the production of non-resonantly interacting hadrons. More- 
over, resonances, observed in terms of the invariant mass signature, are only visible when 
emerging from a more dilute hadron system given the ample potential for rescattering of 
decay products. The combination of experimental invariant mass method with a large 
resonant scattering makes the here presented population study of resonance kinetic freeze- 
out necessary. The evolution effects we find are greatly amplified at low hadronization 
temperatures where greatest degree of initial chemical equilibrium is present. 

Our study quantifies the expectation that in a dense hadron medium narrow res- 
onances are "quenched" [81] that is, effectively mixed with other states, and thus their 
observed population is reduced. Since we follow here the particle density, the effect we 
study is due to incoherent population mixing of A(1520), in particular with E*. This effect 
is possible for particle densities out of chemical non-equilibrium. However, this mixing can 
occur also at the amplitude (quantum coherent) level. As the result the yield suppres- 
sion effect could further increase, in some situations further improving the agreement with 
experiment. 

In first part of this chapter, we have presented master equation governing the evo- 
lution in time of the A, S(1385) baryon resonance yield after QGP hadronization, allowing 
for resonance decay and production process. We have shown considering the properties 
of the master equation that if the yield of hadrons is initially above chemical equilibrium, 
the resonance population increases beyond the initial yield. Conversely, we find that in a 
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physical system in which the particle multiplicities of hadrons arise below chemical equi- 
librium yields, a circumstance expected below threshold to QGP formation, the final yield 
of resonances is suppressed by the dominance of the resonance decay process over back 
reaction resonance production. 

In a quantitative model we evolved the yields after QGP hadronization allowing 
for initial chemical non-equilibrium particle abundances, and volume expansion assuring 
entropy conservation. We found, see figure 15.341 that the thermal freeze-out fractional 
resonance yield differs significantly from the chemical-freeze out SHM expectation, with 
the scenario involving high-T hadronization resonance yield being depleted, and low-T 
hadronization yield scenario further enhanced in relation to the total yield. 

The resonance enhancement effect we presented can only occur when the initial 
state is out of chemical equilibrium, and the decay/formation processes are fast enough 
to compete with the hadron volume evolution. One would thus think that 'narrow', i.e. 
quasi-stable resonances are not subject to the effects considered here. However, a special 
consideration must be given to narrow resonance which are strongly coupled to more mas- 
sive resonances which can decay fast into other channels. An example is A(1520), which 
is considered in next section. Aside of several specific predictions we made here, there are 
three important general consequences of our study: a) the fractional yield of resonances 
A* I A can be considerably higher than expected naively in SHM model of QGP hadroniza- 
tion, b) since there is nearly a factor of two difference in the final thermal freeze-out ratio 
in A/Atot, while the SHM yields a more T independent result, one can imagine the use of 
A/Atot as a tool to distinguish the different hadronization conditions e.g. chemical non- 
equilibrium vs chemical equilibrium a point noted in similar context before [?]; and c) we 
have shown that the relatively high yield of charged S^(1385) reported by STAR is well 
explained by our considerations with hadronization at T = 140 MeV being favored. 

In second model in this chapter our results show that the observable ratio 
A(1520)ob/Atot can be suppressed by two effects. First A(1520) yield is suppressed due to 
excitation of heavy S*s in the resonance scattering process. Moreover, the final A(1520)ob 
yield is suppressed, because S*s, which decay to A(1520), are suppressed at the end of 
the kinetic phase evolution by their (asymmetric) decays to lower mass hadrons, especially 
when dead channels are present (see figure 15^^01) . As a result, fewer of these hadrons can 
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decay to A(1520)ob during the following free expansion. A contrary mechanism operates for 
the resonances such as E(1385), A(1230). These resonances can be so strongly enhanced, 
that in essence most final states strange and non-strange baryons come from a resonance 
decay. 

We note that despite a scenario dependent resonance formation or suppression, 
the stable particle yields used in study of chemical freeze-out remain unchanged, since all 
resonances ultimately decay into the lowest "stable" hadron. Therefore after a description 
e.g. within a statistical hadronization model of the yields of stable hadrons, the under- 
standing of resonance yields is a second, and separate task which helps to establish the 
consistency of our physical understanding of the hadron production process. 

We conclude noting the key result of this study, that we can now understand 
the opposite behavior of A(1520) (suppression in high centrality reactions) and S*^1385) 
(enhancement, and similarly A(1230)) by considering their rescattering in matter. In order 
to explain both, the behavior of the A(1520)ob/Atot and E(1385)/Atot ratios, one has to 
consider T — 95 — 100 MeV as the favorite temperature of final kinetic freeze-out of hadron 
resonances, with Tq = 140 MeV being the favored chemical freeze-out (hadronization, QGP 
break-up) temperature. When there is little matter available to scatter, e.g. in peripheral 
collisions, the average value of A(1520)ob/Atot ratio is higher, approaching the expected 
chemical freeze-out hadronization yield for Tq — 140 MeV. All these findings are in good 
agreement with available experimental data. 
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CHAPTER 6 

RELATIVISTIC e+e-7 PLASMA CREATED BY LASER PULSE 

6.1 Freeze-out condition of relativistic e~, e"*", 7-plasma 
6.1.1 Introduction 

The formation of the relativistic , electron-positron-photon e~, e"*", 7 plasma (EP"^, temper- 
ature T in MeV range) in the laboratory using ultra-short pulse lasers is one of the current 
topics of interest and forthcoming experimental effort [48, 49] . 

For an expanding drop of plasma there is the freeze-out size R where the particle 
density p oc 1/R^ decrease allows the free-out-streaming of all particles, since the scattering 
length I oc 1/ (ap) grows with R^. Here we consider this freeze-out condition for a relativistic 
e~, e"*", 7 plasma. The conventional wisdom from keV temperature 'fusion' domain implies 
that an opaque plasma drop is not possible without lOO's of MJ of energy. 

Here we demonstrate a new temperature domain in which opaque plasma drops 
arc possible for the energy content of = 0.5 kJ with a radius in a range of i? = 2 10 nm, 
at a temperature at the scale of MeV. This new and interesting plasma domain arises since 
for T > rrie the density of electron-positron pairs grows rapidly and the scattering length 
I accordingly decreases rapidly. These physical conditions should become accessible in the 
foreseeable future upon the development of wavelength compression technology employing 
an optical wavelength laser beam reflected from a relativistic mirror, generated by a pulsed 
high intensity laser [51]. 

We evaluate mean free path length / of photon in EP"^ plasma for Compton scat- 
tering and pair production assuming thermal equilibrium. By comparing this length with 
plasma size at constant energy content we determine at what conditions plasma can be 
opaque and therefore in thermal equilibrium. For energy 0.5 kJ we study the limits for 
plasma size and temperatures. Similarly, we can also find chemical equilibration condi- 
tions, considering reactions of pair production and annihilation. We evaluate photon free 
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path using the method of thermal Lorentz invariant reaction rate, which was used before 
for strangeness production [16], [40]. This method allows us to take into account quantum 
effects in dense medium and easy to use in the observer rest frame. Note that corresponding 
equations can also be used in the astrophysical plasma environment. 

To create plasma drop with lager radius for given energy, the plasma temperature 
must be decreased. However we have a limit on the lowest temperature and therefore an 
upper limit on plasma size, defined by opaqeness condition at fixed plasma energy. For 
example, we will show at section 16.21 that fully plasma has maximum radius of 7 nm 
at temperature 2 MeV [and energy E=0.5kJ]. We found that to create low temperature 
(T < 0.5 MeV) opaque EP^ we need to deposit much more energy to a larger volume. This 
is due to the fact that the photon free path and, therefore plasma size, grow exponentially 
for this low temperature, when densities of photons with E > m and electron-positron pair 
are small. 

In susection 16.1.21 of this section we discuss statistical properties of EP'^, includ- 
ing master equation for electron - positron photon densities chemical equilibration under 
assumption that the particles are in thermal equilibrium. In subsection 16.21 we calculate 
the photon mean free path in plasma for Compton scattering and pair production and 
disscuss possible plasma size at given energy. In subsection 16.2.11 there are summary and 
conclusions 

6.1.2 Statistical properties of EP^ plasma 

Up to small QED interaction effects we can use Fermi and Bose momentum distribution 
f ll.27p . respectively to describe the particle content in the plasma 

•^^^^ ^ J^f,{u-p,±u,)/T _^ X' ■^"^ ^ Y^g«-p^/T _ I ' (6.211) 

When the electron chemical potential z/g is small, ^ T, the number of particles 
and antiparticles is the same, Ue- = ne+. Physically, it means that the number of e~^e~ 
pairs produced is dominating residual matter electron yield. Here we will set Ue = , and 
will consider elsewhere the case for very low density plasma where chemical potential may 
become important. Te(-y) is the fugacity of a given particle. 
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If plasma size is large enough, then plasma is opaque for photon electron scattering 
or pair production: 

7 + 7^e+ + e~; (6.212) 
7 + e^^7 + e^. (6.213) 

This plasma lives long enough and electrons and positrons are in thermal and 
relative chemical equilibrium with photons. The maximum photon density (black body 
radiation) is reached when photon fugacity = 1. Plasma is in chemical equilibrium, 
when Tg = 1, and all others particles fugacities are Tj = 1. Under this condition the 
plasma density has maximum value at given T. 

In our model we assume that the relative e+e" pair and photon yields are equi- 
librated by pair production and annihilation reactions. If we assume that thermal equi- 
librium establishes faster than chemical then the photon density evolution equations and 
chemical equilibrium conditions are similar to those for muon production considered in [7] : 

V^-(^e-n)R..^e^e-, (6.214) 



where 



_ 1 dW^^^e+e- _ 1 dWe+e-^^-y 

K^^^e+e- - ^2 ^Vdt ~ Tl dVdt ' 

7 ^ 



dW^^^e+e- / dVdt and dWe+e-^-yy/dVdt are Lorentz invariant rates for pair production and 
annihilation reactions, respectively. Then the EP"^ plasma is in relative chemical equilibrium 



for T-y < 1 when 



Te = = T. (6.215) 



At T ^ 1 we achieve full chemical equilibrium. 

We introduce pair production relaxation time defined by: 



^77 = — ' (6-216) 



then for the simplest case T{t) — const, T{t) — const and R{t) — const the equation for 



is 
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The rates and relaxation times are discussed in depth in [7]. 

If we introduce variable 7 = T^/T, which shows deviation from chemical equilib- 
rium, the equation for 7 is 

^ = (1-7')^. (6.218) 

The choice for the definition of relaxation time is made such that particle multiplicity 
reaches magnitude of equilibrium value during the time interval on the order of relaxation 
time. Note that the relaxation time for particle (7 in this example) production in two to 
two particles reaction increases by factor T^"^, where i is the initial particle in reaction 
(here e^). The physical reason why one introduces into the relaxation time is that 
the collision rate drops by that factor due to reduced density in plasma. 

In simple case, considered here, e^e~ pair production and annihilation are chemi- 
cally equilibrated when 

r^t « V. (6-219) 

where Tpi 0.1 ^ 10 fs is lifespan of plasma. 

The plasma properties, such as particle density and energy density can be evaluated 
using relativistic expressions: 

ni= I g^f^{p)d?p, E = [ Y,9^EiUp)d?pV, (6.220) 



where Ei = mf + p ^ , fi{p) is the momentum distribution of the particle 267, e^,/i^, 
7r°, TT^ and Qi its degeneracy: Qi = 1 for 7r° and Qi = 2 for the other particles, which can 
contribute. 

It is convenient to parameterize the equilibrium electron, positron and photon 
e~,e"'",7 plasma properties in terms of the properties of the Stephan-Boltzmann law for 
massless particles (photons). We present energy of plasma in terms of the effective degen- 
eracy g{T) comprising the count of all particles present at a given temperature T. Energy 
at T = 1 is 



E = giT)aT\ a = ^. (6.221) 

At temperatures T -C me we only have truly massless photons and g{T) ~ 2^. Once 
temperature increases beyond nie we find g ~ 2^+ (7/8)(2e- +2e+) = 5.5 degrees of freedom 
when T = 1. The factor 7/8 expresses the difference in the evaluation of Eq. fl6.220p for 



145 



the momentum distribution of Fermion and Boson Eq. (16.2111) . with Bosons providing the 
reference point at low T, where only massless photons are present. In principle these 
particles acquire additional in medium mass which reduces the degree of freedom count. 
However this effect is compensated by collective 'plasmon' modes. Thus we proceed with 
naive counting of nearly free EP^ components. 

In classical case, T << 1, we have for massless particles {m/T —>■ 0) 

E = 3NT, N = t4tTV, (6.222) 

where g = Q. There is no difference between Bose and Fermi particles. 

The densities and multiplicities of heavy particles {rrii » T) can be calculated 
using relativistic Boltzmann distribution: 

$ ^ = T,^Tm^K2{m,/T), (6.223) 

where subscript i G vr, /i, gi is the degeneracy, V is the volume K2 is the modified Bessel 
functions of integer order '2'. 

6.2 Mean free path of photon in e~^e — 7 plasma 

In order to be in thermal and chemical equilibrium plasma must be opaque for the reactions, 
which establish this equilibrium. The major reactions which may establish thermal and/or 
chemical equilibrium between photons and e~^e~ pairs are Compton scattering and pair 
production and annihilation. 

The mean free path of photon to produce e~^e~ pair is 

^77 = —r^ r = (6-224) 

where v is relative velocity of interacting particles cr^^^e+e- is cross section. For 1+2 —>■ 3+4 
reactions thermally averaged 

{vau^u) = , (6.225) 

nin2 

velocity v = c for photons scattering (we take c = 1). Mean free path length 1^^ is in 
the order of magnitude of r^!^. If size of plasma is on the order of magnitude of CTpi. 
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The condition of opaqueness for pair production is approximately the same as condition 
of chemical equilibration, Eq.( ]6.219|) . 

For Compton scattering mean free path is 

le, = / . = , (6.226) 

where Rg^ is Lorentz invariant Compton scattering rate. The plasma drop is opaque when 

« Rpi. (6.227) 

The equations for cross sections for pairs production and annihilation in center of mass 
frame are [86] 



a. 



77— >e"^e~ o 



{-A/x - 1)^/1 -A/x + 



m?x \ X 



■*77 



27ra2 (_4/a;_ i) 
m^x y/1 - 4/x 



m'^x{l — 4:/x) \ x'^ X J 1 — 1^1 — 4/a;' 

where a; = s/m^, s = {pi +^2)^- Note that there is extra 1/2 factor in equation for pairs 
annihilation because we have two identical particles or symmetrical wave function in final 
state. In rate R we add additional factor 1/2 when there are initial identical particles. In 
backward reaction there is also factor 1/2 from cross section. Therefore rate is symmetrical 
in both reaction directions. 

The cross section for Compton scattering in electron rest frame is 



27ra^ f l + uj 



e 7 9 I % 



2w(l + w) 
1 + 2u; 



-ln(l + 2cj)]) 



, 27ra^ / In 1 + 2^ ^+3^\ 

where uj = Eph/m, Eph is photon energy. 

In figure [632] we show cross sections fl6.228l) - fl6.230p as functions of total energy 
of particles in reaction Etot, Etot = s^^^ for pair production and annihilation and Etot = 
Eph + iTif, for Compton scattering. The decrease of cross sections with particles energy 
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Figure 6.42: The cross sections for pairs production and annihilation in center of mass frame and 
for Compton scattering in electron rest frame shown as functions of total energy of interacting 
particles Efot- 



increase can result that very high energy particle escape from plasma and do not participate 
in heavier particle production. We do not study this question in details here and we assume 
that this effect is small, leaving it for future research. 

The equation for rate R-y-y in Eq. (16.2241) is defined by Eq. fll.36p and (14.1181) in 
Introduction. Similar invariant rate calculations were done before for strangeness produc- 
tion in [16] and [40]. We perform integration similar to in [40]. Here we extend similar 
method to Compton scattering invariant rate Re-y calculation. 

The equation for Compton scattering rate is 



R. 



QeQ-y f d^pl f d^pj f d^pl f d^p. 



X 



'76 

E 

spin 



(27r) 



2E| 



2EI 



2EI 



i^5\pI + pI-pI-pI) 



u 



where s = (ps + pi) , t = (pi - p2 
pairs production s = (pi + P2Y, t = (ps — pi) 



Ps — P2)^, (t + u + s = 2rn?) (compared to 



u 



(Pa ~ P2)^, we cross s and t). ge = 4: 
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is degeneracy. The matrix element for Compton scattering is [38] 

2 2^2 



* m'^ — s m'^ — u 

oof 4m^ 4m^ — u — s\ , , 

levr^a^ ^ + ^ + ^ 6.232 

ym^ — s m'^ — u m'^ — s — u j 

where and (7e are photon and electron (positron) degeneracies respectively. We define: 

q = Pi + P3; p = ^(pi-P3); 

q=P4+P2; P' = ^(P4-P2); (6.233) 

then we have = — = q'q — q'"^ = s > and 

q / 1 

Pi = 2 + P5 P2 = -p + 2q; 

P3 = I - p; P4 = I + p'. (6.234) 

Using 

^ = I rfV(/ - m^)^(po) 
and p| 4 — = and 2 = O5 obtain: 

x/e(p°)exp(go/r). (6.235) 

The integrals from Eq. fl6.23ip can be evaluated in spherical coordinates. The angle 
coordinates are chosen with respect to the direction oi~q + 

q^l = (go,0,0,g), Pf, = {po,psme,0,pcos9), 
v'^jL = (Po,p'sin0sinx,p'sin0cosx,p'cos0). 
Using equations (16.234^ from delta functions we obtain equations: 

pI-P^ + - + Poqo - pq cos (9) = 0; (6.236) 

pI- + ^ - Poqo + pq cos (9) - = 0; (6.237) 

Po-p" + l+ Polo - P'q cos (0) - = 0; (6.238) 

Po-P'^ + I- Poqo + p'q cos (0) = 0; (6.239) 
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This system of equations is equivalent to (add and subtract pairs of equations fl6.236p and 
(K2m . (irasjl and (K2m ): 

2 

Po-p' + |-^ = 0; (6.240) 

9 

TTl 

Poqo -pq cos (9) + — = 0; (6.241) 

2 

Po'-p" + |-^ = 0; (6.242) 

2 

TTl 

Pogo-Acos(0)- — = 0; (6.243) 
then using S(f{x)) = J2i Vl/'(^j)l'^(^ ~ ^i)^ can rewrite integral (similar to [40]) as 

dq / dpo I 

'0 Jqi Jql 



X / dp / dpi / d{cos6) / (i(cos0) / 
Jo io i-i J-i Jo 



X 



2\ 1/2^ 



where 



-H-^-C)Mi+-)/^(i--)-=-^ 

xT~V. (f +Po) /e (f -Po) (lo/n (6.244) 



gi^2 and come from constrains cos^, cos0 < 1 and Eq. (16 .2401) - (16 .243p . 

The integration over p, p', cos 6*, cos0 can be done analytically considering the 
delta-functions. The other integrals can be evaluated numerically. In the order to simplify 
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numerical integration we introduce dimensionless variables: 



g = (g2 _ my'^z; (6.247) 

< z < 1; -1 < x{y) < 1. 

In these variables, using Eg. (16 .2401] - (16 ■243p . we obtain for u and t: 

s Tfl^ s ( Tn^ \ ^ 
u = {p-p'f = ^ 2^ ^ ^ 2\ ^ ~] ^^^^ \/(l -x2)(l -y2) sinx)] 

4 / 2 \ 2 

t = (p + p'f = ^2 _ I _ ^ _ _ !^ j _ ^(i_a;2)(i_^2)sin^^_250) 



Then limits for t are > t > — s + 2m? — m? / s. 

In new variables fl047D - (16:2491) the equation (lOTSj) is 

i?e^^e^ = j ^ dqo exp(go) ^ dzz'^{ql - vr?fl'^ dx j dx J dy (l - 

In figure [6]l3] the thermally averaged products (va) in observer frame, calculated 
using Lorentz invariant rates (Eq. (l6.225p ). are shown for Compton scattering (solid line), 
pair production (dot-dash line) and pair annihilation (dashed line). For pair production 
and Compton scattering v = c = 1. We see that at T << m for Compton scattering (vcr) 
goes to Thompson limit: 

a = = 6.710 mb. 

3m2 

For electrons production (va) starts to decrease at T < m and goes to with T — * 0. The 
density of photons with energy larger than threshold for pairs production drops in the tail 
of Boltzmann distribution with temperature decrease. For electrons annihilation (va) stays 
finite at small temperature, f — > 0, because a diverges as 1/v. In figure [631] the mean free 
photon paths multiplied by T are shown for e~^e~ pair production reaction Tl^^ (Eq. (l6.224p ) 
and Compton scattering Tl^j (Eq. (16.2261) ) are shown as functions of temperature T for 
T = 1 (thick dashed and solid lines, respectively) and T = 0.1 (thin dashed and solid lines. 
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Figure 6.43: The thermal products {va) for Compton scattering (solid line), pair production 
(dash-dot line) and annihilation (dashed line) in observer rest frame shown as a function of 
temperature T at T = 1 




Figure 6.44: IT for Compton scattering and pairs production at T = 1 (thick dashed and 
solid lines) and T = 0.1 (thin dashed and solid lines) as functions of temperature T; radius of 
equilibrium (T = 1) plasma at energy 0.5 kJ (dot-dash). 
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respectively). Cases at T = 1 and T = 0.1 are shown to demonstrate the magnitude of 
the difference between these two cases due to quantum effects. For electron-positron pair 
production l^^ is suppressed slightly at T = 1 due to Bose enhancement of reaction rate. 
Both effects Bose enhancement and Fermi blocking contribute to the Compton mean free 
path . At T << 1, when there is no quantum effects from plasma, the free mean path of 
photon oc We see that l^^ < at T > 8 MeV for T = 1 and at T > 10 MeV for 

T << 1 and therefore the thermal equilibrium in EP'^ is established by reaction fl6.212p 
at T > 8 MeV, approximately at the same time with chemical equilibrium of pairs and 
photons, l^^ drops fast when temperature increases: 

^77 « (6-252) 

At temperature range 1 MeV< T < 10 MeV the temperature dependence of le-y is a little 
slower than 1/T^. In figure [634] we also show the radius of EP^ plasma drop with energy 
0.5 kJ and T = 1 as a function of temperature. We see that at these conditions plasma 
loose opaqueness at temperatures smaller than 2 MeV (thin vertical line). Photon free path 
increases faster with temperature decrease than plasma radius. Corresponding maximum 
radius of equilibrium plasma is about 7 nm {R/le-y ~ 3) at T = 1. 

The pair production or annihilation relaxation time Eq. fl6.216p is approximately 
2 10~^ fs at T = 2 MeV, in agreement with condition (16.2191) . Plasma drop with energy 0.5 
kJ is thermally and chemically equilibrated at T > 2 MeV and largest density with T = 1. 

From figure 16.441 we also see that the mean free path length of photon starts 
to increase exponentially at small temperatures T < 0.5 MeV. Therefore plasma size 
and energy at this temperature also increase exponentially. On the contrary the higher 
plasma temperature is the smaller photon free path length becomes. Higher temperature 
(T > 2 MeV) opaque plasma can have smaller size {R < 7 nm). When / (almost) satisfies 
Eq. (l6.252p . minimum radius of opaque plasma, allowed by opaqueness condition (I6.227P 
is also oc 1/T^. If we focus laser pulse energy in this small volume, we obtain that the 
necessary energy is 

E oc T'^R^ oc l/T^. (6.253) 
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6.2.1 Summary and Conclusion 

In this part we investigated physical conditions suitable to create opaque and, therefore, 
thermally and chemically equilibrated e"*", e~, 7-plasma drop. In order to address this ques- 
tion we evaluated Lorentz invariant rates for the Compton scattering and pair production 
in thermally and chemically equilibrated EP^ plasma. We then used these Lorentz invariant 
rates to evaluate the corresponding mean free path length / of particles. 

Comparing I to plasma drop size we showed that an opaque equilibrium density 
plasma drop can be produced at energy 0.5 kJ in the volume with largest possible radius 
R = 7 nm. This volume corresponds to the smallest possible temperature T = 2 MeV. 
In order to reach higher than 2 MeV temperature, we need to increase energy of plasma 
(which is proportional to T^) or/and decrease plasma size. At higher temperature opaque 
plasma can be created at the total plasma energy smaller than 0.5 MeV, since smaller 
plasma drop size is in agreement with the opaqueness condition Eq. fl6.227p . as seen in 
figure I6.44[ 

On the other hand in order to create opaque plasma with temperature lower than 
2 MeV, the necessary amount of energy is larger than 0.5 kJ. This is so because the plasma 
size has to be large in order to satisfy opaqueness condition Eq. (16.2271) . 

Our main result, perhaps unexpected at the first sight is illustrated in figure [6^441 
For the temperature T > 2 MeV equilibrium plasma production with a relatively small en- 
ergy pulse (compared to lower temperature equilibrium plasma) may be possible. However, 
the challenge here is to focus the energy into the volume of size < 10 nm. 

6.3 Pion and muon production in relatvistic e"'"e~7 plasma 
6.3.1 Introduction 

The elementary properties of EP'^ have recently been reported, see [85], where typical 
properties are explicitly presented for T = 10 MeV. One of the challenges facing a study 
of EP^ will be the understanding of the fundamental mechanisms leading to its formation. 
We propose here as a probe the production of heavy particles with mass m ^ T. Clearly, 
these processes occur during the history of the event at the highest available temperature, 
and thus information about the early stages of the plasma, and even pre-equilibrium state 
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should become accessible in this way. 

We focus our attention on the strongly interacting pions tt'^jTT^ {uIt^c^ ~ 140 
MeV), and muons iM^{m^c^ ~ 106 MeV), (m the following we use units in which k = c = 
h = 1 and thus we omit these symbols from all equations. Both, the particle mass, and 
plasma temperature, is thus given in the energy unit MeV.) These very heavy, compared 
to the electron (mgC^ = 0.511 MeV), particles are as noted natural 'deep' diagnostic tools 
of the EP^ drop. Of special interest is the neutral pion 7r° which is, among all other 
heavy particles, most copiously produced for T <^ m. The 7r° yield and spectrum will 
be therefore of great interest in the study of the EP^ properties. Conversely, the study 
of the in-medium pion mass splitting Am — mT^± — ttIt^o — 4.594 MeV at a temperature 
T > Am will contribute to the better understanding of this relatively large mass splitting 
between 7r°and n^, Am/m = 3.34%, believed to originate in the isospin symmetry breaking 
electromagnetic radiative corrections. 

However, given its very short natural lifespan: 

7r° ^ 7 + 7, r°o = (8.4 ± 0.6) lO'^^s. 

7r° is also the particle most difficult to experimentally study among those we consider: its 
decay products reach the detection system nearly at the same time as the electromagnetic 
energy pulse of the decaying plasma fireball, which is likely to 'blind' the detectors. 

This plasma drop we consider is a thousand times hotter than the center of the 
sun. This implies presence of the corresponding high particle density n, energy density 
e and pressure P. These quantities in the plasma can be evaluated using the relativistic 
expressions: 

= y 9rfi{p)d?P. (6.254) 
e= I Y.SiEifMd'p. E, = ^m\+p (6.255) 

P = \jY.9r{E,-'^^f,{p)dp\ (6.256) 

where subscript i e 7, e~,e+, 7r°, tt"*", 7r~, /i",/!'^, fi{p) is the momentum distribution of 
the particle i and and Qi its degeneracy, for i — e~, e"*", 7, we have gi — 2, and gi — 1 
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for 7r°, TT^TT^. For a QED plasma which hves long enough so that electrons, positrons are 
in thermal and chemical equilibrium with photons, ignoring small QED interaction effects, 
we use Fermi and Bose momentum distribution, respectively, Eq.f l6.2lTi) . 

It is convenient to parametrize the electron, positron and photon e~e~^'~f plasma 
properties in terms of the properties of the Stephan-Boltzmann law for massless particles 
(photons), presenting the physical properties in terms of the effective degeneracy g{T) 
comprising the count of all particles present at a given temperature T: 



we only have in this case truly massless For temperatures T ^ mg photons and g{T) ~ 
g'{T) ^ 2^. Once temperature approaches and increases beyond rrie we find g ~ g'{T) ~ 
2^+(7/8)(2e-+2e+) = 5.5 degrees of freedom. In principle these particles acquire additional 
in medium mass which reduces the degree of freedom count, but this effect is compensated 
by collective 'plasmon' modes, thus we proceed with naive counting of nearly free EP^ 
components. The factor 7/8 expresses the difference in the evaluation of Eq. (16. 2 56 p for 
the momentum distribution of Fermions and Bosons Eq. (16.21 ip . Bosons provides the 
reference point at low T, where only massless photons are present. In passing, we note 
that in the early Universe, there would be further present the neutrino degrees of freedom, 
not considered here for the laboratory experiments, considering their weak coupling to 
matter. 

In figure 16.451 we present both g{T) and g'{T), as a function of temperature T 
in form of the energy density e normalized by crT^, and, respectively, the pressure P, 
normalized by crT^/3 . The g{T) jumps more rapidly compared to g'{T), between the 
limiting case of a black body photon gas at T < 0.5 MeV {g = 2) and the case g = 5.5 
for 7, e~,e~^, since the energy density also contains the rest mass energy content of all 
particles present. The rise of the ratio at T > 15 MeV indicates the contribution of the 
excitation of muons and pions in equilibrated plasma. We note that the plasma produced 
pions (and muons) are in general not in chemical equilibrium. The distribution functions 
which maximize entropy content at given particle number and energy content are [88]: 



V 



e = g{T)aT' 



3P = g'{T)aT\ 




(6.257) 



1 



1 



(6.258) 



n 
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Figure 6.45: On left: the ratios g = e/aT"^ and g' = 2>P/aT'^ as a function of temperature T; 
on right: the equihbrium densities of electrons (blue, solid line), photons (green, dash-dot line), 
muons (red, dashed line), pions (blue dotted line) as functions of temperature T. 

where Tt^oi^t^±^j and are particles fugacities. The chemical equilibrium corresponds to 
T^o(^±) = = 1 used in figure [67i5] on right, since this is the maximum density that can 
be reached in the buildup of these particles, for a given temperature. We occasionally refer 
to /tt f- as the boson distribution function and to fj^ as the Fermi distribution 

function. For Tj ^ the quantum distributions shown in Eq. (16.2581) turn into the classical 
Boltzmann distributions, with abundance prefactor Tj. 

In the case of interest here, when T < m, we also consider consider the Boltzmann 
limit of the quantum distributions Eq. (16.2581) . that is to drop the 'one' in the denominator, 
Eq. (ll.7p and taking the non-relativistic limit Eq. (I1.12p we have: 

^^n^ = T^^TmlK^imjT) ^ f e"'"-/^ + . . . , (6.259) 

V 27r^ \ 2tx J 

The particle densities are shown on right in figure 16.451 The top solid line is the 
sum of 71^.+ + rig- , which is marginally bigger than the photon density (dashed, blue) which 



157 



follows below. We also include in the figure the sum density of muons + n^- (red, 
dashed) , and the density of the neutral pion 7r° (bottom solid line) , both of which appear 
comparatively very small in the temperature range of interest. However, in magnitude 
they rival the normal atomic density (~ 10^/nm^) already at T = 4 MeV, and 5 MeV, 
respectively. This high particle density in the chemically equilibrated plasma explains the 
relatively large collision and reaction rates we obtain in this work. In turn, this opens the 
question how such dense, chemically equihbrated EP^ state can be formed - we observe 
that colliding two ultra intense circularly polarized and focused laser beams on a heavy thin 
metal foil(s) is the current line of approach. Initial simulations were performed [87] . Many 
strategies can be envisaged aiming to deposit the laser pulse energy in the smallest possible 
spatial and temporal volume and this interesting and challenging topic will without doubt 
keep us and others busy in years to come. 

As it turns out, a small drop of EP^ plasma with a size scale of Inm is, given the 
high particle density, opaque. The mean free paths k of particles 'i' are relatively short. 
Where the reference energy values (31.1 and 27.5 MeV) correspond to the mean particle 
energy at T = 10 MeV. Photons are subject to Compton scattering, and electrons and 
positrons to charged particle scattering. In fact these values of li are likely to be upper 
limits, since Bremsstrahlung type processes are believed to further increase opaqueness of 
the plasma [90]. In our considerations plasma particles of energy above 70 MeV are of 
interest, since these are responsible for the production of heavy particles. We see that the 
mean free path of such particles has also nm scale magnitude. 

We note that a EP^ drop of radius 2nm at T = 10 MeV contains 13 kJ energy. 
This is the expected energy content of a light pulse at ELI (European Light Infrastructure, 
in development) with a pulse length of about At — 10~^^s. For comparison, the maximum 
energy available in particle accelerators for at least 20, if not more, years will be in head 
on Pb-Pb central collisions at LHC (Large Hadron CoUider) at CERN, in its LHC-ion 
collider mode, where per nucleon energy of about 3 TeV is reached. Thus the total energy 
available is 200 /iJ, of which about 10%-20% becomes thermalized. Thus ELI will have 
already an overall energy advantage of 10^, while in the LHl-ion case the great advantage 
are a) the natural localization of the energy at the length scale of 10~^ nm, given that the 
energy is contained in colliding nuclei, and b) the high repetition rate of collisions. 
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As a purely academic exercise, we note that should one find a way to 'focus' the 
energy in ELI to nuclear dimensions, and scaling the energy density with up from 
what is expected to be seen at CERN-LHC-ion (T < IGeV), we exceed T = 150 GeV, the 
presumed electro- weak phase boundary. Such consideration lead the authors of Refs. [48,49] 
to suggest that the electro-weak transition may be achieved at some future time using 
ultra-short laser pulses. 

Returning to present day physics, we are assuming here that T near and in MeV 
range is achievable in foreseeable future, and that much higher values are obtainable in 
presence of pulses with At < 10~^^s, cAt < 0.3nm. Hence we consider production processes 
for 77^^,77^,//^ for T < 50 MeV. We study here all two body reactions in EP^ which 
lead to formation of the particles of interest, excluding solely 67 —>■ en^, and the related 
e~e~^ 77r'^. The presence of a significant (1.2%) fraction of 7r° e+e~7 decays implies 
that these related two body processes could be important in our considerations. However, 
these reactions involve the ttq off-mass shell couping to two photons, which needs to be 
better understood before we can consider these reactions in our context. 

We also do not consider here the inverse three body reactions e+e~7 7r°, since 
there is no exponential gain in using n > 2 particles to overcome an energy threshold, here 
mTrO. The independent probability of finding n particles with energy rriT^o/n each is the 
same for any value of n: 

PiP2....Pn oc (e-"*-°/"^)" = e-^. (6.260) 

This resolves the argument that more particles could overcome more easily the reaction 
barrier, n-body reactions with n > 2 are in fact suppressed in EP^ by the weakness 
of the electromagnetic (EM) interaction, since adding an EM-interacting particle to the 
reactions process requires an EM- vertex with a = 1/137. Thus microscopic reactions 
in EP^ involving n > 2 are suppressed by a factor 100 for each additional EM particle 
involved in the reaction. This does not mean that a collective/coherent process of heavy 
particle production by many particles is similarly suppressed: for example fast time varying 
electromagnetic fields provide through E ■ B a collective source of 7r°. We defer further 
study of this production mechanism which requires multi MeV~^ range oscillation to be 
present in EP^. 
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In the following section, we introduce the master equation governing the produc- 
tion of pions and muons in plasma and formulate the invariant rates in terms of know 
physical reactions. In section 16.51 we obtain the numerical results for particles production 
rates and reactions relaxation times which we present as figures. In section 16.61 we discuss 
these results further and consider their implications. 

6.4 Particles production 
6.4.1 production 

TT^ in the QED plasma is produced predominantly in the thermal two photon fusion Eq ll.69[ 
see chapter 14.4.31 Much less probable is the production of ttq in the reaction: 

e- + e+^7r°. (6.261) 

These formation processes are the inverse of the decay process of ttq. The smallness of 
the electro-formation of ttq is characterized by the small branching ratio in ttq decay B = 
Fee/r^-y = 6.2 ± 0.510~®. Other decay processes involve more than two particles, vr" can 
also be formed by charged pions in charge exchange reactions. However, in EP'^ in the 
domain of T of interest we find that at first the neutral pions will be produced. These in 
turn produce charged pions. Therefore we introduce the pion charge exchange process in 
the context of charged pion formation in the subsection 16.4.31 

Omitting all sub-dominant processes, the resulting master equation for pion num- 
ber evolution is: 



V dt dVdt dVdt ' 

where N^^o is total number of 7r°, V is volume of the sysiem^d^W^^^T^o / dVdt is the 
(Lorentz) invariant 7r° production rate per unit time and volume in photon fusion, and 
d'^WTjO_,^^/ dVdt is the invariant 7r° decay rate per unit volume and time. The rates for tx^ 
decay and production can be calculated using Eq. ( 14.1191) and (I4.118p . 

We assume that in the laboratory frame the momentum distribution of produced 7r° 
are characterized by the ambient temperature. Eq. (16.2591) defines the relation of fugacity 
to the yield. This equation allows now to study the production dynamics as if we were 
dealing with a vr*^ in a thermal bath, and to exploit the detailed balance between decay 
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and production process in order to estimate the rate of 7r° production. This theoretical 
consideration should not be understood as assumption of equilibration of 7r°, which could 
upon production escape from the small plasma drop. 

Using the detailed balance relation and R^^^.„o definition (11.361) . Eg. (16.2621) can 
be written in the form: 

(T^ -T,o)i?^^_,o, (6.263) 



V dt 

For T,ro ^ = 1 we reach chemical equilibrium, the time variation of density due to 
production and decay vanishes. 

We introduce the pion equilibration (relaxation) time constant by: 

T^o = — . (6.264) 

Note that when the volume does not change in time on scale of r^o (absence of expansion 
dilution) and thus T is constant, the left hand side of Eq. (l6.263p becomes (in^o/rft. Given 
the relaxation time definition Eq. (16.2641) the time evolution for of the pion fugacity for a 
system at fixed time independent temperature satisfies: 

r.o^ = T^-T.o, (6.265) 

which has for T^o(t = 0) = the analytical solution T^o = (l — 6^*/'^'^°), justifying the 
proposed definition of the relaxation constant. 

The relaxation time Tpjo is calculted in section 16.2641 and shown in figure 14.311 
TpjO ^ r°o for temperatures considered, because the relativistic time dilution effect cancels 
with medium effect. 

The 7r° production rate is thus related to the decay rate 1 / r^o by the simple formula 



where in the last expression we have used Eq. (I6.259P in the limit m » T. It is important 
for the reader to remember that derivation of Eq (I6.266P is based on detailed balance in 
thermally equilibrated plasma, and does not require chemical equilibrium to be established. 

This exact result (blue, solid line) is compared to the approximate result Eq. (16.2661) 
(green, dashed line) in figure 16.461 We note that it is hard to discern a difference on 
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Figure 6.46: The ttq production rate (blue, solid line) and approximate rate from Eq. ()6.266p 
(green dashed line) as functions of temperature T. 

logarithmic scale, especially so for small temperatures where the only (small) effect is the 
relativistic time dilation. This implies that it is appropriate to use the simple intuitive 
result Eg. (16.2661) in the study of 7r° production. 

Before closing this section we note that we can use exactly the same method to 
extract from the partial width of the ttq e^e~ the reaction rate for the inverse process, 
which will be discussed below. All arguments carry through in identical and exact fashion 
replacing where appropriate the Bose by Fermi distributions and using Eq jl.35[ 

6.4.2 Muon production 

In the plasma under consideration, muons can be directly produced in photons or e^e~ 
fusions, reactions (ll.74p and (11.751) . For these reactions the master evolution equation 
developed for the study of thermal strangeness in heavy ion collisions applies [16,17,40,89] 
(compared to these references our definition is changed -Ri2-»34 — > l/('^i^2)-Ri2-^34 ) 




(6.267) 
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Like before for 7r° we consider the master equation in order to find appropriate definition 
of the relaxation time constant for /i^ production. In no way should this be understood 
to imply that muons are retained in the small plasma drop. The n production relaxation 
time is defined by: 

1 dn^/dT^ /fiofifi^ 

= --TB —~R V (6.268) 

where a suitable choice is a = 1,2 for t = 0, oo, respectively (see below). The form 
of Eq. (16.2681) assures that, omitting the volume expansion, i.e. the dilution effect, in 
chemically equilibrated EP^ the evolution of the muon fugacity obeys the equation 

= l-rl, (6.269) 

which has T^(t = 0) = the simple analytical solution [17]: 

= tanht/ar^. (6.270) 

For t oo, near to chemical equilibrium, — 1 — e^^*/'^^^ , while for t — ^ 0, at the onset 
of particle production with small we have = t/(ar^) . Hence, near to chemical 
equilibrium it is appropriate to use a = 2 in definition of relaxation time Eq. (16.2681) , while 
at the onset of particle production, more applicable to this work a more physical choice 
would be a = 1. However, following the convention, in the results presented below the 
value a = 2 is used. 

For invariant muon production rates we use, Eq J1.33[ with photons (bosons) or 
e~^e~ (fermions) in initial state. It is interesting to note that despite including of quan- 
tum effects (Bose stimulated emission and Fermi blocking), using rates as defined, we 
don't change the master population equation form. Only modification is slight fugacity 
dependence of rates presented in Eq. (14. 1181) . 

The ^ |Me+e-^^+;i- P differs from often considered heavy quark production 
^ |Mgg^ccP [91,92] {nic >> rriq) by color factor 2/9, and the coupling — > a of QCD 
has to be changed to a = 1/137 of QED. Then we obtain, based on above references: 



2 



tf + (m^ - uf + 2m^s 



\Me+e--,^i+ii-\ = g^SiT a , (6.271) 

s 

where m = 106 MeV is the muon mass, electron and positron degeneracy Qe = 2, and s, 
t, u are the usual Mandelstam variables: s = {pi + t = {P3 " Pi)^, = (Pa ~ ^2)^, 
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s + t + u = 2m^. For the total averaged over initial states |Mp for photon fusion we have 
\M,,_..^.-\' = gl%^^a^ -4 + +4 -— - + -— — + -— — + 



(6.272) 

where degeneracy = 2. Near threshold s ~ 4m^, with t,u ^ —m? we find 

|M^^^^+^- 12 = UTi^a\ |Me+e-^^v- r = 327r2«l (6.273) 

The e+e^ ^ fx^fi^ reaction involves a single photon, and thus it is more constrained 
(by factor 2) compared to the photon fusion, which is governed by two Compton type 
Feynman diagrams. However, in the rate we compute below, the indistinguishability of the 
two photons introduces an additional factor 1/2, so that both reactions differ only by the 
difference in the quantum Bose and Fermi distributions. 

Integrals in Eq. (]1.33p can be evaluated in spherical coordinates. We define: 

q = Pi + P2] p=^{pi-P2); q' = P3 + P4; p' = ^{ps - Pa); (6.274) 

z-axis is chosen in the direction oi ~q = pi + 

= (go, 0,0,0), p^ = {po,psme,0,pcos9), p'^ = (p[„ p' sin (/) sin sine/) cos cos <^). 
Now we obtain [40]: 



2 „s_ 
2 12 



1 (47r)(27r) f°° ry^i 



X j ^ d(cos 9) J ^ d(cos 0) ^ dxS (^p - [pI + '^'^ 6 (^pf - [pfl - + |) '^'^ 

^rf n qoPo\ r f , ^'oP'o\^\nj i2r^-2f fqo , \r f qo 

X6 {COSO - ^ 1 ^ ( COS0 - j 2^ |Me+e-(77)^^;.| + j - ^ 

xTf'et)/e(7) (f + P'o) fei,) (| " p'o) (qo/T), (6.275) 



where q* = gy 1 f-- The integration over p, p', cos 6', cos(f) can be done analytically 

considering the delta-functions. The other integrals can be evaluated numerically. For the 
case of indistinguishable colliding particles (two photons) there is additional factor 1/2 
implemented by the value / = 1, while for distinguishable colliding particles (here electron 
and positron) / = 0. 
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6.4.3 TT^ production 

Ti"^ can be produced in ttotto charge exchange scattering (11.701) and photons or e+e^ fusion, 
reactions ffLTTj) and (fr72|) 

We find that for vr^ production, the last two processes are much slower compared 
to the first, in case that ttq density is near chemical equilibrium upto temperatures ap- 
proximately 30 MeV. As we mentioned before rate of two production in photons fusion, 
Eq J1.73l is much smaller than rate of one 7r° production at considered temperatures. 

The time evolution equations for the number of vr^ are similar to Eq. (I6.267P : 

1 dN^± 



V dt 



+ (T^-T^±)i?,+e-^.+C6-2r6) 



For the respective three cross sections we use, all results valid in the common range 
s < 1 GeV^ we consider here: 

• The cross section for charge exchange vr^scattering reaction Eq. (ll.70p have been con- 
sidered in depth recently [93]: 



rr- '-^^l^ (0) (2) 2. .. 2771 



where a^^^ — Oo^'' = 0.27/Af^± This is the dominant process for charge pion production, 
subject to presence of 7r°. 



For process Eq. (1 1.71 1) , the cross section of vr^ production in photon fusion we use [94]: 
_ 27ra2 / AmlV^ ( m\ 

^77-+7r+7r~ I 

S \ S 

where my = 1400.0 MeV. As we will see from numerical calculations given the cross 



1 - ^ T-, TTT^, TT , (6-278) 

' ' ' (l/2s + m2,)(l/4s + m2.) ' ' ^ ' 



sections for 77 tt^tt resulting production rates will be smaller than the charge 
exchange vr^vr^ tt+tt" reaction. 

For process Eq. (ll.72|) . the cross section of vr^ production in electron - positron fusion 
we use [95]: 
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The form factor F[s) can be written in the form: 



- s + Tp{m'^p/k^)[k'^{h{s) - h{m?p)) + k'j^h\rn?p){rn?p - s)] - imp{k / kp)^T p{mp/ y/s) 

(6.280) 

where h'{s) = dh/ds and 

k=[-s-m^\ ■ kp=\-m-m^\ ; /i(s) = -^ln' 



mp = 775 MeV, T p = 130 MeV, d = 0.48. Given this cross section we also find that 
the rate of charged pion production is small when compared to TTo-charge exchange 
scattering. 

For reaction (ll.73p we have [96]: 

a^-y/s — 4m^ 



(t(77 n^n^ 
where 



87r2^ 



2 



l + —fs 

s 



a{7v+7v- tt'V), (6.281) 



2 

771 

/, = 2(ln2(z+/^.) - TT^) + -ZL(ln2(^+/^_) + vr^)^, (6.282) 



and z± = (1/2) (1 ± ^s-Aml). 

The cross sections for tt+tt^ pair production, evaluated using Eqs. fl6.277p . fl6.278p 
and (16.2791) are presented in figure 16.481 as functions of reaction energy ^/s . Top solid 
line (blue) is for charged pions production in 7r° scattering Eq. fll.70p . the magnitude of 
this cross section being very large we reduce it in presentation by factor 1000; the dashed 
line is for tt+tt^ production in photon fusion Eq. (11.711) : dash-doted line is for electron 
positron fusion Eq. (ll.72p . The bottom solid line (green) is for photon fusion into two 
neutral pions, Eq. (l6.28ip . The prediction for about 480 nb (data 420 nb) 

at the peak near threshold [96], which is in agreement with calculations presented here. 
The reaction (Eq. (11.731) ) is much smaller than others and we do not consider this 

reaction further. We note that some of these results are currently under intense theoretical 
discussion as they relate to chiral symmetry. For our purposes the level of precision of here 
presented reaction cross sections is quite adequate. 
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6.5 Numerical results 

6.5.1 Particle production relaxation times 

In figure 16.471 we show relaxation time r for the different processes considered as function 
of temperature T G [3, 50] MeV. Because of the large difference in production rates which 
can be compensated by different densities of particles present (magnitudes of fugacities) we 
introduce partial relaxation time for each of the three reactions vr^vr^ 7r"'"7r~, 77 — '7r"'"7r~ 
and e"*" + e~ tt+tt": 

T^o dn^± I c/T^± 1 dn^± / dT^± 1 dn^± / dT^± 

(6.283) 

When T <^ m, we can use the Boltzmann approximation to the particle distribution func- 
tions. Since in this limit the density is proportional to T the relaxation times doesn't 
depend on T. Moreover, even for T ^ 50 MeV, we have for muons e"™/^ ^ 1/3, thus 
quantum correlations in phase space remain small, and the Boltzmann limit can be em- 
ployed. To account for small deviation from Boltzmann limit arising towards the upper 
limit of the temperature range we consider, that is at T ~ 50 MeV, we used the exact equa- 
tions with Tj = 1 to calculate r for each case. In addition to these three cases Eq. fl6.283p 
we show in figure [633 the muon production relaxation time Eq. (16.268) ). the two photon 
fusion into vr" relaxation time Eq. (16.2641) . a nearly horizontal line (turquoise, bottom), 
which is slightly greater than the free space vr'^ decay rate. Finally, the thin dash-dot line 
at about 10^ times greater value of time is the electron-positron fusion into vr", Eq. (16.2611) . 

6.5.2 Rates of pion and muon formation 

In figure 16.491 we show on left solid (blue) line function of fireball temperature 
the rate per unit volume and time for the process 7 -|- 7 — 7r°, the dominant mechanism 
of pion production. The other solid line with dots corresponds to e+ -|- e~ ^ 7r° reaction 
which in essence remains, in comparison, insignificant. Its importance follows from the fact 
that it provides the second most dominant path to ttq formation at lowest temperatures 
considered, and it operates even if and when photons are not confined to remain in the 
plasma drop. 
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Figure 6.47: The relaxation time r for the different channels of pion and muon production (see 
box), as functions of plasma temperature T. 
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Figure 6.48: The cross section a for pion pair production, and pion charge exchange (sohd top 
hne) , as functions of ^/s < 1 GeV^ . 
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Figure 6.49: On left, the invariant pion production rates in units of nm~^fs~^, as a function of 
temperature T. On right the production rate R' per Joule energy content in the fireball, in units 
of MJ~^fs~^, in both cases for reactions shown in the box. 
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We improve the rate presentation on the right hand side in figure [639J considering 
that the formation of a plasma state involves an experimentally given fireball energy content 
£ in Joules, we use Eg. (16.2571) to eliminate the volume V at each temperature T: 



= duT" - gaT^^ dvVi'^ " ^^^^ ^^"^^^^ 

For chemical nonequilibrium, replace a —>■ T^cr(T). Considering the (good) approximate 
Eq. (l6.266p we obtain: 

We use units such that h = c = k = 1 and thus R' is a dimensionless expression. Recalling 
the value of these constants, the units we used for R' derive from MeV s=1.603 10~^ MJ 
is. 

The other lines in figure 16.491 address the sum of formation rates of charged pion 
pairs (dashed, red) by all reactions considered in this work, 7r° + 7r° — 7r+ + tt", 7 + 7 — > 
vr"*" + vr", e"*" + e~ ^ 7r+ + tt^. We also present the sum of all reactions leading to either 
a charged pion pair, or muon pair (dot-dashed, green) lines, adding in 7 + 7 + 
e"*" + e~ fi^ + The rationale for this presentation is that we do not care how a 
heavy particle is produced, as long as it can be observed. The dashed (red) line assumes 
that we specifically look for charged pions, and dot-dashed (green) line that we wait till 
charged pions decays, being interested in the total final muon yield. The 7r° production 
rate (blue, solid line) is calculated using Eq.f l4.4.3i) and yields on the logarithmic scale 
nearly indistinguishable result from the approximation Eq. fl6.266p . For vr^ production we 
refer to section 16.4.31 and for /i^ production we refer to I6.4.2[ 

In table 16.61 we show the values of key reaction rates R and relaxation times r at 
T = 5 and 15 MeV. We note the extraordinarily fast rise of the rates with temperature, in 
some instances bridging 15 - 20 orders in magnitude when results for T = 5 and 15 MeV 
are compared. 

In order to understand the individual contributions to the different reactions en- 
tering the sum of rates presented above, we show as a function of temperature in the figure 
16.501 the relative strength of muon pair (left) and charge pion (on right) electromagnetic 
(7 + 7, e"*" + e~) production, using as the reference the 7 + 7 — 7r° reaction. The /i^ 
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Table 6.6: Values of rates, relaxation times for all reactions at T = 5 MeV and T = 15 



MeV 





T = 5 MeV 


T = 5 MeV 


T = 15 MeV 


T = 15 MeV 


reaction 


r [as] 


i? [nm^'^fs-i] 


r [as] 


R [nm-3fs-i] 


77 ^ TTo 


8.82102 


3.3103 


9.5 102 


1.21012 




1.210^0 


3.2 10-3 


1.9103 


1.51011 


77 ^ jj,^ jj,' 


1.0 10^° 


3.710-3 


1.3103 


2.11011 




2.910^2 


2.1 10-^ 


4.6102 


9.5101° 


77 ^ Tl^TT^ 


6.410^3 


9.710-1° 


5.110^ 


8.7 10^ 


e~^e~ TT+TT- 


7.8 10^5 


7.910-12 


9.5 10^ 


4.6 10^ 



production rates are calculated using Eg. (16.2751) with |M|2 from Eg. (16.2711) and Eg. fl6.272p 
respectively. This ratio is smaller than unity for T ^ 20 MeV. For larger T, the muon 
direct production rate becomes larger than vr^ production rate. Charged pions (on right in 
figure I6.50p can be produced in direct reaction at a rate larger than neutral pions only for 
T > 35 MeV. The photon channel dominates. 

6.6 Discussion and Conclusions 

We found that the production of vr" is the dominant coupling of electromagnetic radiation to 
heavy (hadronic) particles with m 3> T, and as we have here demonstrated that noticeable 
particle yields can be expected already at modest temperatures T G [3, 10] MeV. In present 
day environment of 0.1 -1 J plasma lasting a few fs, our results suggest that we can 
expect integrated over space-time evolution of the EP3 fireball a 7r° yield at the limit of 
detectability. For T 15 MeV the 7r° production rate remains dominant and indeed 
very large, reaching the production rate R' ~ 10i^[MJ-ifs-i]. Charge exchange reactions 
convert some of the neutral pions into charged pions which are more easy to detect. 

In this situation it is realistic to consider the possibility of forming a chemically 
eguilibrated fireball with 7r°, vr^, fi^ in chemical abundance eguilibrium. The heavy par- 
ticles are produced in early stages when temperature reached is highest. Their abundance 
in the fireball follows the fireball expansion and cooling till their freeze-out, that is de- 
coupling of population eguation production rates. The particle yields are than given by 
the freeze-out conditions, specifically the chemical freeze-out temperature Tf and volume 
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T [MeV] T [MeV] 

Figure 6.50: On left: Muon and on right charged pion production rates in electromagnetic 
processes normalized by vr'^ production rate. Solid line (blue) for 77 , dashed line (green) for 
e"'"e~ induced process. 
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Vf, rather than the integral over the rate of production. In this situation the heavy par- 
ticle yields become diagnostic tools of the freeze-out conditions, with the mechanisms of 
their formation being less accessible. However, one can avoid this condition by appropriate 
staging of fireball properties. 

The present study has not covered, especially for low temperature range all the 
possible mechanisms, and we addressed some of these issues in the introduction. Here we 
note further that the production of heavy particles requires energies of the magnitude m/2 
and thus is due to collisions involving the (relatively speaking) far tails of a thermal particle 
distribution. If these tails fall off as a power law, instead of the Boltzmann exponential 
decay [97], a much greater yield of heavy particles could ensue. There could further be 
present a collective amplification to the production process e.g. by residual matter flows, 
capable to enhance the low temperature yields, or by collective plasma oscillations and 
inhomogeneities. 

These are just some examples of many reasons to hope and expect a greater particle 
yield than we computed here in microscopic and controllable two particle reaction approach. 
This consideration, and our encouraging 'conventional' results suggest that the study of 7r° 
formation in QED plasma is of considerable intrinsic interest. Our results provide a lower 
limit for rate of particle production and when folded with models of EP^ fireball formation 
and evolution, final yield. 

It is of some interest to note that the study of pions in QED plasma allows explo- 
ration of pion properties in electromagnetic medium. Specifically, recall that 1.2% fraction 
of TT^ e"'"e~7 decays, which implies that the associated processes such as e^ + e~ —^'j + 'iiq 
are important. We cannot evaluate this process at present as it involves significant chal- 
lenges in understanding of ttq off-mass shell 'anomalous' couping to two photons. 

The experimental environment we considered here should allow a detailed study 
of the properties of pions (and also muons) in a thermal background. There is considerable 
fundamental interest in the study of pion properties and specifically pion mass splitting 
in QED plasma at temperature T > Am and in presence of electromagnetic fields. We 
already have shown that due to quantum statistics effects, the effective in medium decay 
width of 7r° differs from the free space value, see figure 14.311 In addition, modification of 
mass and decay width due to ambient medium influence on the pion internal structure is 
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to be expected. Further we hope that the study of pions in the EP^ fireball will contribute 
to the better understanding of the relatively large difference in mass between vr^and vr^. 
The relatively large size of the PE^ environment should make such changes, albeit small, 
measurable. 

The experimental study of 7r° in QED plasma environment is not an easy task. 
Normally, one would think that the study of the decay into two 67.5 MeV 7 (+ thermal 
Doppler shift motion) produces a characteristic signature. However, the vr" decay is in time 
and also in location overlapping with the plasma formation and disintegration. The debris 
of the plasma, reaches any detection system at practically the same time instance as does 
the 67.5 MeV 7. The large amount of available radiation will disable the detectors. On 
the other hand we realize that the hard thermal component of the plasma, which leads to 
the production of 7r° in the early fireball stage, is most attenuated by plasma dynamical 
expansion. Thus it seems possible to plan for the detection of 7r° e.g. in a heavily shielded 
detection system. 

The decay time of charged pions being 26 ns, and that of charged muons being 
2.2 fis it is possible to separate in time the plasma debris from the decay signal of these 
particles. Clearly, these heavy charged particles can be detected with much greater ease, 
also considering that the decay product of interest is charged. For this reason, we also have 
in depth considered all channels of production of charged pions and muons. Noting that 
practically all charged pions turn into muons, we have also compared the production rates 
of 7r° with all heavy particles, see dot-dashed (green) line in figure 16.501 This comparison 
suggests that for plasmas at a temperature reaching T > 10 MeV the production of final 
state muons will most probably be by far easier to detect. On the other hand for T < 
5 MeV it would seem that the yield difference in favor of 7r° outweighs the detection 
system/efficiency loss considerations. Future work addressing non-conventional processes 
will show at how low T we can still expect observable heavy particle yields. 

An effort to detect 7r° directly is justified since we can learn about the properties 
of the plasma (lifespan, volume and temperature in early stages) e.g. from a comparative 
study of the and vr^ production. We have found that at about T > 16 MeV, the pion 
charge exchange tt^tt^ tt+tt^ reaction for chemically equilibrated 7r° yield is faster than 
the natural decay, and the chemical equilibration time constant, see the dot-dashed line 
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in figure I6.47[ Thus beyond this temperature the yield of charged pions can be expected 
to be in/near chemical equilibrium for a plasma which lives at, or above this temperature, 
for longer than 100 as. 

In such an environment the yield of vr'' is expected to be near chemical equilibrium, 
since the decay rate is compensated by the production rate, and, within 100 as, the chemical 
equilibrium yield is attained. Moreover, the thermal speed of produced tt can be obtained 
from the nonrelativistic relation |m(f^) = |T, thus v oc VT and, for T = 10 MeV, v ~ 
0.5c. This is nearly equal to the sound velocity of EP'^, Vg — c/ Vs = 0.58c. Thus the heavy 
tt" particles can be seen as co-moving with the expanding/exploding EP^, which completes 
the argument to justify their transient chemical equilibrium yield in this condition. 

The global production yield of neutral and charged pions should thus allow the 
study of volume and temperature history of the QED plasma. More specifically, since with 
decreasing temperature, for T < 16 MeV, there is a rapid increase of the relaxation time 
for the charge exchange process, there is a rather rapid drop of the charged pion yield 
below chemical equilibrium — we note that charge exchange equilibration time at T = 10 
MeV is a factor 10^ longer. We note that the study of two pion correlations provides an 
independent measure of the source properties (HBT measurement). 

The relaxation time of electromagnetic production of muon pairs wins over vr" 
relaxation time for T > 22 MeV, see dashed line, red, in figure \6A7\ the direct electro- 
magnetic processes of charged pion production (thin green, solid line for 77 tt+tt^ and 
dashed, blue for e~^e~ — > tt+tt^) remain sub-dominant. Thus for T > 22 MeV we expect, 
following the same chain of arguments for muons as above for charged pions, a near chem- 
ical equilibrium yield. If the study of all these 7r'',7r^,/i^ yields, their spectra and even 
pion correlations were possible, considerable insight into e~,e"'",7 plasma (EP^) plasma 
formation and dynamics at T < 25 MeV can be achieved. 



176 



CHAPTER 7 
PION AND MUON IN EARLY UNIVERSE 



In this chapter we begin to apply methods considere in previous chapters to early Universe. 
These all reactions of muon and pions production, considered in e"'"e~7 plasma take place 
in early universe. Here we show that vr*^ is in chemical equilibrium with photons at all 
temperatures of interest. 

In expanding universe in metric [75] 

ds" = de - R\t) ( + r^de^ + r^sin'edclA , (7.286) 

\1 — kr'^ J 

a =- (7.287) 

1 — fcr^ 

The Eq. (I4.113p . which describes ttq evolution, has dilution term: 

+ mn^ = (1 - T^)^, (7.288) 

dt 

where A is defined by Eq. (14.1261) 

driTT driTT dn.,, ■ 

Then dividing both sides of equation (I7.288P by dn^^/dX and using Eq. fl4.125p we obtain 

- -T. + 3^ 7^(1: = (1 - T.)-, (7.290) 
Tt driT,/ dlj, r 

where 

^ = (7.291) 

and 

1 drij^/dT T 



(7.292) 



Tt driT^/dTT^Tj, 
We put '-' sign to this equation to have tt > 0. 

The temperature can be defined from entropy conservation for radiation dominated 

epoch we have 

f R , , 
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Now we estimate how large is time scale I/tj-. For uIt^/T « 1 the n^^ oc and 
using ( ]7.293p we have in chemical equilibrium (T^o = 1): 

— = 3H. (7.294) 

Tt 

Therefore for ultrarelativistic particles dilution rate compensate the rate of density 
decrease with temperature change in expanding universe. So density of pions may stay as in 
chemical equlibrium for all temperatures while pions are ultrarelativistic. Dilution doesn't 
have direct effect on solution. However, when interaction rate becomes small compare to 
the expansion rate, if pions aren't equilibrium with photons for some reasons it takes large 
time compared to universe age at that moment to get to equilibrium density. So they 
become decoupled from radiation. This decoupling takes place when 

?,H > -. (7.295) 

T 

The value of R/R can be find from Friedmann equation [75] 
For radiation dominated epoch, if /c = 1 we have 

^'■''''> 

where g* is number of degrees of freedom. Then 

H = l.eeVr (7.298) 

rripi 

If we assume t ^ tq = 8.4 10"^'' s we have from (17.295^ condition for freeze-out temperature 

T ^ lO^GeV (7.299) 

For T < 10^ GeV the expansion rate is small compare to vr^ production and decay rates. 
If m/T > 1 for Boltzmann distribution 



1 / TCTTI'^ T"""^ 

= TT^^y ^—^M-mJT) (7.300) 

ZTT^ V 2 



we have 

1 
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For small T « we have 1/tx » H. If l/r^ exceeds decay rate the pions 
may lose chemical equilibrium for small T when universe is matter dominated. For matter 
dominated universe 

Here we used equation (17.2931) . From Eq. (l7.296p we have 



T 



\ 3/2 



H = Ho./no\^-j (7.303) 
- = Ho^/no^^^^, (7.304) 



where 



and To ~ 10^^^ GeV. 1/tt is decreasing as ^/T and it can't reach value of pions decay 
width Fo, 



-p rp3/2 , piXg 

Therefore for considered temperature range the vr" is in chemical equilibrium with 
photons because of their fast decay rate. This is not always that decay is so fast to exceed 
universe expansion rate. For example the decay n — > p + e~ + z/e is much slower r = 885.7 
s. The dilution rate exceeds neutron decay rate at T > 0.1 MeV. 

The relaxation time for fi^ and vr^ in reactions (11.741) and (11.701) respectively 
become many orders of magnitude larger than r^o at temperatures about few MeV, where 
these reaction have to freeze out. Therefore these particles are in chemical equilibrium and 
their densities are also relatively high (about nucleons density or higher) up to temperatures 
of few MeV. 

This process is important to understand how the hadronic component diminish 
with the expansion of the Universe and the possible effects of hadronic relics in the cosmic 
blackbody radiation spectrum, such as its fluctuations and correlations. 
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CHAPTER 8 
SUMMARY AND CONCLUSIONS 

In the first part of dissertation we studied lieavy particles production at liadronization, 
resonance evolution in thermal hadronc gas after hadronization. In the second part e"^e~7 
plasma equilibrium conditions were considered and pion and muon production in this type 
of plasma. 

8.1 Summary of heavy flavor production 

In chapter [3] I considered heavy flavor (charm, bottom) hadrons production within sta- 
tistical hadronization model. The new feature compared to the others studies is that we 
assume entropy and strangeness conservation during hadronization. 

While I compare the yields to the expectations based on chemical equilibrium 
yields of light and strange quark pairs, I present results based on the hypothesis that the 
QGP entropy and QGP flavor yields determine the values of phase space occupancy 
i = q,s,c,b, which are of direct interest in study of the heavy hadron yields. 

For highest energy heavy ion collisions the range of values discussed in literature 
is 1 < 7^ < 1.65 and 0.7 < 7^/7^ < 1.5. However 7^ and 7^ values which are much 
larger than unity arise. This is due to the need to describe the large primary parton based 
production, and considering that the chemical equilibrium yields are suppressed by the 
factor exp(— m/T). 

Our work is based on the grand canonical treatment of phase space. This approach 
is valid for charm hadron production at LHC, since the canonical corrections, as we have 
discussed, are not material. On the other hand, even at LHC the much smaller yields of 
bottom heavy hadrons are subject to canonical suppression. The value of the parameter 
7^ obtained at a fixed bottom yield Nj,, using either the canonical, or the grand canonical 
methods, are different, see e.g. Eq. (15) in [71]. Namely, to obtain a given yield A^^^ in 
canonical approach, a greater value of 7^ is needed in order to compensate the canonical 
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suppression effect. However, for any individual single-6 liadron, the relative yields, i.g. 
B/Bs do not depend on 7^ and thus such ratios are not influenced by canonical phase 
spase effect. Moreover, as long as the yield of single-6 hadrons dominates the total bottom 
yield: A^^t ~ i? + i?s + A;, + . . ., also the A*";, scaled yields of hadrons comprising one 6-quark 
i.e. ratios such as B/Nf,, B^/Ni,, Bc/Nh, etc, are not sensitive to the value of 7^ and can 
be obtained within either the canonical, or grand canonical method. On the other hand 
for bb mesons and multi-6 baryons the canonical effects should be considered. Study of the 
yields of these particles is thus postponed. 

I have addressed here how the yields of heavy hadrons are influenced by 7^/7^ 7^ 1 
and 7g 7^ 1. The actual values of 7^/7^ we use are related to the strangeness per entropy 
yield s/S established in the QGP phase. Because the final value s/S is established well 
before hadronization, and the properties of the hadron phase space are well understood, 
the resulting 7^/7^ are well defined and turn out to be quite different from unity in the 
range of temperatures in which we expect particle freeze-out to occur. We consider in some 
detail the effect of QGP hadronization on the values of 7^ and 7^. 

One of first results I present (figure I3.13P allows a test of the statistical hadroniza- 
tion model for heavy flavor: I show that the yield ratio cc ss/ {cs cs) is nearly independent 
of temperature and it is also nearly constant when the is allowed to freeze-out later 
(figure [3TT^ . provided that the condition of production is at the same value of strangeness 
per entropy s/S. 

I studied in depth how the (relative) yields of strange and non-strange charmed 
mesons vary with strangeness content. For a chemically equilibrated QGP source, there 
is considerable shift of the yield from non-strange D to the strange Dg for s/S = 0.04 
expected at LHC. The expected fractional yield Dg/N^ — Bg/N^ ~ 0.2 when one assumes 
7? = = expected enhancement of the strange heavy mesons is at the level of 

30% when s/S = 0.04, and greater when greater strangeness yield is available. 

As the result we find a relative suppression of the multi-heavy hadrons, except 
when they contain strangeness. This suppression depends on both factors 7s and gammaq. 
When phase space occupancy of light and strange quark is relatively high the probability 
for charm quarks to make hadrons with strange quarks increases and probability to find 
the second charm quark among light and strange quarks decreases. Therefore the cc yield 
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suppression increases when '~is/lq ratio increases for constant 7^. This result is quahtatively 
in agreement with experimental results obtained for SPS energies [36]. 

On the other hand, the yield of cc/N'^ ^ 210^'^ is found to be almost independent 
on hadronization temperature when entropy at hadronization is conserved. That is be- 
cause for larger T 7^ decreases. The suppression effect decreasees, compared to SHM and 
become even negative for T > 200 MeV, resulting to the cc yield almost independent on 
temperature. We don't know exactly equation of state in QGP and so the value of 7^ which 
is needed to conserve the entropy may be different. If 7^ is larger for higher temperatures, 
suppression of cc is larger for a fixed s/ S. The same result is found for i?c ^ 5 — 6 lO'^A'^cA'";,, 
that yield remains considerably larger (by a factor 10 — 100) compared to the scaled yield 
in single nucleon nucleon collisions. 

I have shown that the study of heavy fiavor hadrons will provide important infor- 
mation about the nature and properties of the QGP hadronization. The yield of Bc(6c) 
mesons remains enhanced while the hidden charm cc states encounter another suppression 
mechanism, compensating for the greatly enhanced production due to large charm yield at 
LHC. The results are published in [1] 

8.2 Summary on Chemical Equilibration Involving Decaying Particles at Fi- 
nite Temperature 

In chapter HlI examined in detail the kinetic master equation for the process involving 
formation of an unstable particle through the reaction Eq. (11.481) in a relativistically co- 
variant fashion. Assuming that all the particles in the process are in thermal equilibrium, 
we calculate the thermal averaged decay and formation rate of the unstable particle based 
on the BUU equation. Using the time reversal symmetry, we show that the time evolution 
of the density of the unstable particle as Eq. (14.1081) . Therefore in chemical equilibrium 
particles fugacities are connected by Eq. (I5.16ip as expected. We have explicit the thermal 
decay rate of unstable particle, obtaining Eq. (l4.140p . which is our principal result. 

Using the formalism developed above, I examined the general properties of the 
thermal particle decay/production rate. We see that for T <^ rrii where the Boltzmann 
limit can be applied, the decay width is reduced to Ti/tq and production width is T2T3/ro. 
For larger values of T but Tj ^ 1 so that the Boltzmann approximation is valid, then 
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decay width and production width tend simply to Ti/r and respectively, where 

r is essentially proportional to average Lorentz factor and doesn't depend on Tj. When 
some of rrii/T and Tj are about unity or larger we see dependence of r on Tj. 

I applied our formalism to 3 examples, p *^ tt + tt, S(1385) ^ and 7r° 7 + 7. The 
first and second processes can take place both in a hot hadronic gas created by the heavy 
ion collisions and in the expanding early Universe. In particular for the heavy ion reaction 
case, our analysis, coupled to the hydrodynamical expansion of the system will furnish 
additional information of the dynamics of the system. We will study baryon resonances 
evolution in heavy ions collisions in next chapter. The relaxation time for 7r° decay remains 
close (within 50%) to relaxation time in vacuum for large temperature range. In chapter [6] 
we will apply this for 7r° evolution in e~''e~7 plasma, created by the intensive laser pulse 
and in early universe. This part is going to be publised in [3]. 

8.3 Summary of resonance production in heavy ions collisions 

In chapter [5] I apply equations derived in chapter H] to baryon resonance densities evolution 
in thermal hadron gas after QGP hadronization. The goal is to explain ratios S(1385)/A° 
and A(1520)/A° reported by RHIC experiment and also predict A(1232)/A^ ratio. 

The resonant hadron states, considering their very large decay and reaction rates, 
can interact beyond the chemical and thermal freeze-out of stable particles. Thus the 
observed yield of resonances is fixed by the physical conditions prevailing at a later breakup 
of the fireball matter rather than the production of non-resonantly interacting hadrons. 
Moreover, resonances, observed in terms of the invariant mass signature, are only visible 
when emerging from a more dilute hadron system given the ample potential for rescattering 
of decay products. The combination of experimental invariant mass method with a large 
resonant scattering makes the here presented population study of resonance kinetic freeze- 
out necessary. The evolution effects we find are greatly amplified at low hadronization 
temperatures where greatest degree of initial chemical equilibrium is present. 

Our study quantifies the expectation that in a dense hadron medium narrow res- 
onances are "quenched" [81] that is, effectively mixed with other states, and thus their 
observed population is reduced. Since we follow here the particle density, the effect we 
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study is due to incoherent population mixing of A(1520), in particular with E*. This effect 
is possible for particle densities out of chemical non-equilibrium. However, this mixing can 
occur also at the amplitude (quantum coherent) level. As the result the yield suppres- 
sion effect could further increase, in some situations further improving the agreement with 
experiment. 

Our results show that the observable ratio A(1520)ob/Atot can be suppressed by 
two effects. First A(1520) yield is suppressed due to excitation of heavy S*s in the res- 
onance scattering process. Moreover, the final A(1520)ob yield is suppressed, because 
S*s, which decay to A(1520), are suppressed at the end of the kinetic phase evolution 
by their (asymmetric) decays to lower mass hadrons, especially when dead channels are 
present (see figure [530]) . As a result, fewer of these hadrons can decay to A(1520)ob during 
the following free expansion. A contrary mechanism operates for the resonances such as 
S(1385), A(1230). These resonances can be so strongly enhanced, that in essence most 
final states strange and non-strange baryons come from a resonance decay. 

We note that despite a scenario dependent resonance formation or suppression, 
the stable particle yields used in study of chemical freeze-out remain unchanged, since all 
resonances ultimately decay into the lowest "stable" hadron. Therefore after a description 
e.g. within a statistical hadronization model of the yields of stable hadrons, the under- 
standing of resonance yields is a second, and separate task which helps to establish the 
consistency of our physical understanding of the hadron production process. 

We conclude noting the key result of this study, that we can now understand 
the opposite behavior of A(1520) (suppression in high centrality reactions) and S*^1385) 
(enhancement, and similarly A(1230)) by considering their rescattering in matter. In order 
to explain both, the behavior of the A(1520)ob/Atot and S(1385)/Atot ratios, one has to 
consider T = 95 — 100 MeV as the favorite temperature of final kinetic freeze-out of hadron 
resonances, with Tq = 140 MeV being the favored chemical freeze-out (hadronization, QGP 
break-up) temperature. When there is little matter available to scatter, e.g. in peripheral 
collisions, the average value of A(1520)ob/Atot ratio is higher, approaching the expected 
chemical freeze-out hadronization yield for Tq = 140 MeV. All these findings are in good 
agreement with available experimental data. 

This part is published in [4] and [5]. 
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8.4 Results for relatvistic e+e 7 plasma created by laser pulse 

In chapter [6] I consider e"'"e~7 plasma. We study the freeze-out condition for a relativistic 
e~e"'"7 plasma, when plasma density becomes so low that particles begin to stream freely. 

In order to address this question we evaluated Lorentz invariant rates for the 
Compton scattering and pair production in thermally and chemically equilibrated EP^ 
plasma. We then used these Lorentz invariant rates to evaluate the corresponding mean 
free path length / of particles. 

Comparing I to plasma drop size we showed that an opaque equilibrium density 
plasma drop can be produced at energy 0.5 kJ in the volume with largest possible radius 
R = 7 nm. This volume corresponds to the smallest possible temperature T = 2 MeV. 
In order to reach higher than 2 MeV temperature, we need to increase energy of plasma 
(which is proportional to T^) or/and decrease plasma size. At higher temperature opaque 
plasma can be created at the total plasma energy smaller than 0.5 MeV, since smaller 
plasma drop size is in agreement with the opaqueness condition Eq. (I6.227p . as seen in 
figure I6.44[ 

On the other hand in order to create opaque plasma with temperature lower than 
2 MeV, the necessary amount of energy is larger than 0.5 kJ. This is so because the plasma 
size has to be large in order to satisfy opaqueness condition Eq. (16.2271) . 

Our main result, perhaps unexpected at the first sight is illustrated in figure [6311 
For the temperature T > 2 MeV equilibrium plasma production with a relatively small en- 
ergy pulse (compared to lower temperature equilibrium plasma) may be possible. However, 
the challenge here is to focus the energy into the volume of size < 10 nm. 

These results are in preparation [6]. 

I also study heavy particles (pion, muon) production in e~^e~'~f plasma. We found 
that the production of 7r° is the dominant coupling of electromagnetic radiation to heavy 
(hadronic) particles with m ^ T, and as we have here demonstrated that noticeable 
particle yields can be expected already at modest temperatures T G [3, 10] MeV. In present 
day environment of 0.1 -1 J plasma lasting a few fs, our results suggest that we can 
expect integrated over space-time evolution of the EP^ fireball a vr*^ yield at the limit of 
detectability. For T — > 15 MeV the 7r° production rate remains dominant and indeed 
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very large, reaching the production rate E! ~ lO^^fMJ-^fs"^]. Charge exchange reactions 
convert some of the neutral pions into charged pions which are more easy to detect. 

In this situation it is realistic to consider the possibility of forming a chemically 
equilibrated fireball with tt*^, tt^, in chemical abundance equilibrium. The heavy par- 
ticles are produced in early stages when temperature reached is highest. Their abundance 
in the fireball follows the fireball expansion and cooling till their freeze-out, that is de- 
coupling of population equation production rates. The particle yields are than given by 
the freeze-out conditions, specifically the chemical freeze-out temperature and volume 
V/, rather than the integral over the rate of production. In this situation the heavy par- 
ticle yields become diagnostic tools of the freeze-out conditions, with the mechanisms of 
their formation being less accessible. However, one can avoid this condition by appropriate 
staging of fireball properties. 

The present study has not covered, especially for low temperature range all the 
possible mechanisms, and we addressed some of these issues in the introduction. Here we 
note further that the production of heavy particles requires energies of the magnitude m/2 
and thus is due to collisions involving the (relatively speaking) far tails of a thermal particle 
distribution. If these tails fall off as a power law, instead of the Boltzmann exponential 
decay [97], a much greater yield of heavy particles could ensue. There could further be 
present a collective amplification to the production process e.g. by residual matter fiows, 
capable to enhance the low temperature yields, or by collective plasma oscillations and 
inhomogeneities. 

These are just some examples of many reasons to hope and expect a greater particle 
yield than we computed here in microscopic and controllable two particle reaction approach. 
This consideration, and our encouraging 'conventional' results suggest that the study of 7r° 
formation in QED plasma is of considerable intrinsic interest. Our results provide a lower 
hmit for rate of particle production and when folded with models of EP^ fireball formation 
and evolution, final yield. 

It is of some interest to note that the study of pions in QED plasma allows explo- 
ration of pion properties in electromagnetic medium. Specifically, recall that 1.2% fraction 
of 7r° e"'"e~7 decays, which implies that the associated processes such as e'^-\-e~ — > 7+7ro 
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are important. We cannot evaluate this process at present as it involves significant chal- 
lenges in understanding of ttq off-mass shell 'anomalous' couping to two photons. 

The experimental environment we considered here should allow a detailed study 
of the properties of pions (and also muons) in a thermal background. There is considerable 
fundamental interest in the study of pion properties and specifically pion mass splitting 
in QED plasma at temperature T > Am and in presence of electromagnetic fields. We 
already have shown that due to quantum statistics effects, the effective in medium decay 
width of tt" differs from the free space value, see figure 14.311 In addition, modification of 
mass and decay width due to ambient medium influence on the pion internal structure is 
to be expected. Further we hope that the study of pions in the EP^ fireball will contribute 
to the better understanding of the relatively large difference in mass between vr^and vr^. 
The relatively large size of the PE^ environment should make such changes, albeit small, 
measurable. 

The experimental study of in QED plasma environment is not an easy task. 
Normally, one would think that the study of the 7r° decay into two 67.5 MeV 7 (+ thermal 
Doppler shift motion) produces a characteristic signature. However, the vr" decay is in time 
and also in location overlapping with the plasma formation and disintegration. The debris 
of the plasma, reaches any detection system at practically the same time instance as does 
the 67.5 MeV 7. The large amount of available radiation will disable the detectors. On 
the other hand we realize that the hard thermal component of the plasma, which leads to 
the production of 7r° in the early fireball stage, is most attenuated by plasma dynamical 
expansion. Thus it seems possible to plan for the detection of vr" e.g. in a heavily shielded 
detection system. 

The decay time of charged pions being 26 ns, and that of charged muons being 
2.2 /is it is possible to separate in time the plasma debris from the decay signal of these 
particles. Clearly, these heavy charged particles can be detected with much greater ease, 
also considering that the decay product of interest is charged. For this reason, we also have 
in depth considered all channels of production of charged pions and muons. Noting that 
practically all charged pions turn into muons, we have also compared the production rates 
of 77° with all heavy particles, see dot-dashed (green) line in figure 16.501 This comparison 
suggests that for plasmas at a temperature reaching T > 10 MeV the production of final 
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state muons will most probably be by far easier to detect. On the other hand for T < 
5 MeV it would seem that the yield difference in favor of vr*^ outweighs the detection 
system/efficiency loss considerations. Future work addressing non-conventional processes 
will show at how low T we can still expect observable heavy particle yields. 

An effort to detect 7r° directly is justified since we can learn about the properties 
of the plasma (lifespan, volume and temperature in early stages) e.g. from a comparative 
study of the vr" and vr^ production. We have found that at about T > 16 MeV, the pion 
charge exchange tt^ti^ —>■ tt+tt^ reaction for chemically equilibrated 7r° yield is faster than 
the natural 7r° decay, and the chemical equilibration time constant, see the dot-dashed line 
in figure 16.471 Thus beyond this temperature the yield of charged pions can be expected 
to be in/near chemical equilibrium for a plasma which lives at, or above this temperature, 
for longer than 100 as. 

In such an environment the yield of vr" is expected to be near chemical equilibrium, 
since the decay rate is compensated by the production rate, and, within 100 as, the chemical 
equilibrium yield is attained. Moreover, the thermal speed of produced tt can be obtained 
from the nonrelativistic relation |m(i;^) = |T, thus v oc VT and, for T = 10 MeV, v ~ 
0.5c. This is nearly equal to the sound velocity of EP'^, Vs — c/ = 0.58c. Thus the heavy 
tt" particles can be seen as co-moving with the expanding/exploding EP^, which completes 
the argument to justify their transient chemical equilibrium yield in this condition. 

The global production yield of neutral and charged pions should thus allow the 
study of volume and temperature history of the QED plasma. More specifically, since with 
decreasing temperature, for T < 16 MeV, there is a rapid increase of the relaxation time 
for the charge exchange process, there is a rather rapid drop of the charged pion yield 
below chemical equilibrium — we note that charge exchange equilibration time at T = 10 
MeV is a factor 10^ longer. We note that the study of two pion correlations provides an 
independent measure of the source properties (HBT measurement). 

The relaxation time of electromagnetic production of muon pairs wins over 7r° 
relaxation time for T > 22 MeV, see dashed line, red, in figure I6.47[ the direct electro- 
magnetic processes of charged pion production (thin green, solid line for 77 tt+tt" and 
dashed, blue for e~^e~ tt+tt^) remain sub-dominant. Thus for T > 22 MeV we expect. 
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following the same chain of arguments for muons as above for charged pions, a near chem- 
ical equilibrium yield. If the study of all these , n"^ , fi"^ yields, their spectra and even 
pion correlations were possible, considerable insight into e~,e"'",7 plasma (EP^) plasma 
formation and dynamics at T < 25 MeV can be achieved. This part is publised in [7] 

In chapter [7| we studied the pion equilibration in early universe. In early universe 
for temperature range of interest the 7r° are in chemical equilibrium with photons because 
of their fast decay rate. This is not always that decay is so fast to exceed universe expansion 
rate. For example the decay n p + e~ + Uf. is much slower r = 885.7 s. The dilution rate 
exceeds neutron decay rate at T > 0.1 MeV. 

The relaxation time for fi^ and vr^ in reactions (11.741) and (11.701) respectively 
become many orders of magnitude larger than r^o at temperatures about few MeV, where 
these reaction have to freeze out. Therefore these particles are in chemical equilibrium and 
their density is also relatively high upto temperatures of few MeV. 

These processes are important to understand how the hadronic component dimin- 
ish with the expansion of the Universe and the possible effects of hadronic relics in the 
cosmic blackbody radiation spectrum, such as its fluctuations and correlations. This part 
is going to be published in [3]. 
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